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Kinetisches Durchschlagen des schwach gekriimmten Stabes 
Von E. Mettler und F. Weidenhammer 


Herrn K. von Sanden zum 75. Geburtstag gewidmet 


1. Fragestellung. Zwischen den Stabilitatsproblemen der Elastostatik und der Elastokinetik 
besteht eine enge Analogie. In der Statik ist es bekanntlich méglich, da8 ein ,,diinner“ elastischer 
Korper (Stab, Platte, Schale usw.) einer statischen Last von bestimmter Hohe, die ein indifferentes 
Gleichgewicht bezeichnet, auszuweichen beginnt und bei weiterer Steigerung der Last die alte, 
instabil gewordene Lage verlaBt, um eine neue, ,,iiberkritische Gleichgewichtslage einzunehmen. 
Es entsteht dann also Mehrdeutigkeit des Gleichgewichts. Man unterscheidet zwei Typen: von 
solchen Ausweichproblemen!, nimlich Verzweigungsprobleme und Durchschlagprobleme. Diese 
Phanomene erhalt man allerdings nur, wenn man die streng linearisierte Elastizitatstheorie verlaBt 
und nichtlineare Glieder beriicksichtigt. Die entsprechenden Stabilitatsprobleme der Elasto- 
kinetik entstehen einfach dadurch, da8 man an die Stelle der statischen Lasten zeitlich verander- 
liche Lasten setzt. Dabei beschrankte man sich bisher ganz iiberwiegend auf harmonisch pulsierende 
Lasten. Statt mit statischen Gleichgewichtslagen hat man es jetzt mit erzwungenen Schwingungen 
des elastischen Kérpers zu tun. Indifferenzpunkte und iberkritische Zustinde treten analog zu 
den statischen Erscheinungen auf, wobei natiirlich die Verhaltnisse komplizierter werden, weil 
die Zeit als weitere Variable hinzutritt. Insbesondere spielen die iiberkritischen Zustande als nicht- 

_lineare Schwingungen im kinetischen Fall eine gréBere Rolle als in der Statik. 

Kinetische Verzweigungsprobleme findet man in der Literatur in weitem Umfang behandelt; 
man denke nur an das wohlbekannte Beispiel des axial pulsierend belasteten Stabes, der das Ana- 
logon zum Eulerschen Knickstab bildet. Das kinetische Durchschlagproblem ist daneben nur 

ve wenig bearbeitet worden. Vor kurzem hat der eine von uns” auf die in Rede stehende Analogie 


der Stabilitatsprobleme allgemein hingewiesen, konnte aber das Durchschlagproblem nur am a 
Beispiel des schwach gekriimmten Stabes mit einer sehr speziellen Belastung, namlich einer sinus- ee 
formig verteilten Streckenlast, behandeln. Es diirfte deshalb von Interesse sein, dieses Problem pe 

a erneut aufzugreifen und fiir allgemeinere, in der Praxis vorkommende Belastungen zu lésen. a 
2. Ausgangsgleichungen. Wir betrachten einen an beiden Enden in festen Gelenken gelagerten ee 

- Stab nach Abb. 1 mit der kleinen spannungslosen und zeitlich konstanten Voraushiegung W(x) S 
: a 
. 3 

mw p(at) 3 

om 

Z | : y 

Abb. 1. Der Stab im Koordinatensystem. e. 

Die Streckenlast p(x,t) ist der Deutlichkeit halber nach oben verschoben gezeichnet. a 


und der zusitzlichen, von der Zeit t abhangigen kinetischen Querverschiebung w(x, t). Die Lings- 
_-verschiebung sei u(x,t). Der Stab trage die orts- und zeitabhangige Querlast je Langeneinheit 


‘ Vel. C. B. Biezeno und R. Grammel, Technische Dynamik Band 1, S. 562, 2. Aufl., Berlin/Géttingen/ 2 


- Heidelberg 1953. 
L 2 EF. Matter: 1X¢ Congr. Intern. de Mec. appl., Actes t. V., p. 5, Bruxelles 1957. 
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p(x, t) (Abb. 1). Seine Bewegungsgleichungen und Endbedingungen lassen sich zweckmaBig aus 
dem Prinzip von Hamilton 


5 f(T—V;—V,) dt + [6A dt =0 (1) 


herleiten. Darin bezeichnet V, die potentielle Energie der Formanderung, V, das Potential der 
4uBeren Belastung, T die Bewegungsenergie und 6A die virtuelle Arbeit der Dampfungskraft. 
Die Formanderungsenergie berechnet sich nach K. Marguerre? und R. Kappus? in einer Naherung, 
welche die ersten nichtlinearen Glieder beriicksichtigt, mit Hilfe des Ausdruckes fiir die Dehnung 


der Stabachse 


x w2 
é=u, ++, W, (2) 
zu 
2 2 EJ 
y= (4 +3 +.) de +7 | whe de. ~ (3) 
0 0 


Fiir die weiteren Energieanteile hat man in (1) die aus der Schwingungslehre gelaufigen Ausdriicke 


l 


1 I 
V 33>, [ p(x, t) w dx, T= 5 | (ul +ut) de. 64 = — | Bu, du dx (4) 
J 
6 6 


einzusetzen. Die Arbeit der als geschwindigkeitsproportional vorausgesetzten Dampfungskrafte 
wird nur bei der wesentlich interessierenden Querverschiebung w(x, t) beriicksichtigt. In (3) und (4) 
bezeichnen J das fiir Biegeschwingungen maBgebende, konstante Flachentragheitsmoment, F die 
konstante Querschnittsflache des Stabes, w die Stabmasse je Langeneinheit, | die Stablange, E den 
Elastizitatsmodul und 6 die Dampfungskonstante der Biegeschwingungen. 

Durch Ausfiihrung der Variationen im Hamiltonschen Prinzip (1) ergeben sich die in der Langs- 
verschiebung u und der Querverschiebung w gekoppelten, nichtlinearen und inhomogenen partiellen 
Differentialgleichungen 


2 
EF(u +9 +1, W.) — ft =0, (5) 


zur Beschreibung der Stabbewegung. Hierzu treten noch die Endbedingungen fir den beidseitig 
quer- und langsunverschieblich und momentenfrei gelagerten Stab, namlich 


uO; 2) =a, 4)-—= 0); (5’) 
w(0, t) = w(l, t) = w,,(0, t) = w,,(1, t) =0. (6’) 


Zur naherungsweisen Lésung der Gleichungen (5) und (6) verwendet man zweckmaBig die 
Kirchhoffsche Annahme’, da diese es gestattet, aus dem gekoppelten System (5) und (6) eine Integro- 
differentialgleichung fiir die wesentlich interessierende Querverschiebung w allein herzuleiten. 
Nimmt man zu dem Zweck an, daB die Dehnung der Stabmittelfaser (2) nicht von der Ortskoordi- 
nate x und nur von der Zeit abhangt, so erhalt man nach Integration iiber die Stablange | bei 
Beachtung der Randbedingungen (5’) fiir diese Dehnung genahert 


a(x, 1) me E(t) => { le 40, w.) de. (7) 


Wie man aus (5) abliest, ist diese Naherung immer dann berechtigt, wenn die von den Querschwin- 
gungen geweckten Langsschwingungen nur Tragheitskrafte vernachlassigbarer GréBe besitzen. 
Dies diirfte schon fiir technische Stabe mittlerer Schlankheit stets zutreffen. 


1 K. Marguerre, Jahrb. 1938 der Deutsch. Luftf.-Forsch. S. I 433. 
* R. Kappus, Z. angew. Math. Mech. 19 (1939), S. 350. 
3 Vgl. G. Kirchhoff, Vorles. iiber math. Physik, 4. Aufl., Bd. 1, S. 444, Leipzig 1897. 
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Mit der durch (7) ersetzten Dehnung (2) geht (6) in die Integrodifferentialgleichung 


l 
EF 2 
EJ W222 —-- | (w. W.. zi 2 dx - (We si We x) sik Bw, + Wp = — p(x, t) (8) 
; 0 ; 


liber, aus welcher sich die Querverschiebung unabhangig von der Langsbewegung naherungsweise 
berechnen la8t. Die anfangliche Vorausbiegung des schwach gekriimmten Stabes wird im weiteren 
stets als sinusf6rmig angenommen, so daB sich die Lage der Stabmittelfaser durch 

TTX 


W(x) = wy sin ae (9) 


beschreiben laBt und w, die gréBte Abweichung von der die Lager verbindenden Geraden bezeichnet. 
Die Querbelastung des Stabes kann sich in beliebiger Weise iiber den Stab verteilen. Doch 


werde die Zeitabhangigkeit der schwingenden Lasten stets als rein harmonisch vorausgesetzt, 
so da8 man fiir p(x, t) mit der Kreisfrequenz w einfacher 


P(x, t) = po(x) cos wt (10) 


schreiben kann. 


3. Rickfiihrung auf ein System gewéhnlicher Differentialgleichungen. a) Reihenansatz fir 
die Lésung. Die nichtlineare Integrodifferentialgleichung (8) 14Bt sich fiir spezielle Schwingungs- 
formen streng auf eine gewdhnliche Schwingungsgleichung zuriickfiihren, wenn keine auBere Be- 
lastung vorhanden ist (p(x, t) = 0). In diesem Falle befriedigen die Produktansitze sin (j/I) - T;(t) 
die vier Randbedingungen (6) und reduzieren auSerdem (8) auf nichtlineare Einzeldifferential- 
gleichungen zur Bestimmung der Zeitfunktionen T;. Es liegt daher nahe, eine allgemeine Lésung 
der vollstandigen Gleichung (8) dadurch zu versuchen, dafi man zunachst die Belastungsverteilung 
P (x) nach den Funktionen sin(j7x/l) entwickelt und dann fiir w(x, t) zu einem Summenansatz 
von der Form 


ee oe UX 
u(x,t) =w,, >) sin] ote T;(t) (11) 
j=l 


iibergeht. In diesem Ansatz bedeutet w,, einen Vorfaktor mit der Dimension der Lange, der bis 

- auf kleine GréBen héherer Ordnung die gréBte vorkommende Querauslenkung der erzwungenen 
Schwingungen bezeichnet. Dieser Faktor ist also geeignet, die GréBenordnungen der vorkommen- 
den Nichtlinearitaten erkennbar zu machen!. Die Ortsfunktionen erfiillen samtlich die Rand- 
bedingungen, so daf die Bestimmung der Zeitfunktionen aus einem gekoppelten System nicht- 
linearer, jedoch gewohnlicher Differentialgleichungen die verbleibende Hauptaufgabe ist. Um diese 
Differentialgleichungen aufzustellen, entwickeln wir auch die Belastung 


pola) = 3 my sin j (12) 


und tragen dann die Entwicklungen (11) und (12) in (8) ein. Da die Vorfaktoren der sin(jx/I) 
fiir jedes j einzeln verschwinden miissen, wird man so auf ein Differentialgleichungssystem fir die 


- Zeitfunktionen T; gefihrt. rps 
Um das System in einer fiir die weitere Rechnung geeigneten Form zu erhalten, fihren wir mit 


der Kreisfrequenz @ der schwingenden Belastung eine dimensionslose Zeitzahlung 
T=ot (13) 


ein. Dabei wird auch die im weiteren sehr wesentliche Erregerfrequenz w zweckmaBig in einem 


dimensionslosen MaBe x gemessen, wobei die Bezugnahme auf den Dimensionsfaktor VE Jiu lt 
der Biegeeigenschwingungen des Stabes naheliegt. Dann erhalt man fiir die Erregerkreisfrequenz 


= 14 
re at (14) 


1 Zur Pinfubrang und Verwendung des Faktors w,, vgl. R. Grammel, Publ. sci. techn. Min. Air. Paris 281 


(1953), S. 49. 
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und fir die Eigenkreisfrequenzen der Biegeschwingungen des Stabes 


= fz 12s); 
i VE SIu E 0 
so daB sich die Frequenzverhialtnisse 
@ zx “ 
= — = — = 1,2, 15 
qj 0; xj 6] 1, ) ( ) 


als Quotienten zweier reiner Zahlen berechnen lassen. Auch die Dampfungskonstante f 1aBt sich 
mit dem logarithmischen Dampfungsdekrement und bei Bezugnahme auf die Eigenfrequenzen 
dimensionslos messen. Berechnet man das logarithmische Dekrement 1}; jeweils mit der Eigen- 


kreisfrequenz @;, so findet man fir die in der Technik vorkommenden kleinen Dampfungen 


reo Wels = 1,2) 55)3 (16) 


wo; Me It 


so daB auch die Dampfungskrafte in dimensionsloser Form in die Gleichungen eingehen. Da die 
Entwicklungsbeiwerte 7; von (12) von der jeweiligen Form der Belastung abhangen, nehmen wir 
diese Koeffizienten allgemein in die Rechnung auf und geben sie erst im Abschnitt 2b) fiir zwei 


interessante Falle explizit an. 
Mit den Abkiirzungen (13) bis (16) erhalt man nun nach formaler Rechnung mit den Ent- 


wicklungen (11) und (12) aus (8) das unendliche Differentialgleichungssystem 
oo le, @) co 
qj T; + d; qj T; =: es 2 Siti T, i aa TD Vik T=—1; COS T (a4; Pt (17) 
I=1 k=1 = 


Hierin sind die Faktoren der nichtlinearen Glieder durch 


Wo Wm 11 &jk 


Wo Wm I €11 1k ay 
Cr = os pt eae (7 “kT 3 ). Stk, = #272 
ie +3 (ze) ) 3 
w?, ke w2, Ke see Ce 
i = Vik G2 ap (j = 2,3,...) 


1 jwo\?\ 4° 
(1+ (= 
i see & ) 
gegeben, wobei mehrfach zur Abkiirzung das Kroneckersymbol 
¢;, =0 fir [7 k-und.ey = Lia fa 6 
verwendet wurde. Die Gré8enordnungen der ¢), ;, yj, und damit der Nichtlinearitaten sind durch 
die Verhaltnisse der Vorauslenkung wy und der zusatzlichen Maximalamplitude w,, zum Flachen- 


tragheitsradius ip bestimmt. Auch die Vorfaktoren der inhomogenen Glieder lassen sich durch 
Bezugnahme auf die Eulersche Knicklast P,; =x? E J/I? in der dimensionslosen Form 


1 Vale Spe Fes 3 
i — : r= = > == Dedareet 19 
: 72 ( aia eS (52) | Pr Wm . me j* Pr Wm (7 ) ( ) 
74 Ua 


schreiben. Die Frequenzverhiltnisse g; in (15) sind mit den dimensionslos gemessenen Eigen- 
kreisfrequenzen 


i 2 4 ; 


zu bilden. Offenbar wird nur x, und damit nur die erste Biegeeigenfrequenz von der Vorkriimmung 
des Stabes beeinfluBt, was in der speziellen Wahl der Vorverformung (9) begriindet liegt. : 
b) Entwicklung der Belastungsfunktion. Von den schwingenden Querbelastungen 
P(x, t) eines Stabes beanspruchen zwei Formen besonderes Interesse. 
Zunachst ist dies der einfache Fall, bei welchem die Entwicklung (12) der Lastverteilung sich 
auf ein einziges Glied reduziert: 


PAX) = 7; sing “- ; 


Diese Lastannahme, die schon am SchluB von Abschnitt 1 erwahnt wurde, hat groBe rechnerische 
Vereinfachungen zur Folge und besitzt vorerst mehr theoretische als praktische Bedeutung. Doch 
wird sich spater zeigen, daB die mit dieser einfachsten Annahme gefundenen Ergebnisse auch in 
allgemeineren Fallen fast ungedndert bestehen bleiben. 


if 


ee ee AS pee 
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Die in der Technik haufigste Belastung diirfte wohl durch pulsierende Einzellasten gegeben 
sein. Die Entwicklung (12) einer solchen Einzellast lat sich in bekannter Weise durch Grenz- 
tbergang aus einer gleichformig iiber ein kleines Streckenintervall verteilten Belastung gewinnen!. 
Fiir eine im Abstand x) vom linken Auflager wirkende Last vom Betrag P ergeben sich die Ent- 
wicklungskoeffizienten 


PANEL canoes 
scan icc ac a 
und daraus die in (17) benétigten Koeffizienten 


1 Pel EU Ne gees reo PUES ere : 
sin : r= sin SSO in oc) ie 20 
3 ( z (=F) Pr wn I eet Pew, Saeed ) 2o 


i 
ip 

Fir den Spezialfall einer in Stabmitte (x) = 1/2) angreifenden Last verschwinden offenbar alle 

Koeffizienten mit geradzahligen Werten von j. Dies hat ersichtlich seine Ursache in der Tatsache, 

da8 die Ausschlagfunktionen sin j z x/I fiir gerades j in Stabmitte einen Knoten haben und daher 

die zugehérigen Teilschwingungen durch diesen speziellen Kraftangriff nicht erregt werden kénnen. 


4, Berechnung der erzwungenen Schwingungen. Durch das System (17) werden erzwungene 
Schwingungen beschrieben, deren naherungsweise Berechnung im weiteren die Hauptaufgabe ist. 
Zu dem Zweck wollen wir von dem unendlichen System die ersten vier Gleichungen verwenden, 
was sicherlich genau genug ist, und auBerdem miissen gewisse Kleinheitsvoraussetzungen ein- 
gefiihrt werden. So soll die Rechnung auf Schwingungsamplituden w,, beschrainkt werden, die 
wesentlich kleiner als der fiir die Biegeschwingungen maSgebende Flachentragheitsradius sind. 
Dann bleiben die in (17) vorkommenden Nichtlinearitaten stets klein, wie man aus den Vorfak- 
toren (18) erkennt. Fiir die Zwecke der formalen Rechnung wird voriibergehend ein kleiner Para- 
meter ¢ mit den folgenden Abkiirzungen eingefiihrt, aus denen man die GréB®enordnung der Koeffi- 
zienten (18) abliest: 


Wo Wm Win 
YoE = ° 2 as 2 (21) 
1 1 2 
Lye gated Le = i poy (22) 
ie 1 (wo\’ te gear AC os 
a Nap 2 \ip 


Die Erregeramplituden r; werden allgemein, aber ebenfalls klein von der Gré®enordnung € an- 


genommen, so daf fiir die Rechnung 


geschrieben werden kann. Da die Dampfungskrafte im wesentlichen von der Werkstoffdampfung 
herriihren, sind die zugehérigen Dekremente klein, was durch die Schreibweise 


| She 4) 
zum Ausdruck gebracht werden soll. Die ersten vier Gleichungen des Systems (17) lauten dann 
¥ 3 9 1 v) 

GT; +4,eq, T+ T+Tye(q ogee cayman Ti) + Bel T? +737, 45, 727,447? T,) 
=—7,ecost,, (25.1) 

7 , 1 1 9 ue 
qa Ts + dz € 42 T2 + eee T, T,+€ (t6 Ti T, Lays T; tig fs Ts Ti Ty) = — Tye cost, 
(25.2) 

Sa eT 1 1 l l 4 s 
q2 T; + d3€ q3 T; +7; + 7g Yoe T, T+ 6 (35 fhe Ty ls let a Loar Tr: T,)=—Tyecost, 
: (25.3) 

7 7 1 1 1 9 i} i) 3S 
‘ee , (25.4) 
1 Vgl. z.B. J. Széb6, Héhere Technische Mechanik, S. 118, 2. Aufl. Berlin/Géttingen/Heidelberg 1958. 
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Diese Differentialgleichungen sollen zur Gewinnung einer Naherungslésung mittels der Stérungs- 
rechnung ausgewertet werden. Zu dem Zweck entwickeln wir die gesuchten Zeitfunktionen T;(t) 
nach dem kleinen Parameter ¢ und setzen 


® @) (2) 
T= 1; Tye ee j= 1, 2,3; 4) (26) 


Um die fiir nichtlineare Schwingungen charakteristische Frequenz-Amplituden-Abhangigkeit 
aufzudecken, miissen wir auch die Erregerfrequenz x in der Nahe einer vorgegebenen, festen Fre- 
quenz entwickeln. Dabei ist es zweckmaBig, die Quadrate der Frequenzverhaltnisse q; unmittelbar 
zu entwickeln. Soll z. B. eine Erregerfrequenz x in der Nahe einer festen Eigenfrequenz x; be- 
trachtet werden, so schreiben wir 


\2 (1) (2) : 
gta Sat ae ta (27.1) 
wahrend fiir die tibrigen Quadrate der Frequenzverhaltnisse dann 
x \? ny 3 xy 2 (1) (2) . . . 
(H=g =(Z) ato =(Z) +getget  G+9 (27.j) 


gesetzt werden muB. In der weiteren Rechnung kann nur iiber die a geeignet verfiigt werden, 
wahrend die qi (j + i) dann durch (27.j) davon abhangig zwanglaufig gegeben sind. 

Bei der Berechnung der erzwungenen Schwingungen bieten zwei Erregerfrequenzbereiche beson- 
deres Interesse. Dies ist zunachst natiirlich der Frequenzbereich in der Nahe der tiefsten Eigen- 
frequenz. Die zugehérige Schwingungsgleichung (25.1) mit r, + 0 nimmt auBerdem deswegen eine 
bevorzugte Stellung ein, weil die Vorverformung des Stabes hierin von besonders grofem Einflu8 
ist. Der nachsthéhere Eigenfrequenzbereich, der nicht durch die Art der Vorverformung bevor- 
zugt ist, ist der Bereich in der Nahe der zweiten Eigenfrequenz, und diese soll daher als reprasen- 
tativ fiir eine nicht ausgezeichnete Erregerfrequenz ebenfalls betrachtet werden. Dabei suchen 
wir solche erzwungenen Schwingungen auf, welche die gleiche Periode t = 27 wie die erregende 
Kraft haben. Andere Werte der Erregerfrequenz sowie Schwingungen mit anderen Perioden 
(subharmonische Schwingungen, mit deren Vorkommen man grundsatzlich rechnen muf,) sollen 
in dieser Arbeit ausdriicklich beiseite gelassen werden. Da bei den folgenden Rechnungen wesent- 
lich die Frequenz-Amplituden-Beziehungen interessieren, soll zur Vereinfachung die Dampfung 
vernachlassigt werden. Bei einem spdter notwendig werdenden Stabilitatsnachweis muf die 
Dampfung allerdings in die Rechnung aufgenommen werden. 


a) Erregerfrequenz in der Nahe der ersten Eigenfrequenz. Tragt man die Ent- 
wicklungen (26), (27.i) und (27.j) in die nichtlinearen Differentialgleichungen (25, 1—4) ein und 
ordnet nach Potenzen des Parameters ¢, so erhalt man fiir jede Potenz von ¢ jeweils vier inhomogene 
lineare Differentialgleichungen, aus denen die mit den Stérgliedern gleichperiodischen Lésungen 
rekursibel zu berechnen sind. 

Fiir die von ¢ freien Glieder erhalt man in dem jetzt betrachteten Fall, daB die Erregerfrequenz 
in der Nahe der ersten Eigenfrequenz liegen soll (i = 1), die Differentialgleichungen 


(0), (0) x¢,\2 (9) (0) ¢ 
f+ T, =0., =) 2a + 7T;=0 Gf = 2535029) (28) 
J 


deren periodische Lésungen der Periode 2 7 aufzusuchen sind. Als periodische Lésung der ersten 
Gleichung (28) hat nur 

(0) 

T= cos t (29.1) 


Interesse, da ja insgesamt erzwungene Schwingungen berechnet werden sollen und die Erreger- 
glieder in (26, 1—4) das gleiche Zeitverhalten wie (29.1) aufweisen. Man iiberzeugt sich leicht, 

. . . (°) 
daB die formal gleichfalls mégliche Ausgangsfunktion T, = sint im folgenden Rechnungsgang 


zu keiner periodischen Lésung fiihren wiirde. Die verbleibenden Gleichungen (28) lassen als perio- 
dische Lésungen der verlangten Art nur 


(0) 


Zu. 


Tj=0 (jf =2,3,4) (29.j) 
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Weiter findet man nun fiir die Glieder erster Ordnung in ¢ die inhomogenen Differentialglei- 
chungen 


re Q) 3 (0) (1) (0) 
ee — he Fi cost — 9, T,, (30.1) 


we, : 
a Tj + T; =— 7; cost fea 2/354.) , (30.j) 
wenn man (29.j) beachtet. Durch Einsetzen von (29.1) in (30.1) kommt 


na (1) 
T,+ T, ee +(—7 +4) cos7—3 7, cos2T, 


wovon ey eine periodische Lésung aufzusuchen ist. Es mu also tiber den Entwicklungs- 

koeffizienten i, so verfiigt werden, daB das Resonanzglied mit cos t verschwindet, da dieses in die 
(1) 

Lésung T, nichtperiodische Bestandteile hineintragen wiirde. Man erhalt daher 


(h) 7, Pt 
qi =) (31) 


und als Lésung 


(1) 3 1 
1, = — =z [y+ 7 Loc0s 27. (32.1) 


(1) 
Aus (30.j) berechnen sich die T; als erzwungene Schwingungen zu 


(cae ean 


1) 
[; = Rae, cos tT = a4 COS T (Gj = 253.4); (32.}) 
() 


so daf} also wegen ,/x; == 1 in den Schwingungen T; Resonanz nicht eintreten kann. 

Zu den Lésungen (32.1) und (32.j) kénnte man noch jeweils die Lésungen der homogenen 
Gleichungen mit beliebigen Integrationskonstanten hinzufiigen. Uber diese Konstantéh muB 
daher durch eine zusdtzliche Festsetzung verfiigt werden, damit die Rechnung eindeutig wird. 
Es ist zweckmaBig, auf die Mitnahme der freien Schwingungen in jedem Rechenschritt iitberhaupt 
zu verzichten und also alle Integrationskonstanten gleich Null zu setzen. Dadurch wird die Rech- 
nung formal am kiirzesten und der Nachteil, da sich spater der Anfangsausschlag des Stabes zur 
Zeit t = 0 als Reihenentwicklung ergibt, kann in Kauf genommen werden. 


Fur die Glieder zweiter Ordnung in « folgen die Gleichungen 


(2) (2) 3 (0) @) l (0) qd) @ @. @: 


Pe eo ey, Le, T,—q T;., (33) 
et, hy, THT, 9, Ty, (34) 
ea ty Ph TL (35) 
Bee a a | (36) 


die sich mit den Lésungen (29) und (32) leicht auswerten lassen. Aus (33) ergibt sich zundchst 
mit nach (31) bekanntem q, 


(2) (2) 2 3 1 
fa Tole ra r+ cos +1 Tyoos2t +(—4 P35 Z| cos 3T , 


woraus zuerst 


(2) 15 3 
Ge aay Fa a, F: (37) 
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(2) 
folgt und sodann T, als streng periodisch berechnet werden kann: 


(2) Fal 3 1 
mala fe bores) | 
UT g 008 2T + laee 128 


r+ r| cos 3T. (38) 


Da verabredungsgem4B freie Schwingungen niemals zu den erzwungenen Schwingungen zugefigt 


(n) 
werden sollen, ist also auch (38), wie alle T,(n = 2,3,...) frei von der ersten Harmonischen 
(0) 
cos T, die also allein in T, vorkommt. 
. Bei der Auswertung der Gleichungen (34), (35) und (36) ist zu beachten, da auch die rechts- 


seitig vorkommenden Frequenzparameter q, (j = 2, 3,4) als bekannt anzusehen sind. Mit den 
Lésungen (29) und (32) sind also die rechten Seiten vollstandig gegeben und man erhalt 


(2) 1) A qa) qd) @ 
C= ee Yo ee Ces feet, cos 2 2 ves , COST , 
16(1—4(2)) ve) 
Ho, a9 
(2) (1) (1) qd) dd) 
= 35 Yo ga — og cos 2 t + 2 AE, cost , 
36(1—4/ ) Fee 
3 Hy 
(2) (1) a) qd) qd) 
T,= Se Vy 44. — eel cos2T + A, COST . 
o4(1—4(2)) i 
4 %4 


Aus diesen Lésungsteilen erkennt man, dai iiber die bereits ausgesprochene Bedingung 
x,/%; (J = 2, 3, 4) hinaus auch noch %/x; + 1/n (n ganz, j = 2, 3, 4) gelten muB, so daB sich darin 
die bekannten mathematischen Schwierigkeiten der Stérungsrechnung mit den kleinen Divisoren 
anzeigen. Fir die hier nur interessierenden Effekte iiberwiegender GréBenordnung ist diese mathe- 
matische Problematik jedoch belanglos und die Berechnung des Quotienten q?(w,,) soll mit der 
nachfolgenden Bestimmung von q{3) abgeschlossen werden. Zu dem Zweck wird von den Gliedern 
mit dem Faktor e° lediglich 


sid (3) (3) 3 (72 (0) T,) (1) 9 (1) (1) 
een a 2. ba ay meat? et fer es ra Pr 
3 al, OO. 0:0 er Saag! 
ae oy: dal 1s, T-—¢, T; (39) 
. . . (9) 
gebraucht. Auf der rechten Seite ist allein der bei T; = — cost als Vorfaktor stehende Frequenz- 


(3) : 
parameter q, unbekannt und nur dieser soll berechnet werden. Ausfiihrlicher geschrieben lautet 


Gleichung (39) 


ieee 3 aiken ¥ 
T+ 7, =—27(—27, +L ryeos 22] 
\ 


3 wire! Be 3 1 
— 4 Toco 7(— 5° cos 21 + (35613 + a5g 1) c08 3 7] 


eb) 9 q@) @ 
ee 2 2 2 2 2 
158 a2 , cos* tT — dy a, cost 4 I az , cos? t 


3 1 1 \ 
i | 5 cos2 1)(— 5 Ty + Lycos 27} 


(3) QO gg @ 
+q cost—q, T;—4q, T,, 


ey: EN (3) od 
ea jedoch fiir die Berechnung von q, nur die cos t-Terme von Bedeutung sind, so daB T, und 


Tj nicht ausfiihrlich niedergeschrieben zu werden brauchen. Durch Ordnung nach den Harmoni- 
schen erkennt man, daB rechtsseitig als erste Harmonische allein 


| pee 3) 
Gi ae +4) cos T 
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. . . . 3 
erscheint, und diese mu8 auch hier durch Wahl von 7. zum Verschwinden gebracht werden, damit 
die Lésung rein periodisch bleibt. Mit 


(3) Le 
Gee eae (40) 


ist also die Frequenzgleichung (27.i) fiir i = 1 
(1) (2) (3) 
GM@=l+qe+que+ya é° + O(e4) 
bis auf Glieder der GréBenordnung e* bekannt. Mit den Werten (31), (37) und (40) fiir die Koeffi- 


zienten wird 


; ip tae = 185 es ; 
@=l+ne4 (ig? 3g F8) 8 — Th + O64). (41) 


In den urspriinglichen Bezeichnungen und in der obigen Naherung kann (41) auch in der Form 


S oar 1 me oe 1 15 
aw)! qos 1 a) 7? Pp Wy, | Cs 1 yea 32 
2; ip 


1 "| we, 
1 /w)\? ip \i 
Sore eee, 


ZY UF 
Ire 1 n, 2 iL Wo |2 Wm 
= = 5 : = 0 42 
8 peat w Prip ere | ip (42) 
2 Up 2 \ ip 


geschrieben werden, so daB man diese Gleichung fiir w,, als eine Frequenzgleichung vom Duffing- 
schen Typ erkennt, die wir sogleich diskutieren. 

Zu je einem der drei aus (42) zu berechnenden Wurzelwerte von w,, gehért eine spezielle er- 
zwungene Schwingung des gekriimmten Stabes, gekennzeichnet durch einen bestimmten Schwin- 
gungsausschlag an einer beliebigen Stabstelle x. Dieser Ausschlag zur Zeit t = 0 werde im weiteren 
als Amplitude bezeichnet, da die zugehérige Auslenkungsgeschwindigkeit w,(x, t = 0) verschwindet, 
denn alle Zeitfunktionen T;(t) sind durch reine Sinusreihen gegeben. Die Schwingungsamplitude 

_an der Stelle x ist also in unserer Naherung durch 


Ry ee G (a) 2 
want = 0) = Wm 2 sin j T;(0) + e T;(0) + é T;(0) (43) 


a 


bestimmt, wenn man sich auf den durchgefiihrten viergliedrigen Produktansatz beschrankt und 
die Zeitfunktionen unter Vernachlassigung kubischer und héherer Potenzen von ¢ einsetzt. Z. B. 
erhalt man fir eine in Stabmitte angreifende schwingende Einzellast zuerst aus (20) mit x) = I/2 
die benétigten Koeffizienten r; und dann mit (32.1) und (32.)) fiir die Amplitude in erster Naherung 


w(x,t = 0) = w,, oe i ee 
aaa . 


TOK 2 eat th 
+ sin 3 i [re ee pap] ) 


Die Amplitude (43) 14Bt sich also, wie aus den gegebenen Formeln und dem vorstehenden Beispiel 
im einzelnen hervorgeht, eindeutig aus w,, berechnen. 

Charakteristisch fiir die vorliegende Schwingung ist somit die in (42) liegende Mehrdeutigkeit 
beziiglich w,,. Die algebraische Gleichung (42) hat die Form 


EAN r 


1—(S) tog te Buh te C tty = 9, (44) 
Wy Wm 

wenn A, B, C abkiirzende Konstanten bezeichnen. Man erkennt aus (42) unmittelbar, daB stets 
eA>0 und 2 C>0 ist. Der bei w®, stehende Vorfaktor ¢? B ist positiv fiir wy< ip/V2 und 


negativ fiir wy) > bel 2: Nehmen wir beispielsweise den letzteren Fall an, so erkennt man aus einer 
einfachen Kurvendiskussion leicht, daB (44) graphisch durch eine Kurve dargestellt wird, die in 
Abb. 2 qualitativ skizziert ist. Dabei ist, was man zunachst fiir selbstverstandlich halten wird, 


’ 


A RENE NWP 
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n, = const vorausgesetzt. Abb. 2 entspricht qualitativ véllig der Amplituden-F requenz-Kurve 
des einfachen (Duffingschen) Schwingers mit unterlinearer Federcharakteristik. Hatten 
Wy < ip/V2 angenommen, so hatte sich e? B > 0 ergeben, und statt der Abb. 2 hatte man die 
Amplitudenkurve der iiberlinearen Schwingung erhalten. Entsprechend der Verwandschaft mit 
dem Duffing-Schwinger entspricht der punktierte Ast der Kurve von Abb. 2 einer instabilen 
Schwingung und nur die ausgezogenen Aste gehéren einer stabilen Schwingung zu, was allerdings 
noch nachzuweisen ist (vergl. Abschnitt 5). Wie man aus den Kurven ferner abliest, stoBen der 
, stabile“ und der ,,instabile“’ Ast in Punkten mit lotrechter Tangente zusammen. Die Schwingung 
muB also das vom einfachen Schwinger her bekannte Amplituden-Kipp-Phanomen zeigen. 


AVm, 


stabil 
——-— instabil 


fs) 


a 
£ 


0 
te 
/ 
/ 
: / 
Ba 
/ ——— instabil 
t 
/ 
74 
4 
Abb. 2. Graphisches Bild der Frequenz-Amplituden-Beziehung (42) Abb, 3. Graphisches Bild der Beziehung (42) 
im unterlinearen Fall w.) > ip/V2 (Lastamplitude konstant). fir w/w, = konst. < 1 und w,> ip/)2. 


Die enge Verwandschaft des vorliegenden Schwingungsproblems mit dem statischen Durch- 
schlagproblem wird jedoch erst aus einer anderen Darstellung deutlich. Setzen wir namlich in (42) 
baw. (44) (w/w,)? = const < 1 und lassen dafir 7, und damit nach (12) die Amplitude der auBeren 
Belastung variieren, so wird (42) durch eine Kurve nach Art der Abb. 3 wiedergegeben. Wieder 


ist dabei wy > ip/V2 gewahlt. Und auch in Abb. 3 gehéren die ausgezogenen Teile der Kurve zu 
stabilen und die punktierten zu instabilen Schwingungen. Abb. 3 ist so zu deuten: Steigert man 
die Lastamplitude von Null an, so wachst der Schwingungsausschlag, fiir den wir das Ma’ w,, 
benutzen, zundchst stetig mit; man erreicht dann aber einen Indifferenzpunkt mit horizontaler 
Tangente, von dem aus sich die Last-Ausschlag-Kurve ,,instabil“ nach abwarts senkt, wahrend 
die Schwingungsamplitude nicht dieser Kurve folgt, sondern auf einen stabilen Ast der Kurve 
iiberspringt. Natiirlich ist das nichts anderes als das oben schon erwahnte Kipp-Phanomen in 
einer anderen Koordinatenebene gesehen. Die Ahnlichkeit von Abb. 3 mit der Last-Amplituden- 
Kurve des statischen Durchschlagens! springt in die Augen, woraus wir die Berechtigung ableiten, 
auch im kinetischen Fall von Durchschlagen zu sprechen. 


Allerdings hat die Analogie zwischen dem-statischen und dem kinetischen Problem auch ihre 
Grenzen. Diese bestehen hauptsachlich darin, daf8 das Phanomen des kinetischen Durchschlagens 
oder Kippens auch dann noch, und zwar iiberlinear bestehen bleibt, wenn w, = 0, der Stab also 


gerade ist®, wahrend das Phanomen andererseits verschwindet, wenn w) = ip/V2 wird. Dieser 
spezielle Wert der Vorauslenkung bestimmt die Grenze zwischen dem unterlinearen und dem 


tiberlinearen Fall der Zwangsschwingungen und daher erfolgt dann iiberhaupt kein Kippen der 
Amplituden. 


_b) Erregerfrequenz in der Nahe der zweiten Eigenfrequenz. Da die erste Kigen- 
frequenz und die zugehdrige Schwingungsgleichung (25.1) durch die gewahlte spezielle sinusférmige 
Vorverformung ausgezeichnet ist, soll die Frequenz-Amplituden-Gleichung noch fiir eine Erreger- 
frequenz in der Nahe einer anderen, nicht ausgezeichneten Eigenfrequenz berechnet werden. Fir 


* Vgl. Abb. 3 der in FuBnote 1 auf S. 301 genannten Darstellung von C. B. Biezeno und R. Grammel. 


* Dieser Sonderfall wurde ohne Hinweis auf die Verwandtschaft mit dem Durchschlagen getrennt be- 
handelt: E. Mettler, Z. angew. Math. Mech. 31 (1951), S. 263. 
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diese Rechnung wird die Erregerfrequenz als in der Nahe der zweiten Higenfrequenz gelegen voraus- 
gesetzt und im iibrigen die formale Stérungsrechnung ganz entsprechend der bereits im Abschnitt a) 


-beschriebenen durchgefiihrt. 


In (27.i) ist jetzt i = 2 zu setzen, und wenn man dementsprechend in (28) die Indizes 1 und 2 


vertauscht, findet man die von ¢ freien Glieder einer mit 27 periodischen Lésung zur Beschreibung 


der erzwungenen Schwingungen fiir 7, + 0 
(0) (0) 


P5608, % 5 Pa (j= 1, 3, 4).. 


Die Glieder erster Ordnung in ¢ verlangen entsprechend (30.1) und (30.j) hier etwas anders 


42 ®, Or : 
a f+ f=, Ti — 1, cost, (45.1) 
M 7 1) 
Pie Tees recoan~ Ge Te (45.2) 
i \e Gee ol) is 
(2) T; + T; =—1T;, cost (G34). (45.j) 
oy 
Aus (45.2) kann mit T) = — cos t zunachst 
Nha 
q2 =T (46) 
und sodann 
(1) 
T,= 0 


(1) 
gefunden werden. Fir T, gewinnt man aus (45.1) 


2 (1) (1) rT iB 
iy ean 8 Ley — cos 2T 
(2) ob emt ke 5 Ty COST 3 C08 ; 
~ so daB 
(1) T 
I COS T — I cos 2 4h 


wird. Die Lésungen der Gleichungen (45.j) lauten 


(1) ere (1) 3 
T.= J cost = bj,, cost (j= '354) 


| ( ) 
J 


Die Glieder zweiter Ordnung in ¢ fiihren wegen T, = 0, T; = 0 (j = 1, 3, 4) auf die Differential- 


: gleichungen 


(2) (2) (1) (0) 1 (2) (0) 
T,+T, =—*2 1, T,—> Tl—-G Te (47.2) 


2 (2) (2) qa) @ ; : 


ed) 


deren rechte Seiten wiederum durch die vorhergehende Rechnung bekannt sind. Fiir j = 1 findet 


man 
2 (2) (2) 1) Tr ele 
ry zs a T Seay Cos T aay Be 
Bpe Ne 


\ 44, 


mit der Lésung 
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Aus (47.2) laBt sich der fiir die Frequenzgleichung besonders interessierende Entwicklungskoeffi- 
zient i bestimmen. Man erhalt zunachst 


2) ef) ; r Tr Ih 
” Yo ov eee ees pe OEE CORT. 
T, + WES Se lee a2 COST 3 


Lys 1 (2) 
=e cos T + | cos 3 2) + Jz COST , 


und da auf der rechten Seite das Glied mit cos t zu verschwinden hat, muB also 


% = le es (48) 


, = ee a I \2 
16 16 2fr—a(%) 
1 


(3) Y (2) (0) ; 
Der Wert von q, l48t sich schon mit Kenntnis von T, und T, allein berechnen und soll daher 


noch angegeben werden. Aus 


sein. 


3) @) , 2 (3 «ay @ 
T, +7, =—* T, T,— q@ T,— % T; 
folgt namlich 
(3) (3) (1) — 9 (1) (3) (1) (2) 
T,+7T,=— Ye [a ceed +9 COS 7 ee le cos2T\ cost + g, cost — q, Ty. 
pee 
1—/— | 1—4 
ass) ta 
(3) 
Also hat q, den Wert 
(3) ber iis % Ty Vo (49) 


Die Frequenzgleichung (27.i) fiir.i = 2 lautet demnach mit (46), (48) und (49) 


is i%s 3 E. 1 (1) r. 
g@=lt+rmRe4 (e~e"( + 2-4) ae sees is aaa & + Ole") (50) 


ape “| 
Ae: en 
und entspricht der Gleichung (42). 

Die fiir Erregerfrequenzen w in der Nahe der zweiten Eigenfrequenz jw, geltende Gleichung (50) 
ist wiederum vom Duffingschen Typus und kann zur Bestimmung von w,, verwendet werden. 
Zu jedem Wurzelwert gehért dann wieder eine spezielle erzwungene Schwingung, die durch den 
Stabausschlag an einer beliebigen Stelle x zur Zeit t = 0 gekennzeichnet werden kann. Fiir eine 
Berechnung dieser Auslenkung in der Naherung 0(c*) benétigt man aufer dem bereits angegebenen 


(2) (2) 
T, noch die Lésungen T; fiir j = 2, 3, 4: 


(2) 
T, ort 25 “ont = cos 27 + oe vol al cos 3T, 
16 ( ‘| 48 1— (2) 256 (3) 
4 oy 1 
(2) GF 
T; Z cos T G = 1,3, 4) 


Damit kann dann, wie fiir w © «, geschehen, der Stabausschlag an der Stelle x zur Zeit t = 0 
mit der Reihenentwicklung (43) angegeben werden. Da auch hier w(x, t = 0) = 0 ist, kann man 
durch diese Angaben die zugehérigen erzwungenen Stabschwingungen auch fiir m © wy geniigend 
genau kennzeichnen. Die Formeln sollen jedoch nicht im einzelnen niedergeschrieben werden, 


da Besonderheiten gegeniiber der im Abschnitt a) ausfiihrlich wiedergegebenen Rechnung nicht 
auftreten. 
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5. Stabilititsuntersuchung. Aus den Frequenzgleichungen (42) und (50) kénnen fiir eine vor- 
gegebene Erregerfrequenz jeweils entweder ein Wert oder drei Werte des Amplitudenfaktors w,, 
berechnet werden, so da man mit diesen Werten auf einen mehrdeutigen Schwingungszustand 
gefiihrt wird. Es ware nun noch fiir jeden solchermaBen berechneten Wert von Wp» nachzuprifen, 
ob der zugehérige Schwingungszustand stabil ist und somit physikalische Realitat besitzt. Eine 
derartige Stabilitatsuntersuchung ist mathematisch auSerordentlich weitlaufig, so daB im folgenden 
nur eine erste Naherung durchgefiihrt werden soll. Diese Naherungsrechnung verwendet die von 
J. Haag! auf Schwingungssysteme mit endlich vielen Freiheitsgraden erweiterte Methode der 
langsam veranderlichen Phase und Amplitude und gestattet in der damit nur erreichbaren nieder- 
sten Naherung die Riickfiihrung der Stabilitatsfrage auf die Stabilitat der erzwungenen Schwin- 
gungen des Duffingschen Schwingers mit nur einem Freiheitsgrad. Da dieses einfachere Schwin- 
gungs- und Stabilitatsproblem eingehend untersucht ist, wird mit dieser Riickfiihrung die Aufgabe 
naherungsweise gelist. 

Fiir die Stabilitatsuntersuchung ist das Abklingen einer eingeleiteten Stérung durch kleine 
Dampfungskrafte so wesentlich, da die bisherige Vernachlassigung dieser tatsachlich stets vor- 
handenen Krafte aufgegeben werden mu8. Die Frequenz-Amplituden-Gleichungen werden sich 
allerdings durch die Mitberiicksichtigung einer sehr kleinen Dampfung nur unwesentlich andern, 
so daf die bisherige Vernachlassigung der Dampfungen eine gewisse Berechtigung hatte. 

Zur Durchfiihrung der formalen Rechnung wird das allgemeine Gleichungssystem (17) mit 
einem kleinen Parameter A in der Form 


u" 1 A = , wv N= Nee 2 — A F al 
T; Te itg = 54; eat 2 Sibi T, Tie Uh Ti, — 7; cos t — @ j (51) 


qj J 
Peo ee, N) 
geschrieben, worin die rechten Seiten durch die Bezeichnungen 
Gj= dA, buy =Sngd>  Vye =k, r= TA 


als klein von mindestens der GréBenordnung ( gekennzeichnet sind. Mit den teilerfremden ganzen 


Zahlen m; und M und den Konstanten ¢; werden die Frequenzverhaltnisse durch 1/q; = ;/x 


= m; yl +A e;|M, (j =1,2,...,.N) angenahert und dementsprechend die Lésungen Tj in der 
Form 


Ty =yjcos(Mztaj) = 1,2,3,-2-50) 


mit wenig von der Zeit abhangenden Amplituden y,; und Phasen x; angesetzt. Diese Langsam- 
_ veranderlichkeit soll durch die zusatzliche Forderung 


< m: m: . 
Ty = — "yy sin(T 2 + 2) =— Ey; sin @, (52) 


mit @; = a tT + %; zum Ausdruck kommen, so daB also die Nebenbedingungen 
— y; x; sin ®; + y; cos O; = 0 Gj eel} ds Bey) (53) 
Pestehon miissen. Durch weitere Differentiation findet man aus (52) den Wert T; und kann damit 
an Stelle von (51) 
—y; x} cos ®; —y; sin Bj =A J (—e; y; cos ®; + F;) jee IN) (54) 


schreiben. Die Gleichungen (53) und (54) lassen sich sodann nach y; xj und y; auflésen, was auf 


m: 
Hi Pigg N08 Pia F208 P; » 65) 
yj =A 52 (ey; 208 ®; — F}) sin ®, (Gj =1,2,..-,N) 
fihrt. 


i i i i i Arbeit 
1 M t bratoires, Bd. 2, S. 210, Paris 1955. Die Methode wurde in der 
F =A eset ids ne ate 26 (1956), S. 53, ausfiihrlicher dargestellt, weshalb wir uns hier kurz fassen. 
kénnen. Vergl. auch F. Weidenhammer, Z. angew. Math. Mech. 36 (1956), S. 235. 
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Aus diesem Differentialgleichungssystem erster Ordnung kann man durch Mittelung der rechten 
Seiten iiber eine Periode 2 M ein System ableiten, in dem die Zeit nicht mehr explizite auftritt 
und aus dem sich in definierter Naherung die Frequenz-Amplituden-Gleichungen herleiten lassen 
und auch die Stabilitatsfragen entschieden werden kénnen. Fiir diese nun vorzunehmende Mitte- 
lung ist zu beachten, daB die Erregerfrequenz jeweils in der Nahe einer der beiden Eigenfrequenzen 
liegen soll. Es ist also fiir i = 1 oder t = 2 

SS m.- —— ==, 4 
<=yi+in. aml tite (f = Ay 2p ndeyeen) 
zu setzen, so daB die Mittelwerte der rechten Seiten von (55) leicht zu bilden sind, da man hierbei 
x;,y; wegen ihrer Langsamveranderlichkeit als Konstante ansehen darf. Man erhalt an Stelle 
von (55) die einfacheren Gleichungen 


a =i(—4, é + F sin x), (56.1) 
/ AIS cd —_— 3 — T; BES 
Ni Xi =a(% Di Vin Vi +e Vii Hi + O08 MT E ah (56.2) 
k=1 
+i 
¥ m. do y; 
7 ie fa) (57.1) 


9 iS = oes See fats 2 
95% =A zy B , Vink +e VII + iy (a4, 25 eee (57.2) 


aa) 


worin die letzten Gleichungen erkennen lassen, daB die Amplitudenwerte y; = 0 nicht ausgeschlossen 
zu werden brauchen. Aus den Gleichungen (56) und (57) lassen sich die erforderlichen Schliisse 
auf das Stabilitatsverhalten in dem hier benétigten Umfang ziehen. 

Zunachst berechnet man fiir stationare Schwingungen nochmals, jedoch nun nur noch in erster 
Naherung, aus den Gleichungen (56) und (57) die Beziehungen zwischen den Amplituden y; und den 
Phasen. Denn fiir den eingeschwungenen Zustand miissen die Amplituden y; und die Phasen x; 
simtlich feste Werte annehmen. Man liest zunachst aus (57.1) ab, daB in dieser Naherung alle 
Amplituden y;(j = 1, 2,..., + i) verschwinden, so daB die nur noch verbleibenden Gleichungen 
(56.1) und (56.2) eine Aussage iiber die Amplitude y; und die zugehérige Phase x; allein liefern: 


—d, +4 sin x, = 0, (58.1) 


ta 


Syixt +4 cos x, +e, =0. ", (88:2) 


Diese Gleichungen erhalt man auch, wenn man den einfachen Duffingschen Schwinger allein be- 
trachtet. Tatsachlich ist ja auch in dem System (51) fiir j =i allein und), ; = 0 (I, k = 1, 2,..., N) 
die Duffingsche Schwingungsgleichung enthalten und diese fihrt also auf den gleichen, durch 
(58.1) und (58.2) beschriebenen Zusammenhang zwischen Amplitude und Phase. Aus diesen 
Gleichungen gewinnt man mit 1 — q? =/.e; + 0(/?) nach Elimination der Phase x, die Frequenz- 
Amplituden-Gleichung in der gewohnten Duffingschen Form 
2 
ga 14/o— at td yust. (659) 
- Fir verschwindende Dampfung (d; = 0) miiBte diese Gleichung identisch sein mit der aus- 
fibrlich berechneten Frequenz-Amplituden-Gleichung (41) baw. (42), wenn die Methode der lang- 
sam veranderlichen Phase und Amplitude eine gleichermafen vollstandige Naherung zu berechnen 
gestattete, wie die Stérungsrechnung fiir ungedampfte Schwingungen. Tatsichlich gibt (59) fir 
d;~ 0 nur die ersten drei rechtsseitigen Terme von (41) richtig wieder. Wir miissen uns jedoch 
mit dieser eingeschrankten niedersten Naherung fiir die Stabilitatsbetrachtung begniigen. In 
dieser Naherung wird man also auf ein einfaches Duffingsches Schwingungsproblem gefiihrt und 
man weif daher, daB von den maximal méglichen drei Amplitudenwerten die Schwingung mit der 
Amplitude mittleren Betrages einer instabilen Bewegung zugehért und folglich keine physikalische 
Existenz besitzt. Diese Stabilitatseigenschaften wurden in die Betrachtungen des Abschnittes 4a) 
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bereits aufgenommen, da sie in der hier durchgefiihrten Naherung auch fiir das vollstandige System 
gelten. Dies kann man aus der Tatsache entnehmen, daB die Schwingungen mit den Amplituden 
¥i + 0,¥; =0G = 1,2, ..., N, + i) stabil sind gegen kleine Stérungen. 

Betrachtet man kleine Abweichungen y; und €; der Amplituden und der Phasen von ihren 
stationdren Werten in (56) und (57), so erhalt man fiir diese Stérungen nach Linearisierung in Bezug 
auf die StérungsgréBen die Differentialgleichungen 


(oe 0 
, Ye 6 eve Uh: 

wei =z Viyi ni — = sing, & + ¢; #), (60.2) 
, ey. : : 
Dj aa Ae Ge 1) 2,2 IN (61) 


wenn man y; = 0(j = 1,2,...,N, + i) beachtet und eine leere Aussage iiber die Phasen é, 
sogleich fortlaBt. Aus (61) folgt namlich, daB die Amplitudenstérungen 7;(j =: i) samtlich expo- 
nentiell auf Null hin abklingen, so daB etwaige zugehérige Phasenstérungen belanglos sind. Die 
verbleibenden Stérungsgleichungen (60.1) und (60.2) sind unabhangig von den Stérungen Nj §; 
(j + 7) und enthalten nur die stationare Amplitude y; und die zugehérige Phase x;. Aus diesen 
Stérungsgleichungen kénnte man also das Stabilitatsverhalten des einfachen Duffingschen Schwin- 
gers herauslesen. Diese Kigenschaften sind jedoch bekannt!, und man findet hier in den Grenzen 
der durchgefiihrten Naherung, da diese Stabilitatseigenschaften auch fiir das vollstandige System 
giltig bleiben. 


6. Zusammenfassung. In der Statik ist das Durchschlagen des gekriimmten Stabes ein wohl- 
bekanntes Stabilitatsproblem. Ein ganz entsprechendes kinetisches Stabilitatsproblem entsteht, 
wenn der gekriimmte Stab unter der Einwirkung einer schwingenden Querbelastung steht. Mit 
den Annahmen der technischen Biegelehre und der Kirchhoffschen Naherung fiir die Stablangskraft 
kann man zur Beschreibung der Stabschwingungen eine inhomogene nichtlineare Integrodifferential- 


_gleichung fiir die Biegeschwingung allein aufstellen. Ein Lésungsansatz mit zeitabhangigen Funk- 


-tionen fiihrt diese partielle Integrodifferentialgleichung fiir eine beliebige Lastverteilung auf ein 


_ System nichtlinearer gewéhnlicher Differentialgleichungen zuriick. Aus diesem System lassen sich 


die erzwungenen Schwingungen mit einer Stérungsrechnung berechnen; sie haben im wesent- 


lichen den gleichen Frequenz-Amplituden-Zusammenhang, den man schon vom Duffingschen 


Schwingungsproblem her kennt. Das Stabilitatsverhalten dieser Schwingungen lie sich mit Hilfe 
der Methode der langsam veranderlichen Phase und Amplitude naherungsweise ebenfalls aus den 
entsprechenden Eigenschaften des Duffingschen Schwingers verstehen. Es sind daher Amplituden- 
spriinge zu erwarten, die sich in einer Weise deuten lassen, die erkennen laBt, daB sich das statische 
und das kinetische Problem weitgehend entsprechen. 


Aus dem Institut fiir Mechanische Schwingungstechnik der Technischen Hochschule Karlsruhe. 


(Eingegangen am 21. November 1959.) 


Anschriften der Verfasser: Prof. Dr. E. Mettler, Karlsruhe-Durlach, Geigersbergstr. 12, 
Prof. Dr. F. Weidenhammer, Karlsruhe-West, Hertzstr. 16, Bau 40. 


1 Vgl. z. B. J. Haag, Les Mouvements vibratoires, Bd. 1 S. 170, Paris 1952. 
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Einige Betrachtungen zur Beschreibung von Vorgangen 
in der klassischen Kontinuumsmechanik 


Von Th. Lehmann 


1. Einleitung. Bei der Beschreibung von Vorgingen in einem Kontinuum unterscheiden wir 
die Eulersche und die Lagrangesche Betrachtungsweise. Diesen beiden Betrachtungsweisen liegen 
zwei verschiedene physikalische Standpunkte zugrunde, die sich etwa folgendermaBen skizzieren 
lassen: 

Bei der Eulerschen Betrachtungsweise wird das Geschehen in einem Raume verfolgt, der 
unabhangig von den Vorgangen in dem Kontinuum ist. Als raéumliches Bezugssystem dient ein 
raumfestes Koordinatensystem. Der Ablauf des Geschehens in dem Kontinuum selbst wird 
dabei nur implizit erfaBt. 

Die Lagrangesche Betrachtungsweise richtet sich hingegen primar auf den Ablauf des Ge- 
schehens in dem Kontinuum (Kérper). Als raumliches Bezugssystem wird bei folgerichtiger Durch- 
fihrung ein kérperfestes Koordinatensystem benutzt, das die Formanderungen des deformier- 
baren Kontinuums mitmacht. Ein solches Vorgehen setzt freilich voraus, daB die Vorstellung 
von einem absoluten Raum aufgegeben und durch das ,,Relativitatsprinzip des Raumes in der 
klassischen Mechanik“ ersetzt wird!; denn erst dieses Prinzip erlaubt es, alle Raume mit euklidi- 
scher Metrik als grundsatzlich gleichberechtigt anzusehen. 

Natiirlich kann auch die Eulersche Betrachtungsweise sich beliebig bewegter raumlicher Bezugs- 
systeme (starrer und nichtstarrer) bedienen. Ein innerer Zwang dazu besteht aber nicht, weil das 
raumliche Bezugssystem grundsatzlich unabhangig von den Vorgangen im Kontinuum bleibt. 
Wir wollen jedoch in der Folge beliebig bewegte Bezugssysteme fiir die Eulersche Betrachtungs- 
weise mit zulassen. Der Begriff ,,raumfest“* bedeutet dann nur noch: ,,fest in Bezug auf das vom 
Kérper unabhangige raumliche Koordinatensystem”. In der Regel wird man sich dabei auf starre 
Koordinatensysteme beschranken. 

Bei den physikalischen GréBen, die explizit in die Beschreibung eingehen, treffen wir zweck- 
maBig folgende Unterscheidungen: 


1) Dem Wesen nach unterteilen wir in raumgebundene und raumungebundene Gréfen. 
Das Merkmal raumgebundener Gréfen ist es, dafi sie physikalisch in Beziehung zum (kérper- 
unabhangigen) Raum stehen und deshalb physikalisch nicht invariant gegen beliebige Transfor- 
mationen des raumlichen Bezugssystems sind. Dazu gehéren beispielsweise Verschiebungen, 
Geschwindigkeiten und Beschleunigungen gegeniiber dem Raum, aber auch kinetische Energie, 
potentielle Energie der Lage usw. Raumungebundene GréfBen, wie z. B. Spannungen, Form- 
anderungen und thermodynamische Zustandsgréfen, sind dagegen physikalisch invariant gegen 
solche Transformationen. 


2) Der Zuordnung nach unterscheiden wir lokale und substantielle GréBen. Lokale 
GréBen sind den Raumpunkten zugeordnet, und ihre zeitlichen Anderungen werden bei fest- 
gehaltenem Raumpunkt verfolgt. Substantielle GréBen sind dagegen festen Kérperpunkten 
zugeordnet. 


3) Der Darstellung nach haben wir raumbezogene und kérperbezogene Gréfen, je 
nachdem, ob die MaBzahlen der GréBen auf das raumfeste oder das kérperfeste Koordinaten- 
system bezogen sind. 

Zur Erlauterung dieser Begriffe wollen wir in aller Kiirze einige Beispiele betrachten. 


a) Die Geschwindigkeit eines Kérperpunktes gegeniiber dem Raum ist, wenn wir ihren zeitlichen 
Verlauf bei festgehaltenem Kérperpunkt verfolgen, eine raumgebundene, substantielle GréBe, die 
raumbezogen oder kérperbezogen dargestellt werden kann. Betrachten wir jedoch das Ccschwindl “ 
keitsfeld im Raume, so ist der Feldvektor eine raumgebundene, lokale GréBe, deren MaBzahlen 
18 wiederum nach Belieben auf das raumfeste oder das kérperfeste Koordinatensystem beziehen 

énnen. 


* Vgl. hierzu: Th. Lehmann, Schiffstechnik 6 (1959), S. 49. 
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b) Der Spannungstensor in einem Kérperpunkt ist eine raumungebundene, substantielle GréBe, 
wenn wir seine zeitliche Anderung bei festgehaltenem Kérperpunkt verfolgen. Er ist eine lokale 
GréBe, wenn seine Anderung bei festgehaltenem Raumpunkt verfolgt werden soll. In beiden 


Fallen kénnen die MaBzahlen des Tensors sowohl auf das kérperfeste wie auf das raumfeste Ko- 
ordinatensystem bezogen werden. 


ce) Die Anderung raumungebundener, skalarer GréBen kann lokal oder substantiell verfolgt 
werden. Dabei kénnen wir wiederum raumbezogene und k6rperbezogene Darstellung unterscheiden. 


Aus den Definitionen und den Beispielen erkennen wir: 


1) Die Einteilung in raumgebundene und raumungebundene Groen bezieht sich auf inharente 
Kigenschaften, die als gegeben hinzunehmen sind. 


2) Die Unterscheidung lokaler und substantieller GréBen ist davon abhangig, wie der Vor gang 
betrachtet wird. Fiir die Betrachtung eines Zustandes ist hingegen diese Unterscheidung be- 
deutungslos. Deshalb wirkt sich diese Unterscheidung im wesentlichen nur bei Differentiationen 
nach der Zeit und bei Integrationen iiber die Zeit aus. Die Besonderheiten eines Vorganges legen 
haufig eine bestimmte Betrachtungsweise nahe. Kinematische Randbedingungen oder stationare 
bzw. zeitlich periodische Felder im Raume drangen z. B. nach lokaler Betrachtung, wahrend dem 
Kérper zugeordnete dynamische Randbedingungen sowie die Verfolgung von Zustandsanderungen 
eine substantielle Verfolgung der Vorgange herausfordern, besonders wenn die ,,Geschichte“ des 
Vorganges oder Anisotropieeigenschaften der Kérperelemente dabei ¢ine Rolle spielen. 


3) Die Art der Darstellung, ob raumbezogen oder kérperbezogen, unterliegt der Willkiir. Es 
ist jedoch naheliegend, substantiell betrachtete Vorginge kérperbezogen und lokal betrachtete 
Vorgainge raumbezogen darzustellen. Die Zusammenhange werden dadurch durchsichtiger. Der 


physikalische Unterschied zwischen raumgebundenen und raumungebundenen GréfBen spielt aber 
dabei eine Rolle. 


Diese Erérterungen zeigen, daf} der wesentliche Unterschied zwischen Eulerscher und Lagrange- 
scher Betrachtungsweise in der Zuordnung der GréBen zum Raume (Euler) oder zum Kérper 


_ (Lagrange) besteht. Welche Betrachtungsweise angemessen ist, hangt von den Gegebenheiten des 


Vorganges ab. Im allgemeinen werden wir damit rechnen miissen, daB zunachst einige Gréfen 
(z. B. kinematische Randbedingungen oder raumliche Kraftfelder) dem Raume, andere GréSen 
(z. B. dynamische Randbedingungen oder Temperaturfelder) dem Kérper zugeordnet sind. Es 
ergibt sich somit haufig die Aufgabe, lokale GréBen durch substantielle auszudriicken oder, um- 
gekehrt, substantielle durch lokale. Die hier bestehenden Zusammenhinge sollen die folgenden 
Betrachtungen aufzeigen. Was sich dabei ergibt, ist im wesentlichen bekannt'. Gelegentlich zu 
findende Unklarheiten bei der Formulierung von Problemen der Kontinuumsmechanik sowie die 
Feststellung, da die kérperbezogene Darstellung haufig ttbergangen wird, mégen aber vielleicht 
doch eine zusammenfassende Betrachtung rechtfertigen. 


2. Allgemeine Beziehungen zwischen Lagrangescher und Eulerscher Betrachtungsweise. a) Grund- 
formeln. Wir benutzen ein raumfestes Koordinatensystem x‘ und ein kérperfestes Koordinaten- 
system £', Der im allgemeinen nichtlineare Zusammenhang zwischen diesen beiden Koordinaten- 
systemen sei umkehrbar eindeutig und durch 


& = E(k, i) ° x aie, t) (1) 

gegeben. Die Koordinatendifferentiale transformieren sich linear: 
Sete pe OFF al 2 
dé = 55 dt, dx = 5 1 . (2) 


ee 
Die Transformationsmatrizen kénnen als lokale GréBen 


agi Ont. 
ek (G3t)' gen (*"> t) 


— Rational Mechanics 
1], J. G. Oldroyd, Proc. Royal Soc. London 200 (1950), S. 523. — 2. C. Truesdell, J. 
and Analysis, 1 (1952), S.125; 2 (1953), S.593. — 3. J. G. Oldroyd, Proc. Royal Soc. London 245 (1958), 
S. 278. — 4. R. Hill, J. Mechanics and Physics of Solids. 7 (1959), S. 209. 


22 


318 Th. Lehmann: Beschreibung von Vorgangen in der klassischen Kontinuumsmechanik _Ingenieur-Archiv 


oder als substantielle GréBen a 
Dar aS de 
ak (64) > pee (6 #) 


eingefiihrt werden. 
Die Transformation beliebiger Tensoren! gehorcht dem Schema 


en en ax*! Ax°? 


qiite-- ges, Soe Time e 
ky Kgeepore Ox”! ay ae’ aehe $1 Spee 
und 
i. i Ss 85 
ty tgese dx" Ox? og" 0 bao ens 
ak $1 826 06 


ere gels get ax? ax 


Uberstrichene Tensoren sind hierbei auf das kérperfeste, uniiberstrichene Tensoren auf das raum- 
feste Koordinatensystem bezogen. 


Abgekiirzt schreiben wir hierfiir 


sities aes) » Bee \ terse 
Ls Nec (= AT (=F a 235: ° 


Ox) 5/556 

aoe ore oe°R \i ites: 
oe sara Sas, ee 

fe HT (i) As) 

Die zeitabhangigen Transformationen (3) vermitteln den Zusammenhang zwischen raumbezogenen 

und kérperbezogenen GréfBen fiir beliebige Zeiten t. Wir kénnen raumbezogene GréBen als lokale 

und entsprechend kérperbezogene GréBen als substantielle GréBen ansehen, wenn wir nur einen 

bestimmten Zeitpunkt ¢ ins Auge fassen. Die Ermittlung des Zusammenhanges zwischen den 

zeitlichen Anderungen der raumbezogenen und der kérperbezogenen GréBen erfordert hingegen 

weitergehende Uberlegungen, da sich mit der Zeit nicht nur die GréBen selbst, sondern auch die 


Transformationsmatrizen andern. AuGerdem haben wir hierbei substantielle und lokale Anderungen 
zu unterscheiden. 


Rg 


D 


(3) 


~ 
R 


a 


=~ 


b) Allgemeines zur Differentiation nach der Zeit®. Verfolgen wir bei der Differen- 
tiation die zeitliche Anderung einer substantiellen GréBe, d. h. halten wir bei der Ableitung den 
Kérperpunkt fest (€ = konst.), so sprechen wir von substantieller Differentiation. Fiir diese 
ist das Symbol D/dt gebrauchlich. Die Differentiation kann auf kérperbezogene (iberstrichene 
Tensoren) und auf raumbezogene GréBen (uniiberstrichene Tensoren) angewendet werden. 

Wird hingegen bei der Differentiation der Raumpunkt festgehalten (x' = konst.) und somit 
die zeitliche Anderung einer lokalen Gré8e verfolgt, so wollen wir dieses eine lokale Differentiation 
nennen, fiir die wir das Symbol D/dt benutzen. Sie ist ebenfalls auf kérper- und raumbezogene 
GréBen anwendbar. 

Der substantielle wie der lokale Differentialquotient gelten jeweils fiir einen bestimmten Zeit- 
punkt. Wir kénnen den Differentialquotienten selbst deshalb als substantielle oder lokale GréBe 
betrachten und ihn dementsprechend durch kérperbezogene oder raumbezogene GréBen darstellen. 
Wir erhalten somit insgesamt acht verschiedene Ausdriicke, die wir im folgenden Abschnitt be- 


trachten wollen. Zuvor seien nur noch einige allgemeine Beziehungen abgeleitet, die wir fiir diese 
Betrachtungen bendtigen. 


1) Die Geschwindigkeit des Kérperpunktes gegeniiber dem Raum ist durch 


DB cers 
ae (4) 
definiert. Dabei ist v' auf das raumfeste Koordinatensystem bezogen. 
2) Es ist 
DEF oe pe 
ae reeds) 


1 Wir setzen hier voraus, daB alle physikalischen GréGen als Tens defini i Laie : 
Einschrankung der Allg ewanihere. y ensoren definiert sind. Darin liegt keine 


2 ; . : 
a gae sie der Zeit kann auch jeder andere Parameter treten, von dem der Ablauf des Vorganges ein- — 


XXIX. Band 1960 ‘ : PS : 
Th. Lehmann: Beschreibung von Vorgangen in der klassischen Kontinuumsmechanik 319 


wobei v' derselbe Vektor wie v' i : cages 
Pata besogen. r wie v' ist. Seine MaBzahlen sind jedoch auf das kérperfeste K oordinaten- 


oxm 
= py 


D / dxi 


ai (ack 


Hier bezeichnet v'|,, die kovariante Differentiation, die hier im raumfesten Koordinate t 

auszufiihren ist, weil v' raumbezogen ist. In der Einfithrung der kovarianten Faleitun tn (6) 
kommt zum Ausdruck, daB wir die Anderungen von v' auf dieselbe Basis beziehen, auf fe ak 
bezogen ist, namlich auf die Grundvektoren im Punkte x‘, an dem sich der ecachiete Kor : 
punkt zur Zeit ¢ befindet. Das ist folgerichtig, weil hier v' raumbezogen ist. Fiihrt man d ai 
anstelle der kovarianten Ableitung den partiellen Differentialquotienten dv'‘/0x,, in (6) sie es 
verfolgt man die Anderung der MaBzahlen von v', aber nicht unmittelbar die Adenine des Vek. 


- tors selbst. 


: D /0& 
Zur Bestimmung von a (ast) gehen wir davon aus, daB 


Of! Oxm 
xm oer 


Axi gem 
~ OE" Ox 


= 6 


ist. Die substantielle Differentiation des ersten Ausdruckes ergibt 


ax™ D / ag&i ag D /axm\ __ 0 
ag dt arn) ox™ dt (37) = ; 


a apes art i quoee® 
Multiplizieren wir diese Gleichung mit = , so erhalten wir 


o& ox™ D / 0k a o& o& D /dx™ iF 
dx og dt (se) ax™ Oxk dt (sr) 


Daraus kénnen wir nun leicht ablesen, daB 


: D /aéi a) OE 
en (se) = a (7) 
ist. 
4) In ahnlicher Weise finden wir durch lokale Differentiation 
D / a6! _, aem 
& (a) ig a (8) 
und 
D / xt EO 
ad (ae) = 0", oem * (9) 


_kérperfesten Koordinatensystem 


wo. 


Die kovarianten Differentiationen sind hier im kérperfesten Koordinatensystem auszufihren. 


5) Die Relationen 
far OMe ens 20x ‘amrogt ay osm 
a im 9Ek = OTe BEM? ole Bae = Ohm Buk a 


deren Giiltigkeit mit Hilfe der allgemeinen Transformationsformeln (3) leicht nachzuweisen ist, 
zeigen, daB fiir die Transformationsmatrizen die Begriffe substantielle Differentiation und lokale 
Differentiation, so wie sie hier definiert sind, im Ergebnis zusammenfallen. Darin driickt sich aus, 
daB die Transformationsmatrizen ihrer Bedeutung nach stets sowohl dem raumfesten wie dem 
zugehéren. Wenn wir trotzdem im folgenden bei den durch- 
zufiihrenden Differentiationen aucli bei den Transformationsmatrizen die Unterscheidung in 
lokale und substantielle GroBen beibehalten, so geschieht dies, um den jeweiligen Standpunkt zu 


kennzeichnen und die Einheitlichkeit der Betrachtung deutlich werden zu lassen. 
22* 
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c) Substantielle und lokale Differentialquotienten. Fir eine auf das kérperfeste 
Koordinatensystem bezogene GréBe ist die substantielle Differentiation gleichbedeutend mit der 


partiellen Ableitung nach der Zeit (bet &! = konst.): 


Dri ae pee 11 
‘diakee Fel ie 4 (11) 


Wir kénnen diese Differentiation auch nach Newtonscher Schreibart mit einem Punkt (’) bezeichnen, 
da Verwechslungen nicht zu hefiirchten sind: 


(a:i: Ja ae (12) 
Das Gegenstiick dazu ist die lokale Differentiation emer raumbezogenen GréBe 


DB) eerie, OT aiar 
eee Fee j 1 
dt ea Fe Ot bie (13) 


Die partielle Differentiation ist hier bei x‘ = konst. auszufiihren. Da Verwechslungen nicht zu 
befiirchten sind, kénnen wir auch hier die Newtonsche Schreibart anwenden 


Exe 7 eae 
(oe Ja fee. (14) 


Als substantielle Differentiation einer raumbezogenen GréBe definieren wir 
Dimetra iene : ; ae 
an =a a" OG) 

a o °|m dt ( fas aa -kg- m 


Diese Definition des substantiellen Differentialquotienten eines raumbezogenen Tensors fihrt 
wieder auf einen Tensor, weil die zeitliche Anderung folgerichtig auf die Grundvektoren am Ort x! 
bezogen wird, an dem sich der Kérperpunkt zur Zeit t befindet. Der fiir diesen substantiellen 
Differentialquotienten manchmal zu findende Ausdruck 


(15) 


re aie: 


ist hingegen kein Tensor, weil in den zweiten Term hier auch die Anderung der Grundyvektoren 

langs der Bahn des Kérperpunktes eingeht, und weil darum dieser Ausdruck sich auf zwei ver- 

schiedene Raumpunkte bezieht, die um dx' auseinanderliegen (vgl. hierzu auch Abschnitt 2.b). 
Das Gegenstiick zu Gleichung (15) bildet die lokale Differentiation einer auf das kérperfeste 

Koordinatensystem bezogenen GréBe 

a © res a din ge 1 a 


(16) 


Dém Be 
EN ae 


m at oe 


Wollen wir den substantiellen Differentialquotienten einer kérperbezogenen GréBe (Gleichung 
(11)) durch raumbezogene Groen ausdriicken, so haben wir — 


D —i,. a D “Ty: aga ax"B 
=—@ (09 = —_/a@a a , 

d chins d aSi3 I] r ] 

; é zl B (22 ats 


£ B 


zu bilden. Die Ausfiihrung der Differentiationen auf der rechten Seite ergibt 


De OND ys ek eel acia , [ ax°8 Peed aéia ax'B 
aed Oe pea ra) A ap | 8 spas Dera 
B org et |m| LE V5 rn a NOB pdt | EE \ are JAE \ ack 
° y Sees 3 s i i 
Ga Ye ag hs IT ag'o ,[ ax°B wraaFa: ag"A aé'o *B 
oo i ea hel | | ) ( |- a, it [T 2 IT’ ee 
( e° Bim |e \ axe Jag \ asks 2 es Pho i ax’= JAg \ agks 
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Nach Vertauschung einiger Summationsindizes kénnen wir dafiir auch schreiben 


ie Pe 
gia: sas ote F cai) \| = Ger “M2 .Ty Qe 
dt °.k3. lo" ym sas 4 ag oe 
= ral: agia ax'B 
a3 e A\eg 4 a 5 4 : (17) 
beth, eee II ax" IT aes 


In (17) ist der substantielle Differentialquotient einer kérperbezogenen GréSe durch raum- 
bezogene Gréfen ausgedriickt!. Entsprechend kénnen wir den lokalen Differentialquotienten einer 
raumbezogenen GréBe [siehe Gleichung (13)] durch kérperbezogene GréBen ausdriicken. Das 


fiihrt auf 
nein ig eae , Axia a&°B D — ees a ty oeB 
Seer Ce ea ta* 2) ain ayo | Ox é 

a a Ere } = l ] i RSs Sa , ] , i 
a 8 a sate ( of * ]} “8 ax6 dt B NC ad (os B (a 


ae SD ; axia ae'B =e =. Sie 
+a Re Bees = ( : a m a 
eG. dt I (2° IT | meee: v Osge 
Solis UO ea GG] A ya (ice) ty oe*8 
Bee ate IT Cee ES) 


| (2) (25) 


0g” gata \ OE" ax B 
Set eee oe” axio aes 
— a A , 
De, Pm IT (=) T(E}: 
A 0x4 Ta 0& * ] ga \ Ox B 


Nach geeigneter Umordnung kénnen wir dieses zusammenfassen zu 
ae ye=4 gata: 
m er Snes 
_18r|m > Perey 


aor ea , ax’ ag" 
: 18 
me TE 


ago 


Wir kénnen ferner den substantiellen Differentialquotienten einer raumbezogenen GréBe 
(Gleichung (15)) durch kérperbezogene GréBen darstellen. Dann erhalten wir 


Ditage 2D | airs. 7 [Ox * ae*é ) 
ae Age a (2) m(;)| 


a 


Doe once rion Be Be a oD , [oxi aé°B 
= _—a * —— (i aaa — 
dt “3B. IT (=) IT (5) Ps ses dt I (=) IT (=) (19) 


way =r = 10, : —m aT se : axa aé*F 
=|(a)+ 2% tse tests See 0; | I (=) (25) 


SchlieBlich kénnen wir den lokalen Differentialquotienten einer kérperbezogenen GréBe (Glei- 


_ chung (16)) mit Hilfe von raumbezogenen GréBen angeben. Das ergibt 


a B 


hea ee ert eee ag’a , ( ax‘ 
mega ales (Za) 


ns ag'a , ( ax‘8 gD agia\ _, [ ax‘ 
“265 a(S) (B)--3 3a Sr c 


B \ 


Hep! agi \ 1, [ax‘é 
(4 gig See mend Mis 
Sg A ae ea (255) 


a Cima 


nL Te +Mm.Fr 4° my 
Sli ee ee > v* A) a ee tea ale > v 
( a : _____ my"... 8g. aia 


1 Qldroyd und Hill nennen diesen Ausdruck a.a.O. ,,convected derivative’, Truesdell a.a.O. dagegen 


_ ,,absolute time flux“. Die von diesen Autoren angegebenen Ausdriicke sind nur insofern etwas spezieller, 


F 


als sie voraussetzen, daB kérperfestes und raumfestes Koordinatensystem zur betrachteten Zeit t gerade zu- 
sammentallen. 
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d) Einige erginzende Bemerkungen. Die in Abschnitt a) und c) ae ee 
Beziehungen erméglichen es, lokale Gré®en durch substantielle auszudriicken oder sp: San 
und damit die beiden Darstellungsweisen ineinander zu iiberfiihren. Die Ausdriicke assen sic 
noch ein wenig umgestalten, wenn wir den Gradienten des Geschwindigkeitsfeldes additiv in 
seinen symmetrischen und seinen antimetrischen Anteil, also in den Deformator (dj) und in den 
Rotor (w', bzw. ;') zerlegen: 


vl =F (Oe tol) + > mel) (21) 


= di + OK 


und entsprechend 


1 . : 
n= 5 (all +o) + 5 Cole) 
22 
bg =, (22) 
=id} +a; 
In gleicher Weise erhalten wir 
0. 1 ak eee 
m= sO Hel) +h — vl) (23) 
= dj + ie 
und 
fateh ar (24) 
= dy, + Oy « 


Die Umschreibung der Formeln (17) bis (20) mit Hilfe dieser Beziehungen kann hier tibergangen 
werden. Es werden sich spater jedoch noch einige interessante Zusammenhange ergeben. 


Die Ausdriicke (16), (18) und (19), die zur kérperbezogenen Darstellung gehéren, scheinen bisher 
in allgemeinen Betrachtungen iibergangen worden zu sein. Desgleichen ist anscheinend der Aus- 
druck (20), der zu (16) gehért, bisher ohne Beachtung geblieben. Dafiir werden gelegentlich auch 
noch andere Differentialquotienten angegeben. Genannt seien hier: 


1) ,,Intrinsic derivative“!. Hier werden die zeitlichen Anderungen auf ein Koordinatensystem 
bezogen, das an die Kérperpunkte gebunden ist, aber nur die translatorischen Bewegungen der 
Kérperpunkte mitmacht und nicht die Drehungen und Formanderungen des Kérpers. Dieses 
Koordinatensystem ist im allgemeinen nichtholonom. Die Ableitung stimmt formal mit dem 
substantiellen Differentialquotienten einer raumbezogenen GréBe iiberein, wenn das raumfeste 
Koordinatensystem zur Zeit t mit dem kérperfesten zusammenfallt. 


2) ,,.Material derivative’ 2, Die zeitlichen Anderungen werden hier in einem Koordinaten- 
system verfolgt, das die Translationen und die Drehungen der Kérperpunkte mitmacht, aber nicht 
die Formanderungen des Kérpers. Das — im allgemeinen ebenfalls nichtholonome — Koordinaten- 
system behalt also wahrend des Vorganges an jedem Punkte seine urspriingliche Metrik (MaB- 
tensor) bei. Diese Definition entspricht dem Cosseratschen Kontinuum’. 


Vom Standpunkt der phainomenologischen Kontinuumsmechanik erscheinen diese beiden 


Definitionen von Differentialquotienten als recht willkiirlich. Sie werden deshalb hier nicht 
weiter verfolgt. 


e) Integrationen iither die Zeit. Fiir Integrationen iiber die Zeit gelten ahnliche Uber- 
legungen wie fiir Differentiationen nach der Zeit. Wir kénnen hierbei wiederum Integrationen 


1 vel. Oldroyd und Hill a. a. O, 


® Truesdell a, a. O. benutzt diesen Ausdruck abweichend von Oldroyd und Hill fiir die substantielle Differen- 
tiation einer raumbezogenen GroBe. 


* vgl. Truesdell a.a.O. und W. Giinther, Abh. Braunschweig. Wiss. Ges. 10 (1958), S. 195. 
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bei festgehaltenem Korperpunkt (&' = konst.), d.h. substantielle Integrationen und lokale 
_ Integrationen (x‘ = konst.) unterscheiden. AuBerdem kénnen die Integranden kérperbezogen 
oder raumbezogen sein, so da wir in Analogie zur Differentiation vier verschiedene Integrale 


haben. 


Der Unterschied zur Differentiation besteht darin, daB die Integrale einem Zeitintervall, die 
Differentialquotienten hingegen einem Zeitpunkt zugeordnet sind. Substantielle Integrationen 
iiber die Zeit kénnen deshalb nicht durch lokale Integrationen iiber die Zeit ausgedriickt werden!. 
Unberiihrt von diesem Sachverhalt kann natiirlich das Ergebnis der Integration stets wahlweise 


kérperbezogen oder raumbezogen dargestellt werden, wenn man es einem bestimmten Zeitpunkt 
(z. B. der oberen Integrationsgrenze) zuordnet. 


Die Integrationen iiber die Zeit sind wie gewéhnliche Integrale ausfiihrbar, wenn bei substantiel- 
len Integralen der Integrand kérperbezogen und bei lokalen Integralen der Integrand raumbezogen 
ist. In diesem Falle benétigen wir auch fiir die substantielle und die lokale Integration keine 
unterschiedlichen Symbole (in Analogie zur Newtonschen Schreibweise der Differentiation nach 
der Zeit). Eine solche Darstellung der Integranden ist aber auf Grund der in Abschnitt a) und c) 
zusammengestellten Beziehungen stets méglich. Die Ausfiihrung substantieller und lokaler Inte- 
grationen tiber die Zeit bereitet deshalb keine grundsatzlichen Schwierigkeiten. 


3. Beispiele. a) Verschiebung, Geschwindigkeit, Beschleunigung. Bei Einbeziehung 
endlicher Verschiebungen ist die analytische Definition der Verschiebungen an ein raumfestes 
Koordinatensystem gebunden, dessen MaStensor im ganzen Raum konstant ist. Bezeichnen 
wir dieses Koordinatensystem mit X', so kénnen wir die Verschiebungen definieren durch 


UN(EF, t, to) = X(E*, t) — XE", to) . (25) 


Die Transformation auf beliebige andere raumfeste Koordinatensysteme x’ ergibt die allgemeine 
Definition der Verschiebungen 


wl, tte) = [XU 2) — XU, 1] 25, (P40). 26) 


Damit ist die Verschiebung als substantielle, raumbezogene GréBe eingefiihrt. Ihre MaSzahlen 
sind auf die Grundvektoren des Raumpunktes bezogen, an dem sich der Kérperpunkt &* zur Zeit t 
befindet. Die kérperbezogene Darstellung wird durch 


TCE, t to) = ws ts) ae (FA) 7) 


vermittelt. Da x‘ eine eindeutige Funktion von £* und t ist, steht es uns ferner frei, die Verschie- 
bungen auch als lokale GréSen aufzufassen, indem wir 


ub = u(x", t,i,) bzw. ui = ui(x*, t, t) (28) 
_ setzen. 
Die substantielle Differentiation der raumbezogenen Verschiebung ergibt? 
D , D ax! Dis. 
Is SES aici k eeu /ek 
ae ES ts hy) = ek Et ee (é*, t) = 3 ai(E*, 0) a 
== 0° (E", t) - haw. vi(x*, t) 

1 Kinem Ko6rperpunkt ist in einem Zeitintervall eine stetige und dichte Folge von Raumpunkten zugeordnet 
und umgekehrt. Den Zusammenhang zwischen substantieller und lokaler Betrachtungsweise bei Integrationen 
iiber die Zeit vermitteln deshalb Linienintegrale, auf die hier nicht weiter eingegangen werden soll, 

x ist hier nicht mit zu differentiieren, weil wir definitionsgemaf die zeit- 


oxm 
liche Anderung eines Tensors auf dieselbe Basis beziehen wie den Tensor selbst. 


2 Die Transformationsmatrix 
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oder nach (19) auch 


m w"| as at v e $ f (30) 


= ui(ék, t, to) Be [(w) + oF 


wobei wir die rechte Seite nach Belieben als Funktion von &*, t oder x", t auffassen kénnen. Die 
lokale Differentiation der raumbezogenen Verschiebung fiihrt dagegen auf 


= Ue) = (u') (31) 


bzw. 
=|(e)—or# |, +a] (32) 


Ss 


Hier sind die rechten Seiten Funktionen von x* (bzw. é*), t und t). Die Ausdriicke (31) und (32) 
werden im iibrigen nur fiir sehr spezielle Betrachtungen bendétigt, z. B. etwa wenn die Herkunft 
der Kérperpunkte untersucht werden soll, die an einem bestimmten Raumpunkt vorbeiwandern. 
Das Gegenstiick dazu ist der lokale Differentialquotient nach der Zeit t, 


Di 0 ee 33 
dig tke ts bo) = op Wa te Mo) (33) 


der bei der Ermittlung der ,,Streichlinien“ eine Rolle spielt. 


Die zeitlichen Anderungen der kérperbezogenen Verschiebung fiir sich allein zu verfolgen, ist 
physikalisch irrelevant und kann hier iibergangen werden. Hierin kommt zum Ausdruck, daB die 
Verschiebung eine raumgebundene Grife ist, die deshalb auch raumbezogen einzufihren ist, 
wenn ihre zeitlichen Anderungen verfolgt werden sollen. 


Auf eine andere Méglichkeit der raumlichen Zuordnung der Verschiebungen sei noch kurz 
hingewiesen. Wir kénnen — im Gegensatz zu (28), wo u' dem Raumpunkt x* zugeordnet ist, 
den der Kérperpunkt & zur Zeit t einnimmt — die Verschiebung auch dem Raumpunkt x* zu- 
ordnen, an dem der Kérperpunkt zur Zeit t, war. Das fiihrt zu einer anderen analytischen Definition 
der Verschiebung, namlich auf 

xi 


Ox 
Wi (RE, 0) = [Xm(@, 1) — XG] ew ) - (34) 


Die Abhangigkeit von t) steckt hier implizit in x*. Die lokale Differentiation ergibt bei dieser 
Zuordnung der Verschiebung 


" (35) 


: sey iabuer of ever ao Lee : : 

Die Geschwindigkeit v‘ ist hier wie u‘ auf die Grundvektoren des raumfesten Koordinatensystems 
: ok 

im Punkte x" bezogen. 


Der substantielle Differentialquotient verliert bei dieser Definition der Verschiebung seine 
selbstandige Bedeutung; wir kénnen ihn allenfalls formal dem lokalen Differentialquotienten 


gleichsetzen, weil uw’ stets zugleich einem festen Raumpunkt und einem festen Kérperpunkt zu- 
geordnet ist. Der Ubergang zur kérperbezogenen Darstellung wird im iibrigen hier durch eine 
Abbildung und nicht durch eine Transformation vermittelt, weil bei dieser Betrachtungsweise 
einander entsprechende Raum- und Kérperpunkte im allgemeinen nur zur Zeit fj raumlich zusam- 
menfallen. Die dazugehérigen Uberlegungen laufen infolgedessen ein wenig anders als die all- 
gemeinen Betrachtungen in Ziff. 2. Diese Gedankenginge sollen hier jedoch nicht weiter verfolgt 
werden, weil der Begriff der .,Abbildung“ seiner Natur nach mehr fir den Vergleich zweier Zu- 
stinde als fiir die Verfolgung eines Vorganges geeignet erscheint. 
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Die Geschwindigkeit ist durch 
psec anbe 
i a (36) 


detiniert und kann als substantielle GréBe v'(&*, t) oder als lokale Grée vi(x*, 1) betrachtet werden 
[vgl. Gleichung (29)]}.. Die kérperbezogene Darstellung ist durch 


(37) 


-gegeben. Bei der Verfolgung der zeitlichen Anderungen ist jedoch von der raumbezogenen Dar- 
stellung auszugehen, wenn die Betrachtung physikalisch sinnvoll sein soll, weil die Geschwindig- 
keit eine raumgebundene Gréfe ist. Wir erhalten somit folgende Ausdriicke: 


0 = (a) +o", (38) 
= |e) + 21,0] (39) 
und 
ze=(), (40 
= |(e)— wre, + 2," ge 
- “ 
ag 


Die linke Seite der Gleichungen (38) und (39) kénnen wir als Definition der Beschleunigung be- 
trachten 
Daa. D 


Se es, (42) 


bi 
Die Gleichungen (40) und (41) beschreiben dagegen die zeitliche Anderung des Geschwindigkeits- 
feldes. Die genaue, formale Entsprechung der raumbezogenen und der kérperbezogenen Dar- 
stellung tritt im iibrigen nur bei der Verfolgung des Geschwindigkeitsverlaufes auf, und dort auch 
nur, wenn die Geschwindigkeit kontravariant definiert ist. Die Entsprechung ist also keine al!- 
gemeine Higenschaft. 


b) Spannungen. Es liegt im allgemeinen nahe, die Spannungen dem Kérper zuzuordnen 

und sie deshalb kérperbezogen einzufiihren: o‘*, oj, oder o;,. Wir wollen uns hier der gemischt- 

_varianten Schreibweise bedienen, die verschiedene Vorziige hat. Die substantielle Differentiation 
dieses Tensors ergibt 


=, 4 = (&), (43) 


= [(cz) +» 24), — 271, 0° + 21,00] ne oi : (44) 


Bleiben die beim Vorgang eintretenden Zustandsanderungen ohne Riickwirkungen auf den Vor- 
_ gang wie bei inkompressiblen, viskosen oder ideal-plastischen Medien, oder ist der Vorgang stationar 
im Raum, so kann es zweckmaBig sein, die Spannungen den Raumpunkten zuzuordnen und sie 
deshalb raumbezogen darzustellen (0;,). Die lokale Differentiation ergibt dann 


ar% =(ci) - 


=(() oa, + 1,0" — |, o] se ee (46) 
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Wesentlich seltener wird es vorkommen, daB die lokalen Anderungen des kérperbezogenen 
Spannungstensors oder die substantiellen Anderungen des raumbezogenen Spannungstensors 
verfolgt werden miissen. Deshalb kénnen diese Differentialquotienten hier iibergangen werden, 
zumal die hierbei méglichen Darstellungen aus den allgemeinen Betrachtungen in Ziff. 2 c) sofort 
abgelesen werden kénnen. 


Diese kurzen Ausfiihrungen zeigen, daB bei raumungebundenen GréBen raumbezogene und 
kérperbezogene Darstellung grundsatzlich gleichberechtigt nebeneinander stehen. Es hangt vom 
Einzelfall ab, in welcher Weise man solche GréBen am zweckmaBigsten einfiihrt. 


c) Formanderungen. Betrachten wir die Lage des kérperfesten Koordinatensystems im 
Raume zu zwei verschiedenen Zeitpunkten (f, und t), so wird der Ubergang von der Ausgangslage 
(Zeit t,) zur Endlage (Zeit t) durch einen Abbildungstensor vermittelt, der auf das raumfeste Ko- 
ordinatensystem bezogen ist. Diesen Abbildungstensor kénnen wir auf zwei Weisen multiplikativ 
in einen symmetrischen Dehnungstensor und in einen antimetrischen Versor aufspalten, je nach- 
dem, ob wir uns 


1) die Drehung vor der Dehnung oder 
2) die Dehnung vor der Drehung 


ausgefiihrt denken!. Bilden wir das Quadrat des Dehnungstensors und transformieren wir diesen 
quadratischen Dehnungstensor im Falle 1) auf das kérperfeste Koordinatensystem zur Zeit t, so 
erhalten wir bei gemischt-varianter Darstellung 


Gi: = 8" Brk- (47) 


° 
Hierbei ist g‘” der kontravariante MaBtensor des kérperfesten Koordinatensystems zur Zeit t 
und g,, der kovariante zur Zeit t. Ein geeignetes MaB fiir die Formanderungen gewinnen wir, 
wenn wir den Logarithmus der Wurzel? des quadratischen Formanderungstensors bilden: 


& = (Inq) = (In gi. (48) 


Der Gedanke, einen solchen Tensor einzufiihren, geht auf Hencky® zuriick. Seine Vorteile hat 
besonders Richter* herausgestellt. 


Zu formal gleichen GréBen kommen wir, wenn wir bei Interpretation 2) den quadratischen 
Dehnungstensor auf das kérperfeste Koordinatensystem zur Zeit t, transformieren. Dieses Ko- 
ordinatensystem kénnen wir als einen Sonderfall des raumfesten Koordinatensystems ansehen, 


weil seine Lage im Raum unabhangig von t ist, und wir wollen es mit x‘ bezeichnen. Wir erhalten 
in diesem Falle 


"Bk (49) 


und 
O; a Ovi 
é = > (Inq)k. (50) 


Hier sind gi und & den Raumpunkten zugeordnet, die die Kérperpunkte zur Zeit f, einnahmen. 
Der Ubergang zur kérperbezogenen Darstellung dieser GréSen fihrt wieder tiber Abbildungen. 
Diesen Weg wollen wir aber nicht weiter verfolgen, sondern wir gehen hier von dem kérperbezogenen © 


ee Formanderungstensor (48) aus, den wir natiirlich auch raumbezogen darstellen 
onnen: 


i =, Oxt age 
= es Der Oxk : (51) 


1 Einzelheiten dieser Gedankengange s. bei Leh 
formtechnik, Ing.-Arch. 29 (1960) Sila ee 
Logarithmus und Wurzel sind hier als Matrizenfunktionen zu verstehen. Fiir solche Funktionen wollen 
wir die Schreibweise (In q)j» (q~4)i, usw. benutzen. 
Sei Hencky, 2. angew. Math. Mech. 5 (1925), S. 144, 
* H. Richter, Z. angew. Math. Mech. 29 (1949), S. 65. 


Einige Betrachtungen zu den Grundlagen der Um- 


ere 
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Der substantielle Differentialquotient des kérperbezogenen logarithmischen Formanderungs- 
tensors ist 


wa-), (52) 
=|(é) 4 ome, ore + om on] (53) 


Fihren wir die Differentiation auf der rechten Seite von (52) aus, so erhalten wir unter Beriick- 
sichtigung, daB qj, ein symmetrischer Tensor ist, 


Das ergibt, wie man sofort zeigen kann, wenn man mit Hilfe von (17) Dg,,/dt in einem raumfesten 
Koordinatensystem ausdriickt, das zur Zeit t mit dem kérperfesten zusammenfallt, 


De ge oh ht eee 
ae = ZB (Uke + Pel), = (7h + %!) - 64) 


Vergleichen wir dieses Ergebnis mit den Gleichungen (23) und (24), so finden wir, daB 
D =j ze 
Tale Oban di, (55) 


ist, daB also bei gemischt-varianter Darstellung der substantielle Differentialquotient des kérper- 
bezogenen, logarithmischen Formanderungstensors gleich dem kérperbezogenen Deformator des 
Geschwindigkeitsfeldes ist. Wir kénnen aber auch schreiben 


0& xs 


D =v 1 r r 
a ck = = (1), + el’) Ser BER (56) 
gy O8i xs 
S@xr Qgk * (57) 


Fallen insbesondere zur Zeit t das kérperfeste und das raumfeste Koordinatensystem zusammen 
so wird 


also gleich dem raumbezogenen Deformator des Geschwindigkeitsfeldes. Dieser tiberaus einfache 
Zusammenhang ist eine besondere Higenschaft des hier eingefiihrten Formanderungstensors. Er 
gilt im iibrigen nur fiir die gemischt-variante Darstellung von logarithmischem Formanderungs- 
tensor und Deformator des Geschwindigkeitsfeldes; denn die substantielle Differentiation des 
kovarianten, kérperbezogenen Formanderungstensors ergibt z. B. 


Dx jb 
“a bik = Gp Bir ti): 


= Pa 1 = ae D . 
= (#4, + Yr|,) fk + 9 Bir (“, a Yrl,) (59) 
=2di,& +n, 
=d,, (28 + 5) - 
Die substantielle Integration des kovarianten, kérperbezogenen Deformators fiihrt hingegen auf 
t t 
= Lenk. Joe bie 2 
| a ums | Genet = (g4.—sin)- (60) 


ty ty 
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a Fay es = 5 e 2 : Fa 
Der Ausdruck ce (pissing kann nicht unmittelbar als ein kovarianter Formanderungstensor 
ie} 


gedeutet werden. Die MaBtensoren g;, und g;, gehdren zwar zum gleichen kérperfesten Koordi- 
natensystem, aber zu zwei verschiedenen Zeitpunkten. Transformationen auf andere raumliche 
Bezugssysteme sind deshalb nicht ohne weiteres durchfihrbar. Erst dadurch, daB wir ,,per defi- 
nitionem“ diesen Ausdruck einem Koordinatensystem (dem kérperfesten zur Zeit t) oder dem 
kérperfesten zur Zeit t) zuordnen, wird ein Tensor im tiblichen Sinne! daraus, wobei natiirlich je 
nach der Zuordnung zwei verschiedene Tensoren entstehen. Der Zusammenhang dieser Tensoren 
mit dem logarithmischen Formanderungstensor ist in folgender Weise gegeben. Die Uberfihrung 
in die gemischt-variante Form ergibt in dem einen Falle 


Sls i =4(e'8,—%) => (i—«&) (61) 
und im anderen Falle 
58 (Ge — Bn) = 5 (8% —8"'u) => (— GH). (62) 


Die rechte Seite von (61) lat sich als das erste Glied der Reihenentwicklung von 


r 


on 1 ) 
& = 9 (in C .04F 6), 
deuten, die rechte Seite von (62) dagegen als das erste Glied der Reihe von 
a) - (eae Ley ane : 
Ex = (In qh = —=-(In ieee = 0), : 


Der kontravariante, kérperbezogene Deformator hangt mit dem MaBtensor des kérperfesten 
Koordinatensystems durch folgende Beziehung zusammen: 


il ed Dees Fik 

iter Hee (63) 
Das 148t sich wiederum mit Iilfe von (17) zeigen, wenn raumfestes und kérperfestes Koordinaten- 
system zur Zeit t zusammenfallen. An diesen Zusammenhang lassen sich ahnliche Betrachtungen 
wie an Gleichung (60) anschlieBen. Diese kénnen hier iibergangen werden. Wir wollen aber fest- 
halten, daB bei kérperbezogenen Gré8en zeitliche Differentiation (bzw. Integration) und Heben 
bzw. Senken eines Indexes im allg. nicht miteinander vertauschbar sind. Anmerken wollen wir 
ferner, daB bei kérperbezogenen GréBen kovariante Differentiation nach den Kérperkoordinaten 
und zeitliche Differentiation (bzw. Integration) im allgemeinen nicht miteinander vertauscht 
werden kénnen. Beide Erscheinungen sind dadurch bedingt, daB im allgemeinen 


De eee ee 
ay Bik = Gp Bik AO bzw. en lente a Mies 0 


ist [siehe Gleichung (60) bzw. (63)]. Bei raumbezogenen GréBen bleibt hingegen, sofern das raum- 
liche Bezugssystem starr ist und die sonstigen Voraussetzungen gegeben sind, die Vertauschbar- 
keit erhalten. 


Die bei Truesdell? angegebenen Definitionen des Formanderungstensors (Cauchy-Green, Green- 

St. Venant, Alamansi-Hamel) sind raumbezogen eingefiihrt. In der einen Version wird der Form- 

anderungstensor den Raumpunkten zugeordnet, an denen sich die Kérperpunkte zur Zeit t, be- 

fanden, in der anderen Version den Raumpunkten, die die Kérperpunkte zur Zeit t Gnnebitign? 

Wir wollen hier nur auf die zweite Version, die uns im Rahmen dieser Ausfiihrungen néher liegt 
kurz eingehen. In ihr nimmt imit unserer Bezeichnungsweise der Cauchy-Greensche Fornjanderdngs 

tensor (von Truesdell so genannt) die Form | 

° r s 
Cik = Srs ge st (64) 


ax’ ax" 


1 In dieser Form wird der Formanderungstensor, allerdi f ei 
un W. hie Theoretical Elasticity, Oxford 1954, iene eck a8] eget Pama aie 
a.a.O. 
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an. Der kérperbezogene Cauchy-Greensche Formanderungstensor ist dagegen 


= rh ° ag” 0&* ex” ox” ° 


Cc z ik 9? 
ik Srs ax Qa” ac aek Sik (65) 
und infolgedessen ist 
1D ; 
dt Cir — 0 . (66) 


Die substantielle Differentiation des raumbezogenen Tensors ergibt 


D 
Tap thier —(¢ mk UV v”|, at Cin 


m 


vl), (67) 


wie wir sofort durch Anwendung von (19) ableiten kénnen. Dieser Ausdruck zeigt keinen einfachen 
Zusammenhang mit dem Deformator des Geschwindigkeitsfeldes. Darum ist, wie auch (66) zeigt, 
dieser Tensor nicht sehr gliicklich definiert. Uberfiihren wir hingegen Gen gemischt-varianten 
Formanderungstensor von “Alamansi und Hamel, der bei Truesdell in der Form 


ej. = ob —= ¢ (68) 


angegeben ist, in die kérperbezogene Darstellung, so finden wir 


@, = 6; — Bi" Brk > (69) 
also wieder den Ausdruck (62), der das erste Glied der Reihenentwicklung von é, darstellt. 


Da bei den von Truesdell angegebenen Formanderungstensoren kein einfacher Zusammenhang- 
zwischen dem substantiellen Differentialquotienten und dem Deformator des Geschwindigkeits- 
feldes besteht, unterscheidet er die Begriffe ,,rate of deformation“ (entspricht dem Deformator 
des Geschwindigkeitsfeldes) und ,,rate of strain‘ (entspricht dem substantiellen Differential- 
quotienten des Formanderungstensors). Diese Begriffe fallen jedoch in ihrer Bedeutung genau 


- zusammen, sofern nur der Formanderungstensor in geeigneter Weise definiert wird. Insbesondere 
gilt dies fiir den gemischt-varianten logarithmischen Formanderungstensor, dessen weitere Vor- 


— 


zuge (eindeutige additive Aufspaltungsméglichkeit in Volumendilatation und Formanderungs- 
deviator, lineares Superpositionsgesetz, Invarianz bei affiner Abbildung) ihn fiir allgemeine Be- 
trachtungen besonders geeignet machen. 


Die Verfolgung lokaler Anderungen des Formanderungstensors kénnen wir hier iibergehen. 
Sie bereitet auf Grund der vorausgegangenen Betrachtungen keine Schwierigkeiten. Meist wird 
man auch dabei vom kérperbezogenen Formanderungstensor ausgehen, weil dies aus physikalischen 
Griinden niéher liegt. Bei den Formanderungsgeschwindigkeiten kommt es dagegen wieder sehr 
auf den Hinzelfall an, ob die kérperbezogene oder raumbezogene Darstellung vorzuziehen ist. 


AbschlieBend sollen nur noch zwei Beziehungen hervorgehoben werden, die sich mit Hilfe des 
logarithmischen Formanderungstensors ebenfalls besonders einfach formulieren lassen. Fiir die 
Formanderungsarbeit eines Massenelementes mit der Dichte @ erhalten wir 


{2 ‘ai dt = Spe. 3 (70) 


Die Kompatibilitatsbedingungen fiir die ARE eee SE kénnen wir schlieBlich 
in der Form 


Cee), + (ell, — (eel, — (ee). = ° (71) 
schreiben. Das ist gleichbedeutend mit 
dif, + dnl, — dnl 


und lat sich in dieser Form durch Ausschreiben sofort nachpriifen. Es ergeben sich damit be- 


d,| = 0 (72) 


merkenswerte Analogien zur Theorie infinitesimaler Formanderungen!. Fiir den Formanderungs- 


1 vel. z. B. H. Schaefer, Z. angew. Math. Mech. 33 (1953), S. 356. 
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tensor é, selbst sind die Kompatibilitatsbedingungen sehr viel komplizierter. Sie kénnen nicht 
unmittelbar durch Integration der Gleichungen (71) iiber die Zeit gewonnen werden, weil zeitliche 
Integration und kovariante Differentiationen bei kérperbezogenen GréBen nicht miteinander 
vertauschbar sind. 


4. SchluBbemerkungen. Bei der Beschreibung von Vorgangen in einem Kontinuum wird weit- 
gehend von der Eulerschen Betrachtungsweise Gebrauch gemacht. In neuerer Zeit werden aber 
in zunehmendem Mafe auch Fragen aufgegriffen, deren Beantwortung eine Verfolgung des Vor- 
ganges am Kérperelement erforderlich macht. Damit riickt die Lagrangesche Betrachtungsweise 
wieder mehr in das Blickfeld. Die gegenseitigen Beziehungen zwischen diesen beiden Betrach- 
tungsweisen wurden hier zusammengestellt. Dabei sollte gleichzeitig gezeigt werden, daB die 
kérperbezogene Darstellung die Lagrangesche Betrachtungsweise wesentlich vereinfachen kann, 
und daB die Zusammenhange sehr viel durchsichtiger werden, wenn man bei Definitionen und 
Darstellung immer den physikalischen Sachverhalt im Auge behalt. 


(Eingegangen am 30. November 1959.) 
Anschrift des Verfassers: Dozent Dr.-Ing. Th. Lehmann, Hannover, PodbielskistraBe 93 
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Ausweitung des kreiszylindrischen Hohlraumes 
unter 6rtlichem Innendruck 


Von E. Tremmel 


1. Einleitung und Problemstellung. Spannungen und Verschiebungen eines allseitig ins Unend- 
liche erstreckten, elastisch-isotropen Kérpers mit einer durchgehenden, unter konstantem Innen- 
druck p stehenden kreiszylindrischen Offnung vom Querschnittshalbmesser a sind durch die Airy- 
sche Spannungsfunktion 


F=—pdinr (1) 


bestimmt. Fiir die Radialverschiebungen u des Zylindermantels erhalt man hier den einfachen 


Ausdruck 


cep 
u=3 6%) (2) 


in dem G den Schubmodul bezeichnet, der mit dem Elastizitatsmodul E und der Querzahl y durch 
die bekannte Beziehung G = E/2 (1+) verbunden ist. Wirkt der Innendruck dagegen nur 
értlich, etwa iiber die in der Achsrichtung z des Hohlzylinders gemessene Lange 2c (Abb. 1), wie 
z. B. beim Abpressen von Probestrecken in Stollen und Schachten!, so liegt ein rotationssymmetri- 
sches Problem vor, bei dem die von r und z abhingigen Spannungs- und Verschiebungskomponen- 
ten aus einer Loveschen Verschiebungsfunktion abgeleitet werden kénnen. 


Abb. 1. Kreiszylindrischer Hohlraum unter értlichem hydrostatischem Druck; Koordinaten und Bezeichnungen. 


Eine Naherungslésung fiir diesen Fall hat H. M. Westergaard? angegeben, wahrend die exakte 
Lésung, d.h. die den hier vorliegenden Randbedingungen geniigende Lovesche Verschiebungs- 


 funktion von C. J. Tranter? berechnet wurde. In der zitierten Arbeit stellt Tranter ferner die 


allgemeinen Ausdriicke fiir die Spannungen und Verschiebungen auf und bestimmt die GréBe der 
im Symmetriequerschnitt auftretenden Radialverschiebung des Zylindermantels unter Verwendung 


eines numerischen Verfahrens. Den Fall einer auf den zylindrischen Hohlraum wirkenden kon- 


stanten, einseitig ins Unendliche erstreckten rotationssymmetrischen Druckverteilung behandelt 
O. L. Bowie’, ohne jedoch auf die Verschiebungen einzugehen. 

SchlieBlich seien in diesem Zusammenhang noch die Arbeiten von L. Féppl* und A. W. Rankin® 
genannt, in denen das verwandte Problem des unter rotationssymmetrischen Druckverteilungen 


stehenden, unendlich langen Vollzylinders untersucht wird. 


aaa 


Si i 


is ill ~~ 


1 H, Lauffer, Geologie und Bauwesen 25 (1959) S. 114. 

2 H, M. Westergaard, Karman Anniversary Volume, 1941, S. 154. 
3 C, J. Tranter, Quart. Applied Math., 4, (1946) S.298. 

40. L. Bowie, Quart. Applied Math., 5, (1947) S. 100. 

5 L,, Féppl, Drang und Zwang, Bd. 2, S. 168 ff., Miinchen 1928. 

6 A, W. Rankin, J. Appl. Mech. 11, A 77 (1944). 
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Wir werden im folgenden zunachst die Lovesche Verschiebungsfunktion fiir den unter értlichem 
Innendruck stehenden zylindrischen Hohlraum ableiten, wobei wir im wesentlichen auf die oben 
zitierte Veroffentlichung von C. J. Tranter zuriickgreifen kénnen. Bei der Berechnung der Radial- 
verschiebungen des Zylindermantels werden wir uns jedoch nicht auf die numerische Ermittlung 
der Funktionswerte fiir bestimmte ausgezeichnete Querschnitte beschranken, sondern allgemeine 
Beziehungen fiir den Gesamtverlauf entwickeln. Aus einem Grenziibergang wird schlieBlich der 
von einer im Querschnitt z = 0 konzentrierten Druckverteilung bewirkte Verlauf der Radial- 
verschiebungen, d.h. die FinfluBfunktion erhalten, mit der die Ausweitung des Zylindermantels 
fiir beliebige rotationssymmetrische Belastungen berechnet werden kann. 


2. Berechnung der Loveschen Verschiebungsfunktion. Die Lovesche Verschiebungsfunktion @ 
geniigt der Bipotentialgleichung 


AAG =0, (3) 
in der 
preeos ib Pavs) 2 
A= or? t r Or F Oz” (4) 


den Laplaceschen Operator darstellt; Spannungen und Verschiebungen sind mit ® durch die 
nachstehenden Beziehungen verkniipft! (Abb. 1): 


2G 0 aD 
= i ae x P49 — FAL: (5) 

2G 0 1 a®@ 

Ge = recede oes AP aaa (6) 
2G 0 a@ 

Go, = [ir & (2-48 52 | (7) 
2G ) ad 

Tz = 1975 F) [a —) 40 — Fa]. (8) 

und 

See tea 

1 —2 Or dz ’ (9) 


(10) 


1 o@ 
Ww = L295 a ») A® — FB 


Die fiir den vorliegenden Fall in Betracht kommende Lésung der Bipotentialgleichung (3) hat | 
den folgenden Bedingungen zu geniigen: Am Zylindermantel miissen || 


1) die Radialspannungen entgegengesetzt gleich der gegebenen Innendruckverteilung sein: 


P a —e<xz< +¢ 
—o,= fi = d 
‘ ur  r=a_un Le ede A (11) 


2) die Schubspannungen identisch in z verschwinden: 
tee = 0 fue. f= a und).— 60'<..8 =. 06: (12) 
SchlieBlich besteht die Forderung, daB 


3) sdintliche Spannungs- und Verschiebungskomponenten im Unendlichen beschrankt bleiben 
also 


lim o;, 4 00, (13a) 
lim u; 4 0o (13 b) 
x .—> CO 


J 


gelten muB, wobei hier i, 7, k = 1,2 all i in = 
Si eee rience J allgemein anstelle der Koordinaten x, =r und x, = 2 ge- 


1 Siehe etwa C, B. Biezeno-R. Grammel, Technische Dynamik I, S. 116, Berlin 1953. 


3 
ia 
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Zur Bestimmung der Loveschen Verschiebungsfunktion © gehen wir zundchst von der Rand- 
bedingung (11) aus und stellen die am Zylindermantel wirkende Belastung durch ein Fouriersches 
Integral dar, das fiir die durch (11) definierte, in + ¢ unstetige Funktion bekanntlich durch 


oo 


; ; 
6, = p = Ela cos « z dx (14) 


4 
6 
gegeben ist. Wir werden daher auch die gesuchte Lovesche Verschiebungsfunktion in Form eines 
uneigentlichen Integrals mit der Integrationsvariablen « einzufiihren haben. Da o, gemaB (5) 
durch eine ungerade Ableitung nach z mit ® verbunden ist, mu8 © demnach in der Form 


D = J 90) L(, r) sin « 2 da (15) 


angesetzt werden, wobei L(q, r) sin x z eine den Bedingungen (11), (12) und (13) geniigende Lésung 
der Bipotentialgleichung ist; iiber die nur von a abhangige Funktion g(a) wird noch verfiigt werden. 

Fihren wir den Ansatz (15) in (3) ein, so ergibt sich nach Kiirzen von sina z fiir L(x, r) die 
gewohnliche Differentialgleichung 


d? land aL 1 dL 
(ete ete eee, ea 
deren allgemeine Liésung bekanntlich durch 
L = C, In(awr) + Cyr T(ar) + Cs Ko(x r) + C,r Ky(a r) (17) 


gegeben ist. Hierin bedeuten I), I, bzw. Ky, K, modifizierte Besselfunktionen erster bzw. zweiter 
Art (letztere auch Macdonaldsche Funktionen genannt), nullter und erster Ordnung, die fiir die 
Ordnung n allgemein durch 


Ea) = 2 fom cos nt dt (18) 
; 
baw. 
Rix) = J e-* 60! Cof nt dt (19) 


definiert sind. 

Da, wie aus (18) zu ersehen ist, die modifizierten Besselfunktionen erster Art fiir unendliches 
Argument iiber alle Grenzen wachsen, haben wir im Hinblick auf (13) C, = C, = 0 zu setzen. 
Die verbleibenden Konstanten C, und C, werden so bestimmt, da zunachst die Randbedingung 
(12), mit der das in z identische Verschwinden der Schubspannungen am Zylindermantel gefordert 
wird, erfiillt ist. Diese Bedingung fihrt auf eine homogene Gleichung, bei der es nur auf das Ver- 
haltnis der Konstanten C,; und C, ankommt, so daf wir eine der Konstanten, etwa C, = 1 an- 
nehmen, und wenn C, = « A gesetzt wird, endgiiltig 


L=K,(ar) + Aar K,(x r) (20) 


schreiben kénnen und fiir die gesuchte Lovesche Verschiebungsfunktion mit (15) die Darstellung 


@ =f 90) [K,(« r) + Aar K,(x r)] sin(« 2) do (21) 


erhalten. 

- Wir berechnen nun die im weiteren bendtigten partiellen Differentialquotienten der Loveschen 
Verschiebungsfunktion (21). Mit den aus dem zweiten Klammerterm von (16) rekursiv oder auch 
unmittelbar aus (19) berechneten Ableitungen der Funktionen Kj(a r) und K,(«r) folgen die nach- 
stehenden Ausdriicke: 

lo 2) 


aa | x) « [Ky(a r) + A ar Ky(ar)] sin (x 2) da, (22) 


aD fa 
ra =— | Pla) 
6 


Ce Kix rhe [4 ore a K,(«1)| sin (« 2) do , (23) 


23 


~ 
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= = + | 2) « L(x r) cos (x 2) da, (24) 
Ss = — | p(x) a? L(x r) sin (« 2) dx , (25) 
AD =? f 96 a? A K,(« r) sin (« z) da. (26) 


Wird (25) und (26) in (8) eingefiihrt und die Ableitung nach r gebildet, so erhalt man den Ver- 
lauf der Schubspannungen allgemein zu 


aes | a? gx) {A [a r Ko(o r) — 2 (1 — 9) Ky(«7)] + Ky(or)} sin (a2) dx. (27) 


Fir r — a miissen die Schubspannungen gemaf der Randbedingung (12), in z identisch verschwin- 


den. Die Konstante A ist demnach so zu bestimmen, daf der in (27) innerhalb der geschwungenen .- | 


Klammer stehende Ausdruck fiir r = a Null wird, woraus sich 


a K,(« a) 
= 9 —») K,(@ a) —4 a Ky(a a) 


A(x, a) (28) 
ergibt. Mit (28) ist die Lovesche Verschiebungsfunktion (21) bis auf die von r und z unabhangige 
Funktion g(x) bestimmt. 

Wir berechnen nun die Radialspannungen o, unter Verwendung von (5). Nach einigen Zwischen- 
rechnungen folgt 


le. 2) 


ee [eee P — KEgns 7 Kyaw l= 2») Kyler) or Ky) 
j 
— Kol r) + = Ky(« n}| cos (x 2) dx. (29) 


Dieser Ausdruck muB fiir r = a in (14) titbergehen; wir haben daher die Funktion g(x) so zu wahlen, 
daB mit r =a die mit — 2 G/1— 2,» und 2 p/z multiplizierten Integranden von (29) bezw. (14) 
fiir alle positiven Werte von « iibereinstimmen. Da g(x) fiir feste Werte von « als Konstante zu 
betrachten ist, erhalt man aus dem Vergleich der Koeffizienten von (29) und (14) 


pa pa 2 (1 — v) K,(« a) —a a K,(a a) sin &c 
06) = OP ee Lo) K3(a a) + (« a? [Kia a) Kaa] 


(30) 


Die Einfihrung von (30) in (21) liefert schlieBlich die gesuchte Lovesche Verschiebungsfunktion 
zu 


[e.2) 


P= A—2NGS | Ey RIR Tw aPIRIR KITA 
0 


x {[2 (L—7) Ky(a a) —a a Ky(a a)] Ky(« r) +07 K,(« a) K,(x 7)} oe a On (xz)dx. (31) 
Mit der Festlegung der Loveschen Verschiebungsfunktion ist das Problem im wesentlichen 
gelést; man hatte das uneigentliche Integral (31) auszuwerten und damit @ als eine von r und z 


abhangige Funktion gewonnen, aus der die Spannungs- und VerschiebungsgréBen gemaB (5) 
bis (8) bezw. (9) und (10) berechnet werden kénnen. 


Praktisch kommt dieser Weg wegen der verwickelten Form des Integranden, die eine geschlos- 
sene Integration nicht zulaBt, kaum in Betracht. Man wird daher die zur Bestimmung der Span- 
nungs- und VerschiebungsgréBen erforderlichen Differentiationen am Integranden vornehmen, 
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d. h. die Reihenfolge von Integration und Differentiation vertauschen und die so erhaltenen un- 
eigentlichen Integrale 


CO 
y(r, 2) = f Si (x, r, 2) da (32) 
und 
u(r, 2) = f U,(o, r, 2) do , (33) 
ny 


deren Konvergenz in jedem Fall zu tiberpriifen ist, fiir spezielle Werte von r und z mit Hilfe ge- 
eigneter Naherungsverfahren auswerten. 


3. Radiale Deformation des Zylindermantels. Wir beschrinken uns hier auf die Berechnung der 
Radialverschiebungen u(a, z) des Zylindermantels; in ahnlicher Weise kénnen auch die am Zylinder- 
mantel auftretenden Langsverschiebungen w(a, z) bestimmt werden. 

Fir die Radialverschiebungen ergibt sich aus (21) mit (9) 


u(r, 2) = +75 | a ola) [Ky(a 7) + Aa r Kola r)] cos (« 2) do, (34) 


0 


_ woraus nach Einfiihrung von g(x) und A(x, a) gemaB (30) bezw. (28) fiir u das uneigentliche Integral 


pa ep 2 (1 — v) K,(a a) — a a K,(4 a) sin (4 c) 
CON a Sry | 2 (1 —») KX(a ee (a a)? [K2(a = — Kaa] 4 
0 


K,(a a) 
2(1— 1) K,(a a) —a a K,(a a) 


x Kio r) 4 a1 Kj(«r)| cos (x z) dx (35) 
folgt, das, wie mit Hilfe der asymptotischen Entwicklungen der Funktionen Ky und K, im weiteren 
bewiesen werden soll, fur alle Werte r und z konvergent ist. 

Fiir die am Zylindermantel auftretenden Radialverschiebungen erhalt man speziell aus (35) 
mit rT — a 


r KX%o a) sin(4 6) 
ula, 2) = oa 2(1—») | 2(1 — v) K2(a a) + x a? [K2(aa)— K3(aa)} 
0 


cos(~%2z)dx. (36) 


- Wird anstelle von « die dimensionslose Variable 4 eingefiihrt, die durch 


A106: (37) 
bestimmt sei, und durch K3(« a) = K7(A) gekiirzt, so geht (36) tiber in 
co 
1 sin — A . 
u(a, z) = a 2(1—») : KR) cos aA dj. (38) 
Soy) tA } — Ray 


0 


SchlieBlich liefert die Zerlegung des Produktes sin < A cos < vi 


ee ana R Ee eres} 
sin ieee 
pa 2(1—») a aa (39) 
ula,2)= 56 _ Ki(A)] 2 
20) +F IF 


0 


Behandeln wir den Lastbereich, in dem z < ¢ bleibt, getrennt von dem unbelasteten, in welchem 
z> cist, so gilt fir z< c 


oe) le. @) 


. €+2 pe ets 
sin A sin D 
pa2(1—») g da si = (40) 
ua, 2) = 2a ii 7 K3)| 4 
|) ope)? feat) 
0 


0 
23* 
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und fiir z > c 
Z--¢ ae Bs: 
i A sin 
pa2(l—») ea di a dA . (Al) 
eee Ki) K3() | z 


Tt amr, A 
ae = & via eee 
2(1—»)+# I Kx) 2(1—»)+ | K2() 
0 
Wird in (40), also fiir z << c, die Summe der beiden uneigentlichen Integrale mit U bezeichnet 
und der Ubergang zum unendlichen Lastbereich vollzogen, so muB, wie aus dem Vergleich von (40) 
mit (2) hervorgeht, 
Ss 42 
c aa Z : (1 = 9) ( ) 
sein. Fiir z = c miissen die durch (40) und (41) gegebenen Funktionen stetig ineinander tibergehen, 
d. h. die zweiten Integrale in (40) und (41) verschwinden. 
Die in (40) und (41) auftretenden Integrale 
sind von der Form 


co 


sin pA dA 
a f same ei 
(Ka) 
wobei die mit z veranderlichen Faktoren y(z) 


== AP rie Oe auf Werte y(z) = 0 beschrankt 
+ (¢ —2)/a 
bleiben. 

Da die exakte Berechnung der Integrale (43), 
worauf schon hingewiesen wurde, nicht méglich 
ist und die numerische Auswertung, die fiir jede 
Abszisse 2, also fiir jeden Wert y(z) gesondert 
durchzufiihren ware, einen erheblichen Rechenauf- 
wand erfordert, so werden wir den unendlichen In- 
tegrationsbereich in n Teilbereiche aufspalten, in- 
nerhalb welcher die im Nenner stehende Funktion 


fa) =#)1 — (5) | (44) 


(Abb. 2) durch einfach zu integrierende Naherungsfunktionen ersetzt werden kann. 
Bezeichnen wir die dem k-ten Teilbereich mit den Grenzen w,_1 < 1< @, zugeordnete Nahe- 
rungsfunktion mit f,(A), so soll (43) demnach durch die Summe der Teilintegrale gemaB 


Fe, S) k=n ke 
sin pA dA . sin pa dA 
lcceorson ze ceors of (15) 


> Of 1 


A<<1 
fa) tlAj2¥ 


Abb, 2. Die Funktion f(A) = 2? [ 1— ( za) | ‘ 


ausgedriickt werden, wobei man sich im allgemeinen auf die Aufspaltung in vier bis sechs Teil- 
bereiche beschranken kann. 

Um Aufschlu8 tiber die zur Darstellung der Funktion f(A) in Betracht kommenden Naherungs- 
funktionen zu gewinnen, soll zunachst das Verhalten der modifizierten Besselfunktionen zweiter 
Art untersucht werden. Aus der Integraldarstellung (19) ist zu ersehen, daB die Funktionen K,(x) 


fiir reelles, positives Argument gleichfalls positiv und reell sind, und da® ferner mit J == x fiir — | 


0<A<w 
K,(A) 
| Kp) ~ s) 
sein mu. Mit den in der Umgebung des Nullpunktes fiir p = 0 und p = 1 aus den Reihenent- 
wicklungen folgenden Naherungen! 


Ka) == In A (47) 
La (48) 


+ Jahnke-Emde, Tafeln héherer Funktionen, S. 135, Leipzig 1948. 
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worin Iny = C= 0,577... die Eulersche Konstante bedeutet, folgt weiter 


-  K,(A) =. 
are AOE (49) 


Damit aber wird 


_ [Ke(A)\? 
] f = 0 
ee (Ra) po 
von héherer Kleinheitsordnung. Fir Argumente 4 <1 ist der Verlauf der Funktion f(A) 
K,(A)\2 ; ; 
=}? f —(Ra| daher im wesentlichen durch die Parabel f,(A) = 2 bestimmt (Abb. 2). 


Zur Erweiterung des Giiltigkeitsbereiches der Approximation werden wir f(A) besser durch ein 
Polynom der Form 


f(A) =2# + B44 + Bo AS +e. (51) 


annahern, dessen Koeffizienten 8; unter Verwendung der aus Tafeln berechneten, exakten Werte 
von f(A) aus dem Minimum der Fehlerquadratsumme zu ermitteln sind. Um umstindliche, auf 
komplexe Funktionen fithrende Partialbruchzerlegungen zu vermeiden, wird der in (43) enthaltene 
Bruch gemaf 
1 Besse 2) ge (4) Fi) 
i=) 72d —9)| 20-9) SiGe ed 


in eine rasch konvergierende Reihe entwickelt. Die obere Grenze «, des Intervalls, in dem f(A) 
durch f,(A) angenahert werden soll, ist daher so anzunehmen, daB 


F(@:), — f,(@1) 
2(1— 7) 2(1—») «K ! 


wird. Wahlen wir etwa w, = 0,5, dann wird mit vy = 0,2 und f(0,5) = 0,172 


f(@ ) 215 
AED as 0,108 . 


Der maximale Fehler der nach dem dritten Gliede abgebrochenen Entwicklung (52) betragt 


_ dann etwa 1% des richtigen Wertes 1/[2 (1 — v) + f(A)], wahrend der Fehler des Integrals wesent- 


eye 


eB 


lich geringer ist. : 

In ihrem weiteren Verlauf nahert sich die Funktion f(A), wie aus Abb. 2 zu ersehen ist, asympto- 
tisch einer Geraden; schon fiir kleinere Argumente ist f(A) nur schwach gekriimmt, so daB von 
A =, = 0,5 an eine stiickweise Annaherung durch lineare Funktionen méglich ist. Wir setzen 
hier also 


f(A) =wA+x', (53) 


wobei die Konstanten yu und x’ wieder wie oben aus der kleinsten Fehlerquadratsumme zu berechnen 
sind. Damit folgt weiter 


1 1 


: . "(54 
20 —) +f) — 20—») Ful +H’ A 
Wird die rechte Seite von (54) unter Hinfihrung der abkiirzenden Bezeichnung 
ee ae (55) 
7) 
umgeformt, so erhalt man in 
} 1 1 (56) 


ua) 2(1—») +A) 


einen einfachen Naherungsausdruck, mit dem die Integration (43) geschlossen durchgefithrt werden 
kann. 

Im starker gekriimmten Bereich kleinerer Argumente wird die Funktion f(A) nur innerhalb 
eines kiirzeren Intervalls durch eine Gerade anzunahern sein; fiir gréBere Argumente nimmt die 


’ Kriimmung rasch ab, so daB hier die Intervallange vergroBert werden kann. Wir ersetzen dem- 


zufolge die Funktion f(A) etwa im Intervall 0,5 = a, <> w; = 3,0 durch die Naherungsgerade f, 
und im Intervall 3,0 = w, < w3 = 15,0 durch die Naherungsgerade f5. 
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Um AufschluB iiber das Verhalten der Funktion f(A) im Unendlichen zu erhalten, gehen wir von 
der bekannten asymptotischen Darstellung der modifizierten Besselfunktionen zweiter Art aus; 
fiir groBe Argumente x und die Ordnungen 1 und 2 gilt? 

le oe pe tie bs iG. 19 ones 
K(x) =| Ix ° 1 (@aye Gap alot ay 312 ue |, ee, 
sen Po: Epes 1 30521 
Ix ° fore (ap atat’ @xp 3128 —-. be) 


K(x) =)/ 


Beide Funktionen fallen also mit wachsendem Argument monoton ab und verschwinden im Un- 
endlichen starker als die Exponentialfunktion. 


Bilden wir mit x = / den Ausdruck /? | 1 — (ea) , indem wir die Reihen nach den Gliedern 
Z 
mit A~2 abbrechen, so folgt bis auf Glieder héherer Kleinheitsordnung 
: Kayes pone 9 
yn bees ps = Sey ad Be 
fara SI 72|1 (xa) |=sa By Pate (59) 
als asymptotische Darstellung der Funktion f(A). Fiir die Grenzwerte der in (53) eingefiihrten 
Konstanten gilt 
limp = +1, (60) 
Aco 
lim x’ =— 1, (61) 
4-00 
daher 
limzx = 1—2yr, (62) 
A> co 


und (56) geht iiber in 
1 1 


: Ki\1. 1— iv eA. 
20 +f (EO) 


(63) 


Mit (59) kann nun auch die noch offene Frage nach der Konvergenz des Integrals (43) bzw. (36) 
erledigt werden. Der Integrand ist jedenfalls stetig; fiir groBe Werte A gilt, wenn zur Abkiirzung 
sin pd 1 ee 
z=») +f) zB) 
gesetzt und f(A) gema® der asymptotischen Darstellung (59) ausgedriickt wird, mit 2 > bS>1 


sin pA 


GQ—2n)A+2 ee 


lg(p, a)| = 


wobei M eine positive Konstante ist und damit 


co oo 


ee ciet dae 
[ewau<[=a5=6 
A A 


wird. Da aber ¢ fiir A(e) > b unabhangig von y beliebig klein gemacht werden kann, ist die gleich- 
maBige Konvergenz von (43) und (36) nach bekannten Satzen der Analysis gesichert. 


Von etwa A = 15,0 an wird der bei der Naherung (63) auftretende Fehler unerheblich, so daB 
wir im Bereich 3 = 15,0<A< wy, oo f(A) durch die in (59) definierte Asymptote, die wir hier 
mit f,(A) bezeichnen wollen, ersetzen kénnen. 


oe die durch Einfithrung der Naherungsfunktionen f(A) erhaltenen Briiche (52) baw. (56) 
mit ~ dA multipliziert, so folgen mit k = 1 


41 ra» | A ! se SE GT ee he (64) 
0 


1G. N. Watson, A treatise on the theory of Bessel functions S. 122, Cambridge 1922. 


XXIX. Band 1960 E. Tremmel: Ausweitung des kreiszylindr. Hohlraumes unter ortlichem Innendruck 339 


mmaomit fs=) 2.3. aa 
Of, 


l 
j= | sin pA di 


| ee (65) 


OL—1 


zwei ‘lypen von Integralen, aus denen die Lésung von (43) gemaB (45) aufzubauen ist. 
Das dem ersten Bereich 0 << 4< w, zugeordnete Integral (64) formen wir durch Einfihrung 
von f(A) aus (51) um und ordnen nach sieigenden Potenzen von /: 


ae I ie sinyA[, # 1 44 
Ii aa | H] L 2(1—») F (say A eo dA. (66) 
i) 
Ausmultiplizieren auf der rechten Seite und Trennung der Integrale liefert 


Sei me sin wd 1 ee 
dy 2(1—») l| 7 ar ey [Asinya di 
0 


0 


WO, 


1 By 
+ (cay ai) [# ein dd -—-- ; (67) 


é 


Durch Auswertung der Integrale und Einsetzen der Grenzen ergibt sich schlieBlich 
1 ? 1 Leite 
Jy =a (wy =) Pp (sin yp ©, — Y @; Cos y @;) 


1 ; 1 
+ (gaye ad aa) ge 8 OF —2) sin y os —y 0 (of —6) c08 p olf. (68) 


Dabei ist unter der mit Si bezeichneten Funktion der durch 


x 


i sin t Sy tae tt y2h} 
oe | aoe Ob )-GE ”) 
k=1 


0 


definierte Integralsinus zu verstehen. Fir kleine Argumente gilt die aus der Reihendarstellung 


folgende Naherung x < 1 


Sx, (70) 
wahrend Si (x) fiir groBe Argumente xS>1 durch die asymptotische Darstellung . 
. “13 i 1 
Si (x) = = —* Hi, (x) —“* Hx) (71) 
mit 
2! 4! 
LING Sc 17 te er mea (72) 
und. 
See! 
EG) oe (73) 


angendhert werden kann. Damit folgt weiter 


lim Si(x) ==. 
In diesem Zusammenhang definieren wir ferner den bei der Auswertung von (65) auftretenden 
Integralcosinus durch 


x 


oo cs 2k 
Gifs) =— | dt =Iny— | ee ap 
k=1 


t 
x 0 


1 Jahnke-Emde, a. a. O. S. 3, dort in komplexer Form. 
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mit den Naherungen fiir x < 1 


; 1] 
Ci (x) = —In oar (75) 
daher 
lim Ci (x) = — oo (76) 
x—>0 
und? fir «> 1 
Gi (#) =“S* Hy) — > Hl). (77) 
demnach 
lim Ci (x) =0. (78) 


C = Iny entspricht hierbei der in (47) aufgetretenen Eulerschen Konstanten. 
Mit Hilfe von (70) und den Entwicklungen fiir die Kreisfunktionen erhalten wir fir kleine 
Argumente y < 1 aus (68) 


1 y® w3 1 y w3 1 B wy? w?\ y w? 
aE rah y sreazer ke Women yomere Ts hem er gc Gee iz) | (79) 


und stellen fest, da J, im wesentlichen durch das erste Glied von (68) bestimmt ist und mit y—> 0 
verschwindet. 
Fiir y— oo folgt aus (68) und (71) wegen 


lim H,(x) = lim H,(x) = 1, (80) 
1 1 
Me oa) 5 (81) 


Zur Auswertung der Integrale der Form (65), in denen f(A) durch eine lineare Funktion an- 
genahert ist, setzen wir 


OK eae? 
ee sinyA dA 1 cos pA 
Digatrs ee Aig th  ([ 284) @. e 
0 c 


%—1 , —i 


driicken also das gegebene Integral durch ein Doppelintegral aus, das bei vertauschter Integrations- 
folge 
cos py +2 cos wy 
x | ([Stiaa=s [ore dy (83) 


Of; —1 0 0 Op | 


leicht zu berechnen ist. Fiir das innere Integral erhalt man nach Einsetzen der Grenzen w,__; und O;, 


I cos pA 1 : ; 
it | xa di — <6 {cos p x [Ci p(x + o,) — Ci p(x + ox—1)] 


on | 


+ sin yp x [Si p(x + w,) — Si p(x + ax_i)]} (84) 


Wird die rechte Seite von (84) gemaB (83) nach y integriert, so ergibt sich unter Beachtung der 
Grenzen 0 und p mit (82) 


ae 
J, = a {sin yp [Ci p(x + @,) — Ci ple Eon) 


— cos p x [Si p(x + w,) — Si p(x + w@.—1)] + Sip wo, — Sip Wk —1} + (85) 
Fir Argumente y <1 folgt mit (70) und (75) aus (85) 


ee x + WO, y x+o 1 
hts oe AE apd ca me In a,8, Oe — oe —1) aay (ok of.) (86) 


1 Jahnke-Emde a. a-O. S 3. 


XXIX. Band 1960 EK. Tremmel: Ausweitung des kreiszylindr. Hohlraumes unter értlichem Innendruck 341] 


wahrend man fiir groBe Argumente mit (71) und (77) aus (85) die asymptotische Darstellung 


: ta 1 
y>l I= sey [Oey Oe (pay + a) — crv Or—a/ 1 is 1 ) (87) 


a+ oy ero t wy 


erhalt; fiir y— oo verschwindet das Integral daher. 


Fir den letzten von w,_1 = w,_, bis ins Unendliche erstreckten Teilbereich ergibt sich mit 
On—1<A< @, > 0Q unter Beachtung von (60), (71) und (78) 


1 : ; J 
J, = = = sin y x Ci p(x + w,_1) — cosy x eae Si p(x + on) + +—Si y on} , (88) 


worin x den durch (62) bestimmten Wert 1 — 2 y annimmt. 


Durch Summation der mit (68), (85) und (88) fiir die einzelnen Teilintervalle berechneten 
Integrale gemaB (45) erhalt man schlieBlich das endgiiltige, iiber den unendlichen Bereich 0 <= co 
erstreckte Integral 


J=Jj,4+5J,+---+J,.- 


Wir werden die aus der Integration gewonnenen Funktionen jedoch, nach Teilbereichen getrennt, 
unmittelbar in (40) und (41) einfiihren; hierbei ist y fiir z< c durch die GréBen (z + c)/a und 
(c — z)/a und fiir z > c durch (z + c)/a und (z—c)/a zu ersetzen. 


Mit y = 0 also z =c verschwinden, wie aus (68), (85) und (88) hervorgeht, samtliche Teil- 
integrale; die Funktionen sind daher, worauf oben bereits hingewiesen wurde, stetig. Ebenso 
verschwinden die Teilintegrale (85) und (88) fiir unendliches Argument, wahrend das aus der 
Integration innerhalb des ersten Teilbereiches folgende Integral fiir y—>oo zunachst auf den 
durch (81) gegebenen endlichen Wert fiihrt und die endgiiltigen Grenzwerte hier erst aus der nach- 
stehenden Zusammenfassung der Teilintegrale gemaB (40) und (41) folgen. Fir das Teilintervall 
0<A< a, gilt Folgendes: Die Kinfithrung von (z + c)/a und (-- c — z)/a in (68) liefert fur z< ¢ 


pa2(1—») 1 : z+te : c—Zz. il a \2 , z+te 
CAGE ey epee saa Sie Zot Sie, a ae paste a 


@,a Fae a \2 , Cz @,a Cz 
ae cos W sin @W Sor cos W 
z+e ear +( =) wea c— ee 


Zz 


a (3 al 5 ee 7) [Peet {3 (ot E - =| = 2) sin @, | 
=e ee (o2( : “J — 6) cos Oe = a) aE (x) (3 (otf 2) simon —* 


== 6), —— (ot (J —9) cos on * >| ss (89) 


und fiir z > ¢ 


op ¢.2(1—7) U : Zz oe —ce 1 a : ate 
u,(a, z) = eae eer ae Si Or Si o— idle eee ot 
@, 4 Pie a c Bat Wy ia 
> ese cos eee sin @, Ps COs @, — | 


1 B a \4 af2te® _ i mie 
eo ea (3 (0% . 2) sin @1 — 
we i ; (0% ( 1 As 6 £08 04 - : = — (22) <3(0t( 5) —2) sn 1 ss 
SG zZ—c\% \ peo a0 
oak Do (o%( - Soca al 

Wachst der Lasthereich c iiber alle Grenzen, dann verschwinden in (89), also fiir c > z, alle 
Glieder auBer 


Ce 


ES AAO FS 


Si WM, = a ? 


342 E. Tremmel: Ausweitung des kreiszylindr. Hohlraumes unter 6rtlichem Innendruck _Ingenieur-Archiv 


. . Ce 
und wir erhalten mit (81) wegen lim (si W1~ = A Si Oye )=2 
%—> OO 
; pa 2(1—y») rn _ pa 91 
jimm= 36a T=) 26° ae 


SchlieBlich folgt aus (90) fiir gegen Unendlich anwachsende Werte von z 
liny y= 0., (92) 
Z—>00 
d.h. also, der aus der Integration iiber den ersten Teilbereich folgende Verschiebungsanteil ver- 
schwindet im Unendlichen. ; 
Die aus der linearen Approximation von f(A) in den mittleren Integrationsbereichen w,z_1 << A 
< ©, herriihrenden Verschiebungsanteile u, erhalten wir durch Einsetzen der Argumente (z + c)/a 
und (c — z)/a baw. (z + c)/a und (z — c)/a in (85), namlich fiir z << ¢ 


= y 3 a 
te = 5c — ) fsin Cie + oy) 2 — Ci + on—1) >| 
+ sin x —— Ci Ge + on) S* — Cie + ona) FA] 
Sewn Si (e + on) -E*— Si fe + 1) 4] 
— 08 Si Ge + on) S*— Sie + oi) 
+ Sie, E* — Sioa E* + Si oy F* — Si oni 4} (93) 
und fiir z > c 
wp ee oe {sin v5 Cie + oy) TE — Gi +o) 
— sin x2 —* Cie + om) =F — Cie + ona) 4] 
— cos x2 tS Si fe +o) FS — SiGe + on) 
-+ cos x -— Si Ge + 4) FE — SiGe + om) 7] 
+ Si 5% — Si ona ES — Si T+ Sion I (94) 
Fiir den unendlichen Lastbereich verschwindet (93): 
lin, =0, (95) 
ebenso gilt ; 
dime y7==:0.. (96) 


Lassen wir in (85) und (86) die obere Grenze w, = w,—> 00 anwachsen, so erhalten wir den 
vom n-ten Integrationsbereich herriihrenden Verschiebungsanteil u,,. ‘et 
Mit den durch (60) und (63) bestimmten Grenzwerten der Konstanten yu bzw. x folgt, wenn | 
wir die untere Grenze mit w,_, bezeichnen, fiir z << c 
— par(i—v) 1 [eee 
Una oe ae as ede 


z 


— Ci(L—29 + @,—1) = 


c 
a 


Seal 2 ye = C1 Ow eee 


a a 


a 


— cos (1—2¥)=** fr —S10—2y +05) 2) 


a a 


cos (kD yy ae [z— S129 +o) S| 


a a 


+a — (Sion 2 aie: on") (97) 
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und fiir z > c 


ie ae ae sah 
U, 6 ee eae [= sin (1 — 2) 4 CEL 2 + oy) 


Pei 2p) (d= 9 y tp, 1) 


a a 
2 Ges (le Spy 8 : eo (toy on) | 
a 
+ cos (1 — 2 v) =—* poe (l—=29 42 on) 4] 
—(si a — Si oi} (98) 


Aus (97) und (98) erhalt man offensichtlich gleichfalls die oben fiir (93) und (94) abgeleiteten Grenz- 
werte. 


Die endgiiltigen Radialverschiebungen folgen nun aus der Uberlagerung der einzelnen Teil- 


integrale, und zwar gilt fiir die innerhalb des Lastbereiches liegenden Querschnitte gemaB (40) 
falls z < c ist, 


2 


k=n 


ee Sea | (99) 


a 


und fir die auBerhalb des Lasthereiches gema4B (41), falls z > c ist, 


een es om 


ae 
[=] 
bo 
ce 
ai 
z= 
Mi 
3 


we 2) — 2 
k=l 


_ wobei die aus der Integration iiber die einzelnen Teilintervalle @, 1 <— @, erhaltenen Funktionen 


(89), (93), (97) bzw. (90), (94) und (98), unter Weglassung des herausgehobenen Faktors > — go!) 


anstelle der GréBen U einzusetzen sind. He 4 


Die durch (99) und (100) dargestellte Funktion u(a, z) ist, wie bereits festgestellt wurde, stetig 
und verschwindet im Unendlichen. Mit wachsendem Lastbereich 2¢ strebt sie wegen (91) und 
(95) — wie es sein muf8 — dem aus der Airyschen Spannungsfunktion fiir den ebenen Verzerrungs- 
zustand abgeleiteten Wert (2) zu: ; 


lim u(a, z) = ao" 


€ > CO 


4, Radialverschiebungen zufolge einer in z = 0 konzentrierten Druckverteilung. Durch einen 


- Grenziibergang kénnen wir schlieBlich auch die von einer im Querschnitt 2 = 0 konzentrierten 


Druckverteilung P am Zylindermantel hervorgerufenen Radialverschiebungen up (z) erhalten. 
Mit der auf die Langeneinheit des Umfanges bezogenen Intensitat P kann 


P 
pee ta (101) 


gesetzt werden. La8t man nun bei festgehaltenem Wert von P die Lange c gegen Null gehen, 
dann folgt aus (99), d. h. mit z > 
ee 


. Pa2(l—v wo 1 zte Z—¢ 
up = lim 52°C De [Pel : )—u.( : ) (102) 
k 


= 


und fir c = 0 innerhalb des Summenzeichens die unbestimmte Form 0/0. Die Anwendung der 


de l’Hospitalschen Regel liefert 
k=n 298 
up(z) = lim Pa J ue()+ tee “|. (103) 
=i 


a 
| Me 3G TU a 
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worin die Striche hier und im folgenden Ableitungen nach (2 + c)/a bzw. (2 —)/a anzeigen. 
Vollziehen wir nun den Grenziibergang, dann folgt nach Kiirzen von a schlieBlich 


k=n 
tp) =o Se SU; (z): (104) 
k=1 


Wir bilden nun die den einzelnen Integrationsbereichen zugeordneten Funktionen U;, und 
erhalten so fiir den ersten Integrationsbereich 0 <— @, aus (90) 


Pd Qa as z i 3 2 ‘ ‘ 2 
Ue saw sin Wy — ee E 4 — cos 05 +( 7) —2) sin oy 5 
1 By a \4 z\3 e 
4 Gener lac 4 — 24 0] cos @y —— 


=e < ((o 4)— 12 (1 a + 24) sin ral ; (105) 


Die Funktion verschwindet offenbar mit z—>0o; fiir z—> 0 strebt sie, wie die folgende fiir Argu- 
mente z/a<1 geltende, nach dem quadratischen Glied (z/a)* abgebrochene Reihenentwicklung 
zeigt, einem endlichen Grenzwert zu: 


,[ % iJ if z\2 l w3 3 2\2 
(5) =aa =| 6 (os 3) ee ies fro a; i) 
| 1 B we 1 z \2 Zz 
. (ra el 5 f= 3 (1 zh (z <1). (106) 


Fiir die Integrationsbereiche, in denen f(A) durch Gerade angenadhert wird, erhalten wir die 
Differentialquotienten U; unter Beachtung von (83) unmittelbar aus (84) zu 


, il 
U,, & ss feos ea 
a Mu a 


: Zz 
sin % —_ 
+ a 


Cie + or) = — Ci + ons) a 


Si ( + «) =—Si% + o4—1) z| (107) 


Auch diese Funktion verschwindet, wie es sein mu8, wegen (71) und (77) fiir z—>+0o. Ebenso 
strebt sie fiir z—> 0 einem endlichen Grenzwert zu; die fiir kleine Werte z/a giiltige Reihenent- 
wicklung fiihrt hier auf 


Pf Veen Oo, + # 1/z\2 ae 
vi(=)= : [mn ttt & (on, = — (1 — 9)? + 28 In EE (5 <1) (108) 


und damit auf den Grenzwert 


UX(0) = =. Ga ose (109) 


@O,-1 + % 


Fiir den letzten Integrationsbereich, in dem @, = ,—> 00 gilt, folgt aus (107) mit (71) und (74) 
sowie x = 1 — 2 und w = | die Funktion © 
U,, =— eos x = Ci (@,—1 + %) —— sin % <. ee Si (@,—1 + x) al ; (110) 
die fiir z/a—> oo verschwindet und fiir z/a—> 0 logarithmisch unendlich wird, wie aus der Reihen- 
darstellung hervorgeht: 
Palas z 1 2 
Un(Z)=—In Jy atx) +2 3% +3(3) 


a 


x“? In —¥ (@,—1 +2) 


4 Ona $% 


D) (®,—1 — 3 x) 


G < 1), (111) 


Wir haben nun die durch (105), (107) und (110) bzw. (106), (108) und (111) dargestellten 
Ableitungen U;, in (104) einzufiihren und erhalten damit den gesuchten Verlauf der von einem 
in 2 = 0 konzentrierten Innendruck P bewirkten Radialverschiebungen Up(z). 
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Zur Berechnung der von beliebig tiber z verteilten rotationssymmetrischen Innendriicken 
hervorgerufenen Radialverschiebungen kann (104) mit P=1 als EinfluBfunktion verwendet 
werden. 

Da in den einzelnen Teilfunktionen, aus denen Up aufgebaut ist, als Variable nur z/a = z auf- 
tritt, lassen sich die fiir eine bestimmte Querzahl » und a = 1 berechneten F unktionswerte unter 
einfacher Mafstabanderung fiir beliebige Halbmesser a der zylindrischen Offnung verwenden, 
wobei die fiir z berechneten Funktionswerte iiber den Abszissen z — az aufzutragen sind. 


Abb. 3. Radialverschiebungen zufolge einer beliebigen rotationssymmetrischen Innendruckverteilung: 


u(z) = j Uo(z — c) p(c) do + Tae —2) p(s) dc. 


Ist der mit p(¢) belastete Bereich durch die Querschnitte c, und c, > c, begrenzt (Abb. 3), 
so gilt, wenn wir die Kinflu8funktion abkiirzend mit 


up(2) = uo(2) 


bezeichnen, fiir einen innerhalb des Lastbereiches liegenden Punkt z mit cy << z< cy 


us) = f ule —6) ple) dt + f uglé—2) p(6) ae (112) 
und fiir einen auBerhalb des fires liegenden Punkt mit c, < ¢ < z 
uz) = f gle —6) pl) do (113) 
and anslog fiir z< c, < 6; : : 
u(z) = f “up(£ — 2) plo) a . (114) 


Mit Hilfe der EinfluBfunktion u,(z,¢) lassen sich ferner allgemeine Schliisse iber den Verlauf 
der von beliebigen Belastungen p(¢) bewirkten Radialverschiebungen u(z) des Zylindermantels 
ziehen. Aus der quellenmaBigen Darstellung der innerhalb des Lastbereiches c,< 2 < cy, auf:- 
tretenden Verschiebungen gem4B (112) 


u(a, 2) = f uz —C) p() do + f uo(o —2) p(t) ac, 

in der die Belastungsfunktion p(¢) positiv und samt ihrer ersten Ableitung im offenen Bereich 
c, <&< cy stetig und beschrankt sei, ist zundachst zu ersehen, daB u,(z—C) bzw. u9(C —2) 
fiir Derechaindended Argument nur logarithmisch unendlich werden kann, da andernfalls die 
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beiden positiven Integrale divergent wiirden. Auf diese, fiir die Greensche Funktion in zweidimensio- 
nalen Bereichen kennzeichnende Singularitat fiihrt auch unser Naherungsverfahren, und zwar 
gilt fiir verschwindendes Argument d. h. fiir €—> z mit (104) und (111) 


1 at 2) leet. (115) 
a 


NG a 2G cA 


Der zweidimensionale Bereich ist im vorliegenden Falle durch den Zylindermantel gegeben, wobei 
zufolge der Rotationssymmetrie anstelle der punktformigen hier die linienférmige Quelle tritt. 
Wir berechnen nun die Ableitung von (112), die bei Beachtung der singularen Stellen z des Inte- 
granden unter Einfiihrung der kleinen GréBe ¢ aus einem Grenziibergang gewonnen wird, 


MKC a ears " 

u'(z) = lim | f uglz —€) p(2) do + f uo(S —z) plc) a] 
Be On| er zte 

und erhalten durch Differentiation unter dem Integralzeichen und Anwendung der Leibnizschen 

Regel zunachst 


C flere “1 Cy 


f ER ple) dt + woe) pe—e) + f MSE pe) te — ude) e+ 0) 


u’(z) = lim 


e—>0 ye 


Bei der partiellen Integration fallen die mit e—> 0 gegen Unendlich strebenden Glieder uo(e) heraus 
und wir erhalten nach Durchfiihrung des Grenziiberganges 


w (2) = wl — 4) Plea) + f ugle —£) PC) A — tle —2) ples) + J uolE—2) PO) AE. (116) 


Der durch (116) gegebene Wert des Differentialquotienten bleibt fiir alle innerhalb des offenen 
Bereiches c, < z < cy liegenden Werte z endlich. Fallt z aber mit einer der beiden Lastgrenzen 
zusammen, dann kénnen die Differentialquotienten nur beschrankt bleiben, wenn die Belastungs- 
funktion p(¢) dort verschwindet. 

Bei Lastfunktionen, die an den Grenzen c,, (m = 1, 2) unstetig auf Null abfallen, mu8 der dort 


auftretende Differentialquotient der Funktion u(z) daher singular werden. Die logarithmische _ | 


Singularitaét von w’(z) — es handelt sich hierbei offenbar um die Tangente im Nullpunkt der Funk- 
tion t In t mit t = |z — c,,|/a — wird im allgemeinen schwach sein. In der graphischen Darstellung 
des Gesamtverlaufes der Radialverschiebungen tritt sie daher, zumindest bei kleineren MaBstaben, 
kaum in Erscheinung. 


Die fiir die Radialverschiebungen erhaltene Funktion u(a, z) geniigt, wie schon betont wurde, 
obwohl sie auf Grund einer Approximation des Integranden (39) berechnet wurde, saimtlichen 
Grenzbedingungen. Dies beruht darauf, daB der Funktionsverlauf im wesentlichen durch das 
Verhalten des Integranden in der Nahe des Nullpunktes und im Unendlichen bestimmt ist; in 
unserer Berechnung wird dies durch die aus (50) folgende Naherung (51) bzw. durch die Ver- 
wendung der asymptotischen Entwicklung (59) der Funktion 


K,\2 
fa) =? i a (R) | 
beriicksichtigt. 


Kine absolute Gewahr fiir die Giite des Naherungsverfahrens bietet die notwendige Erfiillung 
der Grenzbedingungen nicht. Will man die numerische Genauigkeit erhéhen, so laBt sich das 
am einfachsten durch die weitere Unterteilung des Intervalles @, = 0,5 <A< ws = 3,00 erreichen. 
Eine andere Méglichkeit der exakteren Erfassung des Integranden, wie sie etwa in der Annaherung 
der Funktion f(A) im Bereich starkerer Kriimmung durch quadratische oder noch besser kubische 
Parabeln bestiinde, bringt eine bedeutende Erhéhung des Rechenaufwandes, da die Integrale, 
die dann anstelle von (82) treten, durch Partialbruchzerlegung auszuwerten waren. Die hierbei 
auftretenden komplexen Funktionen Si (£ + i), Ci(E-+ in) sind tabuliert!, so daB dieser 
Weg in besonderen Fallen gangbar ware. 

SchlieBlich_kénnte die asymptotische Naherung der Funktion f(A) bereits von @,; = 0,5 an 
verwendet und die fiir gréBere Argumente 2 gegen Null gehende Restfunktion numerisch aus- 


1M. Mashiko, Tables of generalized exponential-, sine-, and cosine integrals, Numerical computation 
Bureau, Tokio 1953 
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gewertet, oder durch eine einfachere Funktion approximiert und geschlossen integriert werden. 
Das erstere wiirde allerdings den Verzicht auf einen allgemeinen Uberblick iiber die Funktion 
u(a, z) bedeuten, da die Auswertung fiir spezielle Werte von z durchzufiihren ware. Im zweiten 


Fall erhalt man bei einer halbwegs genauen Annaherung der Restfunktion vielgliedrige Ausdriicke, 
deren numerische Auswertung auBerst miihsam ist. 


Die Erweiterung des Giltigkeitsbereiches der asymptotischen Naherung bis A = 0 fihrt ohne 
Beriicksichtigung der Restfunktion — und das gilt ebenso fiir die Berechnung der Radialverschie- 
bungen des Vollzylinders — auf unrichtige Ergebnisse. 


Die theoretisch in Betracht kommende Anwendung des Residuensatzes bringt schlieBlich in 
der Auflésung von n-> oo Paaren transzendenter Gleichungen, auf die die Bestimmung der kom- 
plexen Nullstellen im Nenner des Integranden fihrt, in der Untersuchung der Konvergenz von 


co 
> Res usw. eine Fiille weiterer, bisher noch kaum behandelter Probleme. 
I 


5. Zahlenbeispiel. Zum AbschluB seien noch die Ergebnisse einer numerischen Berechnung 
mitgeteilt. Der Querschnitt der zylindrischen Offnung habe einen Halbmesser von a = 1,41 m, 
die halbe Lange des Lastbereiches sei durch c = 0,875 m gegeben. Das Material habe die Querzahl 


C=0875m. C=G875m 


ON iN N 
e § 
, 
) 
tte 3 £0 : 
E 
SS Makstab 
9 G2 04.06.08 yom 
e SS he a (a ee | 
aS q 
pa 
Geb da 
fadialverschiebungen pe 
des Zylindermantels : 


(Elastizitatstheoretische Berechnung 
unter Annahme eines isotropen | 
Materials, Querzahl »=G2) 


Lc ke 4e 5 
Pa ee eee 30 


Abb. 4. Zylindrische Offnung im unendlichen Raum unter értlichem hydrostatischem Druck; Zahlenbeispiel. 


by = 0,2, wahrend fiir den Schubmodul G und die Belastungsintensitat p keine speziellen Annahmen 
getroffen werden. Die Radialverschiebungen u werden fiir die Querschnitte 


1 2 4 é “Ds 
3 C, C, ae ? 3 


Oe Cc, Dee, 3c 


berechnet. 


Der unendliche Integrationsbereich wurde in vier Teilintervalle w,—1 <A< «, aufgespalten, 


in denen die Funktion 
eidal | Paat KO) | 
. \K,(A) 
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mit hinreichender Genauigkeit durch die Naherungsfunktionen f,() ersetzt werden kann (Tabelle 1) 


Tabelle 1 
Pee EE A a et 
0= a <A< m= 0,5 | fd) = Bt +2 
O,<A<a,= 3,0 | fo(A) = Hah + Hs 
@,<A< w, = 15,0 | F(A) Wg hae 
W,<A< > 0 lets) == 2 ered 


Zur Bestimmung der unbekannten Koeffizienten f,, 44, und “;, sind vorerst die exakten Werte 
der Funktion : Ki) p 

)—(Ra)| 
fiir eine Reihe von Argumenten zu berechnen. In den Tafelwerken findet man entweder die Funk- 


tionen K,(x) und K,(x) selbst oder die Funktionen e* K(x) und e* K,(x)*, oder schlieBlich die mit 
den modifizierten Besselfunktionen durch die Beziehungen 


fQ) =? 


i H(i x) = = Ky(x) 
und : 
— HP(ix) = — Kya) 


verbundenen Hankelschen Funktionen? tabuliert vor. 


Mit den aus dem Minimum der Fehlerquadratsummen erhaltenen Koeffizienten folgen die in 
Tabelle 2 zusammengestellten Naherungsfunktionen; im ersten Bereich wurde f(A) durch ein 
Polynom vierten Grades approximiert, also zur Vereinfachung 6, = i; = ...0 gesetzt. 


Die aus diesen Fehlern folgenden Abweichungen der Integrale liegen weit innerhalb der bei elastizi- 
tatstheoretischen Untersuchungen im allgemeinen erforderlichen Genauigkeitsgrenzen. 


Fiir die numerische Auswertung der Ausdriicke (93), (94), (97) und (98) reichen die etwa in — 
dem Tabellenwerk von Jahnke-Emde fiir Si(x) und Ci (x) enthaltenen Tafeln nicht aus, da die 
Argumentschritte von x = 5,0 an zu gro sind. Die Auswertung der genaueren Interpolations- 
formeln 


Si (e+ ds) Sin + de 822 + GF (oop x 884) 4 


: : 
Ci (x + Ax) = Cix + Ax oo ei (sin x | Set « 

ist, ebenso wie die Berechnung der Funktionen aus den oben angegebenen asymptotischen Ent- 

wicklungen, ziemlich zeitraubend. Die vorliegende Berechnung wurde unter Verwendung der 

in den Vereinigten Staaten vom National Bureau of Standards? herausgegebenen Tafeln durch- 

gefihrt. 


Die Integration der Naherungsfunktionen ergab die in Tabelle 3 angegebenen Werte u,. 


Man sieht, das der vom letzten, unendlichen Integrationsbereich gelieferte Beitrag zu den Ver- 
schiebungswerten klein gegeniiber den anderen ist. In Abb. 4 ist der Verlauf der Verschiebungen 
dargestellt; die in +c auftretende Singularitat der Ableitung wurde bereits erértert. 


Beim Aneinanderreihen mehrerer Lastbereiche der Lange 2c miissen die Koeffizienten von |] 


p 4/2 Ga gegen x gehen: so erhalt man z. B. fiir die Mittelordinate einer zum Ursprung symmetri- 
schen Lastverteilung der Lange 6 c aus der Uberlagerung der fiir die Belastungslange 2 ¢ berechne- 
ten Werte 


. Ug(0) = uz,(0) + 2 uy,(2¢) , 
1G. N. Watson, a. a. O. S. 699 ff. 
2 Jahnke-Emde, a. a. O. S. 232 ff. 
3 Nat. Bur. of Stand., Tables of sine-, cosine- and exponentialsintegrals, vol. I. and II., New York, 1940; 


Nat. Bur. of Stand., Tables of sine- and cosine-integrals (fiir Argumente von 10,0 bis 100,0) New York 1954; 
Nat. Bur. of Stand., Tables of exponential-integrals for Complex Argument, New York 1951. 
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woraus mit den Zahlenwerten der Tabelle 3 


a 
ite (0) == (2,248 + 2 - 0,346) A — = 2,94 ie 


folgt. Fiir den Gré®twert der Radialverschiebungen zufolge einer in bezug auf den Querschnitt 
z = +c symmetrischen Lastaufteilung der Lange 8c ergibt sich allgemein 
Ug max — 2 (uz,(¢) *E Uz (3 c)) 


und mit den Zahlenwerten 


PS - 238 (4d 


= 2 (1,363 + 0,159) 2 * a 


Ug c¢ max 


Da der Faktor von p a/2 G a den Wert z erst mit lim n c erreichen kann, lassen sich die vorstehenden 


n—> oo 
Beziehungen zur iiberschlagigen Kontrolle der numerischen Rechnung beniitzen. 


6. Zusammenfassung. Fiir den Fall eines unter értlichem Innendruck stehenden zylindrischen 
Hohlraumes im unendlichen, elastisch-isotropen Kérper, wurde die Lovesche Verschiebungs- 
funktion bestimmt, wobei hier auf eine Arbeit von C. J. Tranter zuriickgegriffen werden konnte. 
Die Ausweitung des Zylindermantels wurde sowohl fiir diese Belastung und, um die KinfluB- 
funktion der Radialverschiebungen zu erhalten, auch fiir die in einem Querschnitt konzentrierte 
Druckverteilung allgemein berechnet. Die Auswertung der uneigentlichen Integrale erfolgt durch 
Aufspaltung des unendlichen Integrationsbereiches und Ersatz des Integranden durch geeignete 
Naherungsfunktionen, deren Abweichungen von den wahren Werten in der Umgebung des Null- 
punktes und im Unendlichen verschwinden. 


Das hier geschilderte Verfahren kénnte in analoger Weise zur Bestimmung der Radialver- 
schiebungen des unendlich langen Vollzylinders und gegebenenfalls auch des unendlich langen 
zylindrischen Rohres angewandt werden. - 


(Eingegangen am 3. Dezember 1959.) 


Anschrift des Verfassers: Dozent Dipl.-Ing. Dr. Erwin Tremmel, Salzburg (Osterieich), Schwarzstr. 31. 
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_ Berechnung der Druckverteilung von ebenen Schaufelgittern 
mit stark gewdélbten dicken Profilen bei inkompressibler Stromung * 


Von W. Richter 


1. Einleitung. Eine der Grundaufgaben bei der Berechnung der Stromung in Turbomaschinen 
ist die Bestimmung der Strémung im Schaufelgitter. Hier soll dazu die Strémung im axialen 
Schaufelgitter, das bei der Abwicklung eines koaxialen Schnittes durch das Schaufelrad entsteht, 
berechnet werden. Gelést wird die sog. zweite Hauptaufgabe (Direct Problem): Bestimmung der 
reibungsfreien, inkompressiblen Strémung in einem gegebenen Gitter. Als Lésung ergeben sich 
Geschwindigkeits- und Druckverteilung langs der Schaufelkontur, sowie Schaufelkraft und Ab- 
strémrichtung in Abhangigkeit von der Zustrémrichtung. Aufbauend auf diesen Ergebnissen 
kénnen der Einflu8 der Reibung und damit Strémungsverluste und Wirkungsgrad nach dem 
Berechnungsverfahren von L. Speidel und N. Scholz’ ermittelt werden. Ebenso kann das hier 
beschriebene Verfahren nach einer Arbeit von H. Schlichting und E. G. Feindt? auf die Berechnung 
der kompressiblen Gitter-Strémung erweitert werden. 

Hine grundsatzliche Liésung dieser Aufgabe, die jedoch in der rechnerischen Durchfiihrung sehr 
miihsam ist, liegt bereits vor: Nach W. Traupel® kann man das gegebene Gitter durch mehrfache 
konforme Abbildung in eine kreisahnliche Kontur iiberfiihren. Fiir diese wird dann der Potential- 
verlauf iterativ bestimmt und auf das Gitter zuriick iibertragen. Dabei erhalt man die gesuchten 
Lésungen nach einem numerisch-graphischen Verfahren und mit einem Rechenaufwand, der fiir 
ein Gitter bei fiinf verschiedenen Zustrémwinkeln etwa 200 Stunden betragt (vgl. G. Hubert*). 
Nach H. Schlichting® wird das Strémungsfeld aus einer Singularitatenbelegung auf der Profilsehne 
aufgebaut, wobei das gegebene Gitter durch geschlossene Stromlinien dargestellt wird. Zusammen 
mit der iiberlagerten Translationsstrémung ergeben die von den Belegungen induzierten Ge- 
schwindigkeiten die gesuchte Geschwindigkeitsverteilung langs der Kontur. Bei W. Isay® wird 
das Strémungsfeld mit einer Zirkulationsverteilung langs der Kontur erzeugt. 

Von diesen bekannten Lésungen erfordert das Verfahren von H. Schlichting den weitaus gering- 
sten Berechnungsaufwand. Jedoch ist es auf Profile mit maBiger Wélbung und mafiger Dicke 
beschrankt. Ziel der vorliegenden Weiterentwicklung des Berechnungsverfahrens ist es, diese 
Beschrankungen aufzuheben, ohne den Berechnungsaufwand allzu sehr zu erhéhen. Hierbei soll 
ein potentialtheoretisches Berechnungsverfahren fiir Schaufelgitter mit beliebigen Profilen ent- 
wickelt werden, das mit einfachen Hilfsmitteln (Tischrechenmaschine und Rechenschieber) und 
mit einem ertraglichen Aufwand durchfiihrbar ist’. 

Mit dem nachstehend beschrisbenen Verfahren wird der von H. Schlichting eingeschlagene Weg 
fortgefiihrt. Wie dort werden Singularitatenverteilungen vorgegeben, die in Anlehnung an die Wirbel- 
verteilungen von W. Birnbaum® ausgewahlt wurden. Sie werden auf einer gewélbten Skelettlinie im 
Inneren des Profiles angeordnet. Die induzierten Geschwindigkeiten werden zum Teil wie bei A. Betz® 


* Gekiirzte Fassung der von der Fakultat fiir Maschinenwesen der Technischen Hochschule Braunschweig 
genehmigten gleichnamigen Dissertation 1959. Berichter Prof. Dr. H. Schlichting; Mitberichter Prof. Dr. 
H. Petermann. 

1, Speidel und N. Scholz, Untersuchungen iiber die Strémungsverluste in ebenen Schaufelgittern, VDI- 
Forsch.-Heft 464, Diisseldorf 1958. — Eine kurze Zusammenfassung davon ist gegeben in H. Schlichting, 
Anwendung der Grenzschichttheorie auf Strémungsprobleme der Turbomaschinen, Siemens-Zeitschrift, 
33 (1959), S. 429—438, vgl. auch: Journ. Basic Eng. ASME (1959), S. 543—551. 

2 H, Schlichting und E.G. Feindt, Forsch. Ing. Wes. 24 (1958), S. 19—28, vgl. auch Z. angew. Math. 
Phys. 9b, (1958) S. 274284. 

3 W. Traupel, Die Berechnung der Potentialstromung durch Schaufelgitter, Sulzer Techn. Rundsch. Nr. 1 
(1945) und Nr. 2 (1948). 

4G. Hubert, Abh. Braunschw. Wiss. Ges. 9 (1957), S. 106—114. 

| 5 H. Schlichting, Berechnung der reibungslosen inkompressiblen Stromung fiir ein vorgegebenes ebenes 
Schaufelgitter, VDI-Forsch.-Heft 447, Diisseldorf 1955. 

6 W. Isay, Z. angew. Math. Mech. 33 (1953), S. 397409 und Z. angew. Math. Mech. 37 (1957), S. 321—335. 

7 An der Aerodynamischen Versuchsanstalt Gottingen ist zur gleichen Zeit der Ansatz von W, Isay weiter 
entwickelt worden. Dabei konnte der Rechenaufwand durch Einsatz elektronischer Rechenanlagen auf sehr 
kurze Rechenzeiten herabgesetzt werden, vgl. E. Martensen, Arch. Rational Mech. Anal. 3 (1959) S. 235. 

8 W. Birnbaum, Z. angew. Math. Mech. 3 (1923), S. 290—297. 

9 A, Betz, Ing. Arch. 2 (1931), S. 359—371. 
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mit Diagrammen, zum Teil wie bei H. Schlichting mit Tabellen berechnet. Zirkulations- und Quell- 
starken ergeben sich auch hier aus einem linearen Gleichungssystem. Die so erzielte Lésung wird 
ahnlich wie bei dem Verfahren von H. Kriiger! abschlieBend verbessert. Auf diese Weise ist es 
gelungen, den Rechenaufwand fiir die Ermittlung der Druckverteilung eines Gitters bei finf 
verschiedenen Zustrémwinkeln auf etwa 50 Stunden herabzusetzen. 


2. Grundlagen des Berechnungsverfahrens. a) Schaufelkraft und Umlenkung der Stré- 
mung im Gitter. Das Strémungsfeld der vorgegebenen unendlichen Schaufelreihe ergibt sich 
durch Uberlagerung einer Translationsstrémung W,, mit gewissen Wirbel- und Quellsenkenver- 
teilungen in jeder Schaufel. Diese Verteilungen werden weiter unten naher angegeben. Die Schaufel- 
kraft wird dabei aus der Gleichung von Kutta-Joukowsky 
erhalten. Hiernach ergibt sich nach Abb. 1 fir den 
Auftrieb A einer Schaufel 

A = ob We Tie. (1) 


Dabei bedeutet 0 die Dichte des strémenden Mediums, b 
die Breite der Schaufel (senkrecht zur Strémungsebene), 
I’ die Gesamtzirkulation einer Schaufel und W,, die sog. 
Translationsgeschwindigkeit, die nach Abb. 1 gleich 
dem vektoriellen Mittel der Zustrémgeschwindigkeit Wy 
und der Abstrémgeschwindigkeit Wy, ist. Diese resul- 
tierende Kraft A steht senkrecht auf W,,. Sie mige 
nach Abb. 1 in die Umfangskomponente P, = A sin §,, 
und in die Axialkomponente P, = A cos, zerlegt 


¥ | . werden. Fiihrt man dabei noch die Durchsatzgeschwin- 


——> 
AS § 


cS 


: / 
Gittertront 


digkeit U = W,, sin f,, ein, so ergibt sich fiir die durch 
Py = c,5-U8bl und > --P, = ¢ 4 U2b 19) 


Abb. 1. Schaufelkrifte, Zu- und Abstrémung des Gitters. 
A Auftrieb einer Schaufel mit den Komponenten Py in 
Umfangs- und Px in Achsrichtung; [’ Gesamtzirkulation 


einer Schaufel; Wy Zustrémgeschwindigkeit weit vor, W{] definierten Kraftbeiwerte 


Abstrémgeschwindigkeit weit hinter dem Gitter; U Durch- 
satzgeschwindigkeit (Komponente von Wy und Wy] senk- ee 2r d = 2r 
recht zur Gitterfront); Woo iiberlagerte Translationsge- Cy et: un CE DG ctg B oo °* (3) 


schwindigkeit; ['/2t und —[/2t weit vor bzw. hinter 
dem Gitter induzierte Geschwindigkeiten der Wirbel- s x 
reihe I’; By und Bry Zu- baw. Abstrémwinkel; Boo Winkel | Dabei wurde als Bezugslange die Lange / der Schaufel- 


von Wo gegen die Gitterfront; t Teilung; | Schaufeltiefe. : Re ° ° Ps . 
pee ents f means F Senaweens sehne eingefiihrt. Der Winkel f,, ist der Zustrémwinkel 
der iiberlagerten Translationsstrémung W.,. 


Die Umlenkung 6, — f,, ist der Unterschied der Strémungsrichtungen weit vor und weit hinter 
dem Gitter. Dort wirkt die Wirbelreihe J’ (Abb. 1) wie eine kontinuierliche konstante Wirbel- 
verteilung y = J’/t langs der Gitterachse. Sie induziert weit vor dem Gitter die Geschwindigkeit 
v,=+y/2=+T1/2t und weit hinter dem Gitter die Geschwindigkeit v; =—y/2 =—J/2t 
parallel zur Gitterachse. Damit ergibt sich nach Abb. 1 fiir den Zustrémwinkel B, bzw. den 
Abstrémwinkel £,,: ¢ 


ctg B, = ctgB, —I/2tU und ctgB, =ctgB.+T/2tU =ctgB, 4+ Tu. (4) 


Fiihrt man noch die dimensionslose Gesamtzirkulation 


Dae! BPs Le 
ein, so gilt nach (3): ©) 
¢, =4l", c, =41* (I* It + ctg Bs), (6) 
ctg Py =ctgh, +2I* It, tg B, = ctgB, +I* It. (7) 


Somit ist die dimensionslose Zirkulation ["* die einzige GréBe, die zur Ermittlung von Schaufel- 
kraft und Umlenkung der Strémung im Gitter bestimmt werden muB. Sie kann, wie sich spater 
noch zeigen wird, als lineare Funktion von ctg 8, angegeben werden. 


b) Induzierte Geschwindigkeiten der Einzelsingularitaten. Die gesuchte Ge- 
schwindigkeitsverteilung auf der Schaufelkontur ergibt sich, indem fiir jeden Konturpunkt die 


1 H. Kriiger, Ing.-Arch. 26 (1958), S. 242—267. 
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Summe der induzierten Geschwindigkeiten aller Singularitaten zur tiberlagerten Translations- 
strémung addiert wird. Hierzu sollen zunachst die induzierten Geschwindigkeiten der Einzel- 
singularitaten bestimmt werden. 

Gegeben sei nach Abb. 2 die Einzelsingularitat S=Q+iTl y 
(bestehend aus Quelle Q und Wirbel J’) im Nullpunkt der soup 
z-Ebene. Nach Biot-Savart induziert sie im Punkt z = x +i ¥ 
die komplexe Geschwindigkeit 


eae ee, ee PIT) Sty) 4 NI 
ey LCs i a 2 x (x? + y?) : (8) ae 
Dabei weist die Komponente u; in x-Richtung, die Komponente v; : rg 


in y-Richtung. 


Entsprechend gilt fiir die Singularitatenreihe bei der gemaB 
Abb. 2 auf der y-Achse im Abstand t eine unendliche Reihe von 
5 


gleichen Singularitaten S angeordnet ist, nach Pistolesi} 


Abb. 2. Induzierte Geschwindigkeiten einer 
Singularitatenreihe. S Singularitat bestehend 
aus Quelle Q und Wirbel I’; t Teilung (gegen- 
, : wee seitiger Abstand der unendlich vielen Singu- 
(Q +1 af) (Sin [2 EA x/t] it Sin [2 It y/t]) (9) laritaten in Richtung der Gitterfront) ; 1; in- 
= = : duzierte Geschwindigkeit mit den Komponen- 

2t (oj [2 14 x/t] —— cos [2 x y/t}) ten u; senkrecht zur Gitterfront und v; in 

Richtung der Gitterfront. 


: S Arg 
ww; = u;—t Uv; = 9,010 


Durch Einfiihrung der dimensionslosen Werte 


e Wa e_ % ee reed 
0: ees ian or 3 Q oe oie LO hee owes (10) 
und der EinfluBfunktionen 
‘ lx F ly 
Up = a (x? + y?) 2 Vr = Ga yi) = y) 5 (11a) 
tae Sin (2 2 x/t) SRC sin (2 7 y/t) 
ae Coj (2 a x/t) — cos (2 a y/t) ’ sie Cof (2 a x/t) — cos (2 2 y/t) ” (11b) 
itn = fi— fig tf, ig =F — gill (110) 
ergibt sich fiir die induzierten Geschwindigkeiten u* und v* der Einzelsingularitat aus (8): 
us = Q* tz + I* 92, ve = Q* o,—I* tz, (12 a) 
der Singularitatenreihe aus (9): 
ut =H +I, of =Q*H—T* A) I. (12b) 


Ferner gilt fiir die spater bendtigte sog. Singularitaten-Restreihe, bei der die Einzelsingularitat 
im Nullpunkt fortgelassen ist: 


ut =O" fig + T* iq) Yt, vf =(Q* bp —I™ fig) It. (12¢) 
Die EinfluBfunktionen (11) sind fir die numerische Rechnung in Abb. 3 aufgetragen. 


c) Vorgegebene linienhafte Verteilungen und Belegungslinie. Sowohl beim Einzel- 
fliigel, als auch beim Schaufelgitter sind zur Darstellung des Strémungsfeldes linienhafte Singu- 
laritatenverteilungen besser geeignet als Einzelsingularitaten. Diese Verteilungen werden 

auf einer Belegungslinie (Skelettlinie) vorgegeben. 


Bei unendlich diinnen Profilen stimmt die Skelettlinie mit der Profilkontur iiberein. Bei dicken 
Profilen wird sie so im Profilinneren gewahlt, da8 ihr Abstand von der Saugseite etwa halb so grof 
ist wie von der Druckseite. Kleine Kriimmungsradien und starke Anderung der Kriimmung werden 
nach Mdglichkeit vermieden. Als Anfangs- und Endpunkt der Skelettlinie werden die Mitte zwi- 
schen der Nase und deren Kriimmungsmittelpunkt bzw. der Hinterkante und deren Kriimmungs- 


mittelpunkt gewahlt (siehe Abb. 6 und 7). 


3 Wahrend zur Berechnung unendlich diinner Profile eine linienhafte Wirbelverteilung y(C) 
langs der Skelettlinie z(C) geniigt, ist bei dicken Profilen auBerdem noch eine linienhafte Quell- 


1 E. Pistolesi, Aerotechnica 17 (1937) und 18 (1938). 
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Senken-Verteilung q(¢) lings der Skelettlinie erforderlich. Zur Darstellung dieser linienhaften 
Singularitaétenverteilung s(¢) dient der Reihenansatz 


0) = 40) Hin) =a SOF HiT ot). (13) 


gr 


| 


J 
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Hierbei ist 1; = 1,/2 % eine Teillange der Skelettlinienlange 1,; ferner 1 = C/3 1, der dimensionslose 
Ortsparameter der ebenfalls dimensionsfreien Verteilungsfunktionen o,(z) aad % die Anzahl der 
vorgegebenen Einzelverteilungen o,(t). Diese Einzelverteilungen sollen sowohl als Quell-(Faktor Q*) 
als auch als Wirbelverteilungen (Faktor I ,) verwendet werden. AuBGerdem sollen diese Faktoren 
gemaB (10) die dimensionslosen Intensitaten der Einzelverteilungen darstellen. Es soll also. fiir 
die Intensitat der gesamten Singularitatenverteilung s(£) gelten: 

i l 


s 


{© de = fla) +in@de=2tu YO +i Ty. (14) 
ri 
Damit ergibt sich zugleich fiir die dimensionslose Gesamtzirkulation der Gleichungen (6) und (7): 
ap (15) 
R= 
Durch Einsetzen von (13) in (14) erhalt man fiir die Verteilungen o,(t) die Bedingung 
; I, 2/3 
35 J ar) de = J ont dz = be. fix alles. (16) 
AuBerdem mu die vorgegebene Quellenverteilung q(¢) die SchlieBungshedingung 
I, 
fa) dt =0 (17) 
erfiillen, die auf 
Sor =0 (17) 
ra 


hinauskommt. 

Ausgewahlt werden fiir die Singularitatenverteilungen drei Grundverteilungen (Abb. 4), deren 
besonderes Kennzeichen darin bésteht, da sich jede von ihnen nur iiber ein kurzes Stiick der 
Skelettlinie erstreckt. Die erste Verteilung 
o,(t) erstreckt sich tiber ein Stiick der 

_Skelettlinie an der Profilnase und stellt 
dort im wesentlichen die erste Birnbaum- 
Verteilung (Cotangens-Glied des Glauert- 
schen Ansatzes) dar. Die zweite Vertei- 
lung o,(t) beginnt ebenfalls an der Profil- 
nase (Anfangspunkt der Skelettlinie) und 
stellt mit ihrem ersten Stiick die zweite 
Birnbaum-Verteilung dar. Ihr Spiegelbild 
wird als Verteilung o,(t) am Ende der Ske- 
lettlinie (Hinterkante) verwendet, so daB 
im letzten Stiick der Skelettlinie ebenfalls 
die zweite Birnbaum-Verteilung dargestellt 
werden kann. Die dritte Grundverteilung 
o,(t) 148t Nase und Hinterkante unbeein- 
fluBt und kann beliebig oft (an verschie- 
denen Stellen der Skelettlinie) verwendet 
werden. 

Die Wahl dieser drei Grundverteilungen 
hat zwei Vorteile: erstens geniigen die in 
den Tabellen angegebenen Hilfswerte, um 
die Rechnung mit beliebig groBer Anzahl 
von vorgegebenen Verteilungen durchfiih- -1 -2 - oor ¢F g 
ren zu kénnen; zweitens werden mit zu- 4), 4. Vorgegebene Singularititen- (Wirbel- und Quell-)Verteilungen, 
nehmender Anzahl der vorgegebenen Hinzel- —_ + Koordinate langs der Skelettlinie (siehe Abb. 8); o(t) Verteilungsfunk- 

: Sieur tionen, von denen jede nur lings eines Stiickes der Skelettlinie definiert 
verteilungen deren Intensitaten OF und p lie ist: o,(t) und o,(r) an der Profilnase; o,(r) an der Hinterkante; o,(z) kann 
kleiner, so daB die Genauigkeit bei der Be- dazwischen beliebig oft vorgegeben werden. 
rechnung der induzierten Geschwindigkeiten 
wachst, ohne daB die EinfluBfunktionen genauer bestimmt werden miissen. Es wird dafiir in 
Kauf genommen, da® die formale Darstellung dieser Grundverteilungen etwas umstandlich ist. 
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Auf die praktische Rechnung hat das ohnehin keinen Einflu8, da die Gleichungen fiir die Grund: 
verteilungen dabei nicht benétigt werden. Diese lauten 


0,(t) = 0,735 ctg [0,5 arc cos (1 — 1,449 r)] — 0,605 o,(r) fir 0 2 (a (18) 
o,(z) = 1,138 0,799 + 0,5337,—72 fir —2/3<7,=1—2/35—0,1 
— 0,5 (1 + cos xt) fe 01S eee eee (19) 
o,(t) =0,5(1+ cosat,) verwendet fir 352 rix—1 
im Bereich —1$7,=1—~ +4541, (20) 


2% 2 
o,{z) = 0,5 (1 + cosxt,) fir’ —T1 <S%=t—-Z4+ 5S + 0,1 


2 2 
fir +01 57%,=71—-44fe4 5 CD 


= 1,138 /0,799 + 0,533 t,, — 72 


Diese vier Verteilungen sind in Abb. 4 dargestellt. 

Sowohl bei Profilen mit sehr dicker Nase, als auch bei Profilen mit dicker Hinterkante hat es 
sich fiir die Berechnung als zweckmaBig erwiesen, neben diesen Singularitatenverteilungen noch 
gewisse Zusatzsingularitaten zu verwenden. 

Sofern der Kriimmungsradius der Profilnase gréBer als 0,41; ist (Abb. 12), wird ein Dipol als 
Zusatzsingularitat zur Darstellung der Profilnase verwendet. Dargestellt wird er durch zwei gleich 
starke entgegengesetzt wirkende Einzelsingularitaten. Die eine davon wird in den Anfangspunkt 
der Skelettlinie gelegt, die andere wird so angeordnet, daf der Mittelpunkt der Nasenabrundung 
in der Mitte zwischen diesen beiden Einzelsingularitaten liegt. Die Intensitat dieses Dipoles wird 
durch die dimensionslose GréBe D angegeben. Dabei ist Dp = ['/2 1 U die Zirkulation des Wirbels 
im Anfangspunkt der Skelettlinie zur Kennzeichnung des Wirbelpaares, Dg = Q/2 1 U = Starke 
der Quelle im Anfangspunkt der Skelettlinie zur Kennzeichnung des Quell-Senken-Paares. Da diese 
beiden GréBen weder die Gesamtzirkulation noch die Gesamtquellstarke des Profiles verandern, 
brauchen sie in den vorangehenden Gleichungen nicht beriicksichtigt werden. 

Ist der Kriimmungsradius der Profilhinterkante gréBer als 0,08 1; (Abb. 11), so wird zur Dar- 
stellung des Profiles eine zusatzliche Senke verwendet, die im Endpunkt der Skelettlinie angeordnet 
wird. Ihre Intensitat wird gemaB (10) gekennzeichnet durch die dimensionslose GréBe — Q% 1. 
Sofern sie verwendet wird, muB sie in die Summen der Gleichungen (13), (14) und (17a) mit auf- 
genommen werden. 


d) Induzierte Geschwindigkeiten der vorgegebenen Singularitatenverteilungen. 
Zur Berechnung der induzierten Geschwindigkeiten u*, = u;,/U und v*. = v;,/U einer vorgegebenen 
Singularitatenverteilung o,(t) wird diese gemaB Abb. 5 in zehn gleich starke Einzelsingularitaten 


Abb. 5. Induzierte Geschwindigkeiten einer Wirbelverteilung. 
t Koordinate lings der Skelettlinie; | i Teillange der Skelettlinie; 
t 4 Koordinaten der 10 gleichstarken Wirbel, durch die die Wirbel- 
verteilung ersetzt wird; uj, und v;r Komponenten der im Punkt z 
induzierten Geschwindigkeit in x- bzw. y-Richtung; h Abstand 
des Punktes z von der Skelettlinie; y Winkel der Skelettlinien- 
tangente gegen die x-Achse am Fufpunkt des Lotes von Punkt z. 


Abb. 6. Komponenten der resultierenden Geschwindigkeit im Kon- 

turpunkt A; u senkrecht zur Gitterfront und v in Richtung der 

Gitterfront; we in Richtung der Kontur; w —p senkrecht zur Kontur; 

6 Winkel zwischen Gitterfront und Konturtangente im Punkt A; 

h Abstand des Punktes A von der Skelettlinie; ¢ Winkel zwischen 

Kontur- und Skelettlinientangente in den Schnittpunkten des Lotes 
von A auf die Skelettlinie. 


zerlegt, die im Schwerpunkt des jeweiligen Flachenstiickes der Verteilung angeordnet sind. Die 
Abszissen t der Kinzelsingularitaten (bezeichnet mit 7,) sind in Tabelle 1 fiir die einzelnen Ver- 
teilungen angegeben. Sofern der Abstand h des Punktes z (Abb. 5) gréRer ist als 0,4 1; ergibt sich 
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u;, als Summe der fiir die zehn Einzelsingularitaten gebildeten u;*-Werte von Gl. (12). Entsprechen- 
des gilt fiir v4. Dies sowie die Aufstellung der m- und n-Werte (Tabelle 2 bis 5) fiir Punkte <z, 
deren Abstand h kleiner ist als 0,4 I;, ist im Anhang II abgeleitet. 

Zur Bestimmung der Druckverteilung an der Schaufelkontur benétigt man die Geschwindig- 
keiten wx in Richtung der Kontur und zur Erfiillung der kinematischen Strémungshedingung 
nach Abschnitt 2 e) die Geschwindigkeiten wy senkrecht zur Kontur, vgl. Abb. 6. Es sollen des- 
halb die Geschwindigkeiten u* und v*. in die Komponenten wt = wx/U parallel zur Kontur 
und ws = wp/U normal zur Kontur umgerechnet werden. Dabei werden sogleich die Komponen- 
ten U und V., = — Uctg #.,, der iiberlagerten Translationsstrémung mit erfaBt: 


CS u/USH 1+ > ut, und v* = v/U = — ctg Bf, + cores (22) 
esl 


r= 


Mit dem Konturwinkel 6 nach Abb. 6 ergibt sich fiir jeden Punkt der Kontur 
wy =u* sin 6 + v* cos 0d, wh = u* cos 6 —v* sind. (23) 


Unter Beriicksichtigung von Gl. (7), (22) und (49) ergibt sich hieraus fiir das Gitter (der Einzel- 
fligel wird von hier ab getrennt in Abschn. 4 behandelt): 


we = sin 6 — ctg B, cos 6 +( 3 OF M,— STEN.) ie ie 
Ve | rt 


x ; : (24) 
wy, = cos 6 + ctg fp, sind +( 2 Qe Neale ae M,) 
n't! hae 
mit 
M, = (1 + u,) sin 6 + », cos 0, N, = (1 + 4,) cos 6 — 1, sind. (25) 


Hierbei wurde von der SchlieBungsbedingung Gl. (17a) Gebrauch gemacht.. Es kénnen daher 
mit (24) und (25) nur geschlossene Profile berechnet werden. 


Fiir die praktische Berechnung ist es zweckmaBig, (25) mit den Gleichungen (56) und (57) des’ 
Anhanges zusammenzufassen. Mit ¢ nach Abb. 6 ergibt sich 


10 10R 
M, = i + = Son tea] sin 6 + ss > v (e—2[t]) cosd+ AM, 
: Si 10,71 


(26) 
1 oR 1 oR, 
- N,=|1+— DY a (z—2[r])|cosd —— > v (ez —2[x,]) sin 6 + AN 
10,7) 10,7 
mit 
AM —AN =0 fir h/l; > 0,4, 
bzw. 
AM = on (m, cos € + n, sin €) (27) 
und fir h/l; < 0,4. 


{nee (m, sin € — n, cos €) 
31; 


R 

In (26) ist fiir die Saugseite + AM, fir die Druckseite —AM einzusetzen. Das Zeichen Ba bedeutet : 
Bei Punkten mit h/l; < 0,4 sind die Werte 4 und %, deren 2[t,] auf dem von 2 auf die Skelettlinie 
_ gefallten Lot liegt, durch die Werte fig (2 —- 2[t,]) baw. ?p (2 —2[t,]) zu ersetzen. Man beachte 
auBerdem die unterschiedliche Festlegung des Winkels ¢ fiir Saug- und Druckseite (Abb. 6). Da- 
durch geniigt es, die m- und n-Werte in den Tabellen 2 bis 5 nur fiir eine Seite zu tabellieren. Die 
Koordinaten von z sind dabei durch 7, nach (18) bis (21) und h/l; nach Abb. 6 bestimmt. Tabelle 5 
enthalt zusatzliche m- und n-Werte fiir Punkte der Profilnase nach Abb. 7. 

Zu bestimmen sind nun noch die induzierten Geschwindigkeiten der Zusatzsingularitaten, die 
bei dicker Profilnase bzw. dicker Hinterkante gema® Abschnitt 2c) verwendet werden. Fiir den 
Dipol gilt nach (12b) 


uf = {Do [ai @ — a+ Dl) — i @—A—D))] + Dr [> @ A+ DI) —4 @—A-D) fe, (28) 
vt = [Dg [* (s —2[+ D]) — > (@ —2f— D))|— Dr [# (s —z[+ D]) — 4 (@—2al— D))}} It . 


358 W. Richter: Berechnung der Druckverteilung von ebenen Schaufelgittern Ingenieur-Archiv 


Dabei kennzeichnet 2[-- D] den Anfangspunkt der Skelettlinie (Lage der einen Singularitat des 
Dipoles) und z[— D] die Lage der anderen Singularitat des Dipoles. Setzt man analog zu (25) 
Mo [ie— at DO AOA PU 
Np = [fe @ —2[+ D)) — i @ —2[— D))] cos 6—[F (@ — [+ D]) —# (@ —2[— D))] sino, 
so kénnen die induzierten Geschwindigkeiten des Dipoles in (24) mit erfaBt werden. Dabei wird 
die Gleichung fiir w um den Summanden 
und die Gleichung wf um den Summanden 

Aw}, = (Dg Np + Dr Mp) lit 
erweitert. 


Tabelle 1. Lage und Starke der Einzelsingularitaten auf der Skelettlinie, 
vgl. Abb. 4 und 5 


Anzahl der zugehérigen Werte yi (2 — 2[z,]) 
die fiir die einzelnen Verteilungen o,, o,,...in Gl. (53) 
einzusetzen sind 
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Tabelle 2. m,- und n,-Werte zur Verteilung o,(t) nach Gl. (18) 
(Hilfswerte zur Berechnung der induzierten Geschwindigkeiten fiir Konturpunkte mit h/l; < 0,4; s. Abb. 5 und 6) 
a 


m,- Werte 

Ajj; = 0 | 2,12 0,338 0,192 0,055 0,121 | —0,126 0,029 0,046 | 0,030 

0.02 0,90 0,266 0,165 0,037 0,107 0,135 0,023 
> ; ' ; — 0, : 0,043 0,029 
0,04 —0,10 0,204 0,141 0,023 0,095 | —0,140 0,018 0,041 0,028 
0,06 —0,32 0,154 0,118 0,012 0,084 | —0,142 0,014 0,039 0,026 
0,08 —0,43 0,113 0,098 0,004 0,074 | —0,141 0,010 0,037 0,025 
0,10 —0,47 0,079 0,080 | —0,002 0,065 | —0,139 0,007 0,035 0,024 
0,12 —0,50 0,055 0,065 | —0,007 0,057 | —0,135 0,004 0,033 0,02 
> > 9 9 > > ° 7 ° 3 
0,14 —0,51 0,037 0,053 | —0,011 0,050 | —0,130 0,002 0,031 0,022 
0,16 —0,50 0,022 0,043 | —0,014 0,043 | —0,124 0,000 0,029 0,021 
0,18 —0,48 0,009 0,034 | —0,017 0,037 | —0,117 | —0,002 0,027 0,020 
0,20 —0,44 — 0,003 0,027 | —0,018 0,032 | —0,109 | —0,004 | 0,025 0,020 
0,22 —0,40 —0,012 0,021 | —0,020 0,027 | —0,101 | —0,006 0,024 0,019 
0,24 —0,37 —0,020 0,016 | —0,021 0,022 | —0,094 | —0,008 0,022 0,018 
0,26 —0,33 —0,027 0,012 | —0,022 0,018 | —0,087 | —0,009 0,020 0,017 
0,28 —0,29 —0,033 0,008 —0,022 0,014 —0,080 | —0,010 0,019 0,016 
0,30 —0,25 —0,037 0,005 | —0,022 0,010 | —0,073 | —0,011 0,018 0,016 
0,32 —0,22 —0,040 0,002 —0,022 0,007 | —0,067 | —0,012 0,016 0,015 
0,34 —0,19 —0,042 0,000 | —0,022 0,004 | —0,062 | —0,012 0,015 0,015 
0,36 | —0,17 | —0,044 | —0,001 | —0,021 0,002 | —0,057 | —0,013 0,013 0,014 
0,38 —0,15 —0,045 | —0,003 | —0,021 0,000 | —0,052 | —0,013 0,012 0,014 
0,40 —0,13 —0,046 | —0,004 | —0,020 | —0,001 | —0,048 | —0,013 | 0,011 0,013 
n,-Werte 
3T= 0,05 | 040 =| 0,85 1,20 1,562. ni, 1,85 | 2,15 2,45 2,80 

h/lj = 6,54 1,918 0,984 0,625 0,376 0,216 0,106 0,039 0,004 
0,02 5,97 1,782 0,948 0,559 0,367 0,186 0,101 0,039 0,005 
0,04 Dear 1,666 0,912 0,502 0,357 0,160 0,097 0,039 0,006 
0,06 4,77 1,562 05876 0,448 0,348 0,137 0,092 0,039 0,008 
0,08 4,17 1,463 0,843 0,401 0,339 0,119 0,088 0,038 0,009 
0,10 3,67 iepyies 0,813 0,359 0,330 0,103 0,083 0,038 0,010 
0,12 3.20 L287 0,782 0,320 0,321 0,089 0,078 0,038 0,010 
0,14 2,85 1,208 0,752 0,283 0,313 0,077 0,074 0,037 0,011 
0,16 25 1,132 0,723 0,250 0,304 0,065 0,070 | ~ 0,037 0,012 
0,18 2,34 1,062 0,695 0,219 0,296 0,055 0,066 0,037 0,013 
0,20 2,15 0,993 0,667 0,190 0,287 0,046 0,062 0,036 0,014 
0,22 1,98 0,931 0,642 0,164 0,279 0,037 0,058 0,036 0,014 
0,24 1,83 0,872 0,616 0,141 0,271 9.030 0,054 0,035 0,015 
0,26 1,69 0,818 0,591 0,121 0,263 J 024 0,051 0,035 0,016 
0,28 1,58 0,769 0,567 0,103 0,255 0,018 0,047 0,034 0,016 
0,30 1,48 0,725 0,945 0,087 0,248 0,014 0,044 0,034 0,016 
0,32 1,38 0,683 0,522 0,074 0,240 0,011 0,042 0,033 0,016 
0,34 1,30 0,646 0,501 0,063 0,233 0,009 0,040 0,032 0,017 
0,36 1,23 0,611 0,482 0,053 0,225 0,007 0,038 0,032 0,017 
0,38 1,17 0,580 0,464 0,045 0,218 0,006 0,036 0,031 0,017 
0,40 1,11 0,549 0,447 0,040 0,211 0,005 0,034 0,030 0,018 - 
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Tabelle 3. m,- und ng- Werte zur Verteilung o,(t) nach Gl. (19) fiir t, < 0 und zur Verteilung o,(t) nach Gl. (21) 


fir tT, > 9 : 
Nene eee eee een eee ee ee eee eee ee earn 
Werte fiir m,(t,) = — myz(—tx) 
3 = —1,95 —1,60 —1,15 —0,80 —0,45 —0,15 
h/l; = 0 —0,141 —0,081 —0,021 —0,065 0,024 —0,010 
0,02 —0,105 —0,069 —0,015 —0,059 0,027 —0,008 
0,04 —0,080 —0,059 —0,010 —0,054 0,029 —0,007 
0,06 —0,062 —0,050 —0,006 —0,049 0,030 —0,005 
0,08 —0,047 —0,042 —0,003 —0,044 0,930 —0,004 
0,10 —0,034 —0,035 0,000 —0,040 0,029 —0,003 
0,12 —0,022 —0,029 0,002 —0,036 0,028 —0,002 
0,14 —0,012 —0,023 0,004 —0,032 0,026 —0,001 
0,16 —0,004 —0,018 0,005 —0,028 0,024 0,000 
0,18 0,002 —0,013 0,006 —0,024 0,022 0,000 
0,20 0,008 —0,009 0,006 —0,021 0,019 0,001 
0,22 0,012 —0,005 0,007 —0,017 0,016 0,001 
0,24 0,016 —0,002 0,007 —0,015 0,014 0,002 
0,26 0,018 0,001 0,007 —0,013 0,012 0,002 
0,28 0,020 0,004 0,006 —0,011 0,010 0,002 
0,30 0,022 0,006 0,006 —0,009 0,009 0,002 
0,32 0,023 0,008 0,005 —0,008 0,008 0,002 
0,34 0,024 0,009 0,005 —0,006 0,007 0,002 
0,36 0,024 0,010 0,004 —0,005 0,006 0,001 
0,38 0,023 0,011 0,003 —0,005 0,006 0,001 
0,40 0,022 0,011 0,003 —0,004 0,006 0,001 


Werte fiir n,( 7.) = ny,(—tx) 


ar —1,95 —1,60 1s —0,80 | 0,45 | 0,15 
h/lj = 0 0,199 0,547 0,761 0,870 0,949 0,994 
0,02 0,181 0,528 0,723 0,823 0,888 0,925 
0,04 0,163 0,509 0,684 0,776 0,832 0,859 
0,06 0,145 0,489 0,648 0,730 0,777 0,797 
0,08 0,129 0,467 0,612 0,686 0,724 0,737 
0,10 0,114 0,446 0,576 0,642 0,672 0,680 
0,12 0,099 0,425 0,540 0,598 0,621 0,625 
0,14 0,085 0,404 0,506 0,554 0,570 0,571 
0,16 0,073 0,383 0,474 0,513 0,522 0,522 
0,18 0,062 0,363 0,443 0,475 0,477 0,476 
0,20 0,052 0,344 0,412 0,437 0,437 0,433 
0,22 0,045 0,325 0,383 0,402 0,402 0,394 
0,24 0,039 6,307 0,356 0,373 0.373 0,363 
0,26 0,033 0,290 0,332 0,348 0,347 0,336 
0,28 0,028 0,274 0,312 0,326 0,325 0,312 
0,30 0,024 0,260 0,295 0,307 0,305 0,294 
0,32 0,021 0,248 0,278 0,290 0 
’ > : 5 ,288 0,277 
0,34 0,017 0,237 0,265 0,275 0,273 ; 
Ae é 5 0,262 
; 0,015 0,228 0,253 
; 0,262 0,260 0,249 
0,38 0,013 0,219 0,241 0,250 0,248 0,237 
0,40 0,011 0,211 0,231 0,239 0,237 0,226 
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Tabelle 4. m,- und n,-Werte zur Verteilung 0,(t) nach Gl. (20) sowie zur Verteilung o,(t) nach Gl. (29) furz,->0 
und zur Verteilung o,(t) nach Gl. (21) fiir t,, <0 


Eee 


Werte fiir m,(z;) m,(— t;) = m,(r.) my(— tx) 
37, = 0,15 0,45 0,80 1,20 Pomeen 185° | 2.15 2,45 2,80 
hij =0 | 0,011 | —0,017 | 0,037 | 0,036 | 0,211 | 0,113 | —0,137 | 0,017 | 0,028 
0,02 0,009 | —0,024 | 0,028 | 0,027 | 0,199 | 0,102 | —0,142 | 0,011 | 0,023 
0,04 0,007 | —0,027 | 0,021 | 0,021 | 0,186 | 0,092 | —0.145 | 0,005 | 0020 
0,06 0,005 | —0,030 | 0,015 | 0,016 | 0,174 0,082 | —0.146 | 0,000 0/018 
0,08 0,003 | —0,031 | 0,010 | 0,012 | 0,162 | 0,073 | —0,145 | —0,004 | 0,016 
| | 
0,10 0,002 | —0,032 | 0,006 | 0,008 | 0,150 | 0,065 | —0,141 | —0,008 | 0,014 
0,12 0,000 | —0,032 | 0,003 | 0,005 | 0,139 0,058 | —0,135 | —0,012 | 0,012 
0,14 | —0,001 | —0,031 | 0,000 | 0,003 | 0,127 | 0,051 | —0128 | —0,015 | 0.010 
0,16 | —0,002 | —0,030 | —0,002 | 0,001 | 0,116 | 0046 | —0.121 | —0.018 | 0,008 
0,18 | —0,003 | —0,028 | —0,003 | 0,000 | 0,106 | 0,041 | —0.113 | —0.020 | 0,006 
0,20 | —0,004 | —0,026 | —0,004 | —0,001 | 0,096 | 0,036 | —0,105 | —0,022 | 0,005 
0,22 | —0,005 | —0,024 | —0,005 | —0,002 | 0,087 | 0,032 | —0,098 | —0,023 | 0,003 
0,24 | —0,005 | —0,022 | —0,005 | —0,003 | 0,079 | 0,028 | —0,090 | —0,024 | 0,002 
0,26 | —0,006 | —0,020 | —0,005 | —0,003 | 0,071 | 0,025 | —0,083 | —0,025 | 0/001 
0,28 | —0,006 | —0,018 | —0,005 | —0,003 | 0,064 | 0.022 | —0.077 | —0,026 | 0,000 
| | | 
0,30 | —0,006 | —0,017 | —0,004 | —0,003 | 0,058 | 0,019 | —0,071 | —0,026 | —0,001 
0,32 | —0,006 | —0,015 | —0,004 | —0,003 | 0,052 | 0,017 | —0.065 | —0,026 | —o,002 
0,34 | —0,006 | —0,014 | —0,004 | —0,003 | 0,047 | 0,015 | —0,059 | —0,026 | —0,002 
0,36 | —0,006 | —0,012 | —0,003 | —0,003 | 0,043 | 0,013 | —0,053 | —0,025 | —0,003 
0,38 | —0,006 | —0,011 | —0,003 | —0,003 | 0,038 | 0,012 | —0,048 | —0,024 | —0,003 
0,40 | —0,006 | —0,010 | —0,002 | —0,002 | 0,034 | 0,010 | —0,043 | —0,023 | —0,003 
Werte fiir n,(t;) = n;(—T,) = n.(t2) = ny(— Tx) 
3t; = 0,15 0,45 0,80 1,20 1,55 1,85 2,15 2,45 2,80 

h/l; = 0 0,994 0,946 0,835 | 0,655 | 0,474 0,321 | 0,185 | 0,081 | 0,011 
0,02 0,916 0,878 0,778 | 0,622 | 0,417-| 0,313 | 0,160 | 0,076 | 0,012 
0,04 0,838 0.810 0.726 | 0,590 | 0,364 | 0,304 | 0.133 | 0,072 | 0.013 
0,06 0,767 0,746 0,675 | 0,558 | 0,318 | 0,296 | 0,108 | 0,068 | 0,013 
0,08 0,701 0,688 0,628 | 0.528 | 0,277 | 0,287 | 0,087 | 0.064 | 0.014 
0,10 0,640 0,636 0,587 | 0,499 | 0,242 | 0,279 | 0,070 | 0,060 | 0,015 
0,12 0,590 0,586 0,547 | 0,473 | 0,208 | 0,271 | 0,054 | 0,056 | 0,015 
0.14 0,542 0,540 0.512 | 0,449 | 0,179 | 0.263 | 0,040 | 0,052 | 0,016 
0,16 0,499 0.497 0.477 | 0.426 |° 0,154 | 0,255 | 0,028 | - 0,049 | 0,016 
0.18 0,459 0,457 0.446 | 0,406 | 0,132 | 0.247] 0,018 | 0,045 | 0,017 
0,20 0,425 0,423 0,418 | 0,386 | 0,112 | 0,239 | 0,009 | 0,042 | 0,017 
0,22 0,396 0,394 0,392 | 0,367 | 0,095 | 0,232 | 0,002 | 0,038 | 0,017 
0.24 0,370 0,368 0.368 | 0,349 | 0,082 | 0,225 | —0,004 | 0,035 | 0,017 
0,26 0,347 0,346 0.346 | 0,333 | 0,069 | 0,218 | —0,008 | 0,032 | 0,017 
0,28 0,326 0,325 0.325 | 0,317 | 0,059 | 0,212 | —0,012 | 0,030 | —_ 0,016 
0,30 0,307 | 0,306 0,306 | 0,302 | 0,050 | 0,205 | —0,014 | 0,027 | 0,016 
0,290 0,290 0.290 | 0,288 | 0,041 | 0,199 | —0,016 | 0,025 | 0,016 

EP 0,275 0,275 0.275 | 0,275 | 0,034 | 0,193 | —0,018 | 0,022 | 0,016 
0,36 0,261 0,261 0,262 | 0.262 | 0,028 | 0,188 | —0,018 | 0,020 | 0,015 
0,38 0,247 0,247 0.249 | 0.250 | 0,023 | 0,183 | —0,019 | 0,018 | 0,015 
0,40 0,236 0,236 0,238 | 0,239 | 0,019 | 0,178 | —0,019 | 0,017 | 0,014 
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Tabelle 5. m,- und n,-Werte zur Verteilung 0,(t) nach Gl. (18) fiir die Profilnase 
und zur Verteilung 0,(t) nach Gl. (19) fiir die Profilnase (s. Abb. 7) 


Punkt m,-Werte m,-Werte 


(Abb. 7) 


0 

0 

0 

0 

0 

0 

0 

| 0 
0 1,52 0,032 0,023 0 0,026 
0 1,36 0.026 —0,018 0 0,024 
0 1,23 0.022 —0,014 0 0,022 
0 1,12 “2p 18. | ke OTL 0 0,020 
0 1,02 —0,014 | —0,009 0 0,019 
0 0,95 —0,012 | —0,007 0 0,017 
0 0,88 0,010. | - 6,005 0 0,015 
0 0,83 0.008  —0,004 0 0,014 
0 0,78 0.007 0,003 0 0,013 
0 0,74 —0,006 —0,002 0 0,011 


Analog hierzu wird fiir die zusatzliche Einzelsenke an der Hinterkante (im Endpunkt der 
Skelettlinie z[E]) gebildet 


Myyi=([1+ (s —2[E])] sin 6 + 9 (s —2[E]) cos 6, 


(31) 
Neoi =[1 +h (s —2[E))| cos 6 — ¥ (zs — 2[E]) sind. 
x x z+l1 zxz+l1 
Dann brauchen in (24) die Summen >'Q* M, und >’ Q? N, nur auf > Q* M, bzw. > QF N, er- 
rat eke | r= rf 
weitert zu werden, um die induzierten Geschwindigkeiten der Einzelsenke — Q%., mit zu erfassen. 


Damit kénnen fiir samtliche vorgegebenen Einzelsingularitaten und Singularitatenverteilungen 
die induzierten Geschwindigkeiten w und w} nach (24) fiir das Gitter bestimmt werden, wobei 
zugleich die tiberlagerte Translationsstrémung erfaBt wird. Die entsprechenden Gleichungen fiir 
die Berechnung des Einzelfliigels werden im Abschnitt 4 angegeben. 


e) Gleichungssystem zur Bestimmung der Q?- und J\*-Werte. Bei der exakten 
Lésung miiBte fiir samtliche Konturpunkte die kinematische Strémungsbedingung 


wp = 0 (32) 


erfiillt sein, welche zum Ausdruck bringt, daf langs der ganzen Kontur die Geschwindigkeits- 
komponente normal zur Kontur verschwindet. Da jedoch nur eine beschrankte Anzahl von Singu- 
laritatenverteilungen (x > 3) vorgegeben wird, kann diese Bedingung beim unendlich diinnen 
Profil (QF = 0) nur fiir x-Punkte, beim dicken Profil (Q* +: 0) fiir 2 x-Punkte erfiillt werden. Bei 
Verwendung des zusatzlichen Dipoles erhéht sich die Zahl dieser Punkte um 2, bei Verwendung 
der zusatzlichen Einzelsenke — QZ, 1 um einen weiteren Punkt. 


Ks zeigte sich nun, daf es nicht méglich ist, ein allgemein giiltiges Gesetz fiir die Verteilung 
dieser Punkte anzugeben, wenn zugleich verlangt wird, da die Ergebnisse der Berechnung bei 
beliebigen Profilformen befriedigen sollen. Dies liegt daran, daB die kinematische Strémungs- 
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bedingung w}, = 0 keine Einschrankungen beziiglich des Abstandes von Kontur und berechneter 
Konturstromlinie enthalt. So ist es méglich, daB dieser Abstand von der Profilnase bis zu Profil- 
hinterkante betrachtlich zunimmt, obwohl die Strémungsrichtungen in den ausgewahlten Kontur- 
punkten tibereinstimmen. Es wurde deshalb anstelle der einfachen kinematischen Strémungs- 
bedingung (32) die verscharfte kinematische Strémungsbedingung 


/ 


i) 


9 


wp(S) dé 


=a ae 


(32 a) 


angesetzt. Hierbei ist € die Koordinate lings der Kontur; é, und &, sind die Koordinaten der 
noch naher zu bestimmenden Konturpunkte (siehe Abb. 8). Diese Bedingung verlangt, daB die 


E, bn 5 
res Gre ~t,-Punkte 
Dae Punkt 30=0,05 % rae 
un Kriimmungsmitejounkt S 3 
ee s tintunke— =o 4 
Nase Ss g 


Punkt yg” 


Druckseite 


Abb. 7. Berechnung der induzierten Geschwindigkeiten fiir Kontur- 

punkte an der Profilnase. Punkt 0° wird bestimmt durch den 

Schnittpunkt der Skelettlinientangente im Punkt 37 = 0,05 mit 

der Nase; Punkt 45° durch den Schnittpunkt der Winkelhalbierenden 

zwischen Tangente und Normale im Punkt 31 = 0,05 mit der 

Kontur; h Abstand zwischen Punkt 3+ = 0,05 und Konturpunkt; 
¢ Winkel zwischen Skelettlinien- und Kontur-Tangente. 


Abb. 8. Kontur- und Skelettlinienpunkte. € Koordinate lings der 
Kontur (ausgehend von der Hinterkante tiber Druckseite, Nase zur 
Saugseite); &,, > &, usw. Konturpunkte zur Aufstellung des Glei- 
chungssystems; iibrige Konturpunkte (schwarze Punkte) zur Be- 
rechnung der Komponenten w’- und wp) der resultierenden Ge- 
schwindigkeit; t,-Punkte auf der Skelettlinie bezeichnen die Lage 


der Einzelsingularitaten, durch die die vorgegebenen Verteilungen 
ersetzt werden. 


 berechnete Stromlinie, die durch den Konturpunkt &, hindurchlauft, auch durch den Kontur- 


punkt &, hindurchgeht. Zugleich muB damit an einer Stelle £) innerhalb des Intervalles &s < & 
< &, die kinematische Strémungsbedingung w}, = 0 erfiillt werden. Die verscharfte kinematische 
Strémungsbedingung (32a) geht im Grenzfall £,— é% in die einfache Bedingung (32) iiber. 

Zur Bestimmung der J\*-Werte bei der Berechnung des unendlich diinnen Profiles werden 
x-Gleichungen gemaB (32a) angesetzt. Dabei wird die Kontur von Profilnase bis -hinterkante 
in x-Intervalle unterteilt, von denen das erste und letzte halb so lang gewahlt werden, wie die 


_ibrigen gleichlangen Intervalle. 


Bei der Berechnung von dicken Profilen ohne Zusatzsingularitaten (Abb. 8) wird sowohl die 
Druckseite der Kontur (von Profilnase bis zum Beginn der Hinterkantenabrundung) als auch die 


Saugseite der Kontur (von dem Konturpunkt, der durch die Senkrechte zur Skelettlinie im Punkt 


31, = 0,4 bestimmt wird, bis zum Beginn der Hinterkantenabrundung) in der oben angegebenen 
Weise in x-Intervalle unterteilt. 

Bei Verwendung des Dipoles zur Berechnung von dicken Profilnasen (Abb. 12) werden zunachst 
zwei gleichlange Intervalle (Lange = 1,2 Kriimmungsradius der Profilnase) an der Nase ab- 
gegrenzt. Das eine davon legt man symmetrisch zur Profilnase, das andere schlieBt auf der Saug- 
seite an. Die restlichen Stiicke von druck- und saugseitiger Kontur werden wieder in jeweils x- 
Intervalle unterteilt, wobei hier nur das letzte auf jeder Seite halb so lang gewahlt wird, wie alle 


ubrigen. 


Bei Verwendung der Einzelsenke fiir ein Profil mit dicker Hinterkante (Abb. 11) enden die 
zu unterteilenden Stiicke von druck- und saugseitiger Kontur nicht mehr am Beginn der Hinter- 
kantenabrundung, sondern dort wo die Senkrechte zur Skelettlinie im Punkt 37, = (2 x — 1) 
++ 0,15 (siehe Tabelle 1) die Kontur schneidet. Im iibrigen erfolgt die Unterteilung wie oben 
angegeben. 
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Betrachtet man im letzteren Fall das abgeschnittene Hinterkantenstiick als zusatzliches Inter- 
vall, so stimmt in jedem Fall die Anzahl der zu bestimmenden I*-Werte und Q7-Werte mit der 
Anzahl der abgegrenzten Intervalle iiberein!. Bis auf das erste druckseitige Intervall bei dicken 
Profilen wird fiir jedes der abgegrenzten Intervalle nun die verscharfte kinematische Strémungs- 
bedingung nach (32a) angesetzt. Unter Hinzunahme von (17a) geniigt dieses Gleichungssystem, 
um auch bei den dicken Profilen die unbekannten [\*-Werte und Q#-Werte (einschlieBlich der 
gegebenenfalls vorhandenen Werte D; und Dg) zu bestimmen. Die Summe in (17a) ist dabei 
gegebenenfalls bis x + 1 zu erstrecken. 

Unter Beriicksichtigung von (24) lauten die Gleichungen (32a) dieses Gleichungssystems fiir 
das Gitter (fiir Einzelfliigel siehe Abschnitt 4) 


cos dg + ctg B, sin dg + ( 5. ONS SD TS Ms) Hie (33) 
r= r=I1 
mit je einer Gleichung fiir ? = 1, 2,... zu den oben angegebenen Intervallen. Hierbei ist nach 
(32a): 
f (é) dé 
Fes Ge 34 
COS Og Rare (34) 
f9 


Die Berechnung von sin 69, M,» und N,» erfolgt analog. Beim unendlich diinnen Profil verschwindet 
in Gl. (33) die Summe > QF N,9, da Q7 =0 ist; bei Verwendung von Zusatzsingularitaten treten 
ra 


die entsprechenden Glieder hinzu. 

AbschlieBend mége fiir die Aufstellung des Gleichungssystems noch auf einen Sonderfall hin- 
gewiesen werden: Die Aufstellung von (33) fir das abgeschnittene Hinterkantenstiick bei Ver- 
wendung einer zusatzlichen Einzelsenke an der Hinterkante ist etwas unbequem. Es wird deshalb 
die aus (33) sich ergebende Gleichung ersetzt durch 


I, 
h, wx(h,) — hy wx(hy) = ae q(S) 4S — Qu 41. (35) 


Der linke Teil dieser Gleichung stellt die durch den Schnitt hindurchstrémende Menge dar (h, die 
Schnittlange von Skelettlinie bis Saugseite, hy die Schnittlange von Skelettlinie bis Druckseite, 
wx(h,) die saugseitige und wx(h,) die druckseitige Konturgeschwindigkeit an der Schnittstelle). 
Das Minuszeichen ergibt sich dabei aus der Definition von wx (siehe Abb. 6). Auf der rechten 
Seite der Gleichung steht die Menge, die in die Senken einstrémt, die hinter dem Schnitt [(h) liegen. 
- Mit Hilfe von (13), (21) und (24) laBt sich (35) auf die Form von (33) bringen. Dabei ergibt sich 
fiir die Konstanten anstelle von (34) 


cos dg = (A, sin d[h,] — hy sin d[h,]) 10/t , 
sin 0g = (hz cos 6[h,] — h, cos d[h,]) 10/t , 
Ns = (h, M,[h,] — hy M,[h,]) 10/t + 3,02, + 20,41, 
M, = (hy N, [hg] — h, N,[h,]) 10/t . 


(36) 


Die Indices x und x + 1 in der Gleichung fiir N,s bedeuten, daB die betreffenden additiven Kon- 
stanten nur bei N,9 baw. N,,+1,9 einzusetzen sind. Die Abhangigkeiten [h,] und [h,] bedeuten, 
daB die betreffenden Werte fiir den saugseitigen bzw. druckseitigen Konturpunkt auf der Schnitt-_ | 
linie zu nehmen sind. | 

Da die Konstanten des Gleichungssystems Gl. (33), namlich die GréBen cos dg + ctg B, sindy | 
lineare Funktionen von ctg f; sind, ergeben sich auch fiir die Q*-Werte und die /'*-Werte lineare 
Funktionen von ctg {, in der Form 


QF = OFo + QF, ctg By bzw. TS = Tye Te pete po (37) 
Damit ergibt sich fiir die dimensionslose Gesamtzirkulation 
Mo SI =Tet+ Thos, mit P= ST wd P= ST. (38) 
r= as rea 


i 1 nr te zusatzlichen Intervalle bei der Verwendung des Dipoles dienen zur Bestimmung der Werte 
Q un T: 
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f) Geschwindigkeits- und Druckverteilung langs der Kontur. Mit den Lésungen 
gemaB (37) kénnen aus (24) die Geschwindigkeitsverteilungen wx (&) und wi(é) langs der Kontur 
bestimmt werden. Die Abweichungen der Lésung, die sich durch w}(€) 4 0 ausdriicken, werden 
durch eine zusatzliche Quellverteilung lings der 


Kontur q*(é) = — 2 w3(é) korrigiert. Dazu werden f SEP ER oS co ee 
die von dieser Quellverteilung in Richtung der Kontur al aii en ee 
induzierten Geschwindigkeiten iiberschlagig be- | | wo aN | 
stimmt. Man erhalt hierfiir | 
€ qf [ ieee 
*(E a7 eee Ox 
Aw (é) = | q (é eos (€ x) dy. (39) : Zio | | T +1 
é Oe = <ehnas 

Es geniigt hierbei, nur einen kleinen Bereich « bei- PT Mies 1 
derseits des betrachteten Punktes € zu erfassen. Die -7 _f p39 — 
Konturkriimmung kann daher vernachlassigt werden. i f 


‘ : & a ; r ~E EAECE ES 

Bei der numerischen Rechnung wird (39) durch a a ee 
eine einfache Regel ersetzt: Die Maxima von Aw (é) Abb. 9. Geschwindigkeitsverteilung langs der Kontur. & Ko- 
° + 5 i a Dey Wear rs 
liegen dort, sae w3() rencmittlcrens Wert cwischen ne lange der Kontur nach Abb. 8; Cit w/U und Index 0 

2 Pie ir By = 90°; WE und wp Geschwindigkeitskomponenten nach 
zwei Extremwerten erreicht (Abb. 9). Dabei ist Aw(é) Abb. 6 aus dem ersten Lésungsschritt; Aw induzierte Ge- 
positiv, erent w(€) ni wachsenden E aint peace der zusitzlichen Quellverteilung lings der Kontur, 

e i : : ie wp = 0 bewirkt; we = We + 4uK Geschwindigkeit lings 

(Abb. 8). Die GréBe der Amplituden sind bei Aw #(é) der Kontur aus dem zweiten Lésungsschritt. 
und w3(é) gleich. 

Mit den so bestimmten Werten von Aw X(é) ergibt sich als endgiiltige Verteilung der Kontur- 
geschwindigkeit gemaB Abb. 9 


we(é) = we (E) + AwK(é) - (40) 
Ebenso wie Zirkulations- und Quellverteilung gemaB (37) lineare Funktionen von ctg f, sind, 
ist auch wx(&) in allen Punkten der Kontur eine lineare Funktion von ctg f, in der Form 
; wx(E) = wko(S) + wxa(é) ctg f; - (41) 
Es miissen deshalb bei der numerischen Rechnung die beiden Geschwindigkeitsverteilungen wx 0(&) 


und w%x,(&) nach (24), (37), (39) und (40) getrennt bestimmt werden. 
' Fur den Druck gilt nach der Bernoullischen Gleichung mit wy = we U 


p+< wk =p,+$W}. (42) 
Daraus ergibt sich fiir die Druckverteilung 
— Wye \? 0 
PSM = 1 — (pe) =1— (ub sin By? (43) 


wobei gq, = @ W;/2 den Staudruck der Zustrémgeschwindigkeit W, bedeutet. Diese Druckver- 
teilung mu8 fiir jeden Zustrémwinkel 6; nach (43) aus der entsprechenden Geschwindigkeits- 
verteilung wf berechnet werden. Die Einzelheiten der numerischen Rechnung werden im Anhang I 
an Hand eines Zahlenbeispieles erlautert. 


3. Ergebnisse und Vergleiche. In den Abb. 10, 11 und 12 sind fiir drei verschiedene Gitter, 
deren geometrische Daten in den Bildunterschriften angegeben sind, die Berechnungsergebnisse 
dargestellt. Die beiden Gitter in Abb. 10 und 11 haben Schaufelprofile mit groBer Wélbung 
(f/l=0,15 bzw. 0,20) und maBiger Dicke (d/l = 0,15 bzw. 0,17), wahrend das Gitter in Abb. 12 ein 

-Schaufelprofil mit groBer Wolbung (f// = 0,25) und groBer Dicke (d/l = 0,37) besitzt. Die be- 
rechneten Druckverteilungen und Abstrémwinkel sind soweit als méglich mit Messungen ver- 

-glichen worden, die im Institut fiir Strémungsmechanik der T. H. Braunschweig durchgefiihrt 
wurden!. Dabei ist die Ubereinstimmung zwischen der vorliegenden Rechnung und den Messungen 

_durchweg befriedigend. Bei dem Vergleich mu8 beriicksichtigt werden, da in der vorliegenden 

-Rechnung die Reibung nicht beriicksichtigt wird, so daB keine vollkommene Ubereinstimmung 
von Messung und Rechnung erwartet werden kann. 

Fir zwei der berechneten Gitter (Abb. 10 und 12) sind auBerdem die nach dem Berech- 
nungsverfahren von W. Traupel? ermittelten Druckverteilungen eingetragen. Diese wurden von 


1 Herrn Prof. Schlichting méchte ich dafiir danken, daB er mir diese MeBergebnisse zur Verfiigung stellte. 
2 Siehe FuBnote 3 von S. 351. 
9 25 
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Pt 


t/1=075; B,=60° 


Abb. 10. Berechnungs- 
und MeBergebnisse fiir 
ein Turbinengitter mit 
Schaufelprofilen von 
mittlerer Wélbung und 
Dicke; Teilungsverhilt- 
nis t/1= 0,75; Staffe- 
lungswinkel fs = 60°; 
Xs  Zahlenrechnung hierzu 
z siehe Anhang I. 


60° 50° 100° i 02 6 O6 08 400 2 OF 06 08 400 2 OF 06 8 10 


t/t =05; B,=51° 


> Aa | J" 
7 ce} 


° 


Abb. 11. Berechnungs- 
und MeBergebnisse fiir 
ein Turbinengitter mit 
Schaufelprofilen von 
starker Wélbung und- 
mittlerer Dicke; Tei- 


lungsverhdltnis t//=0,5; 
ee 3 Staffelungswinkel 
20° fs U Bs = 51°. 


60° 90° 120° 09 G2 OF G6 08 1400 G2 OF G6 G8 100 O2 OF 06 o8 40 


SS 


t/l=40; f,=50° 


Abb. 12. Berechnungs- 
und MefSergebnisse fir 
ein Turbinengitter mit 
Schaufelprofilen von 
groBer Wélbung und 
Dicke; Teilungsverhilt- 
nis t/l=1,0; Staffe- 
lungswinkel 6, == 50°. 


9 0 OF 06 0& WO a2 OF O6 08 wo 2 OF O6 8 10 


Abb. 10 bis 12. Schaufelgitter (links) mit drei Zustrémwinkeln bi; Abstrémwinkel fy in Abhangigkeit vom Zustrémwinkel Aj nach Rechnung 
und Messung; Druckverteilungen fiir die drei eingetragenen Zustrémwinkel Ay aufgetragen tiber der Profilsehne | (rechts), dabei gelten untere 
Kurve und volle Punkte fiir die Profilunterseite (Druckseite), obere Kurve und offene Punkte fir die Oberseite (Saugseite). 


Punkte = Messungen des Institutes fiir Strémungsmechanik der T.H. Braunschweig; gestrichelte Kurven = Berechnung von G. Hubert 
_nach W. Traupel; ausgezogene Kurven = eigene Berechnung. N 
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G. Hubert! tibernommen. Da sie teilweise besser mit den Messungen iibereinstimmen, als die hier 


berechneten Lésungen, wird angenommen, daB die Unterschiede Fehler des vorliegenden Ver- 
fahrens sind. 


4. Berechnung des Einzelfliigels. Das entwickelte Verfahren kann auch zur Berechnung der 
Druckverteilung des Einzelfliigels verwendet werden. Dieser Fall ergibt sich durch den Grenz- 
itbergang t—>oo. Dabei erhalt man aus (7) 8, = 8, =, und auBerdem Bo = 90° + a, wobei « 
den auf die Profilsehne bezogenen Anstellwinkel bedeutet. 


Da die Zerlegung des Auftriebes A nach (1) in die Komponenten senkrecht und parallel zur 
Gitterachse hinfallig wird, setzen wir fiir den Auftrieb 


A=cys W2bl. (44) 
Dabei ist 
c4=4I™* cosa 
mit 
i 1p 


es a 
L = 2LU ~21.W,.cosa, 


Bis zur Gleichung (23) sind im Abschnitt 2 Gitter und Einzelfliigel gemeinsam behandelt worden. 
Von da ab schien es sinnvoller, zunachst das Gitter allein zu behandeln, da durch die unterschied- 
liche Definition des Zu- und Abstrémwinkels 6, bzw. 6,, beim Gitter gegeniiber dem Anstellwinkel « 
beim Einzelfliigel und durch die unterschiedliche Bezugslange der EinfluBfunktionen /1 (mit Bezugs- 
lange t) und (1; (mit Bezugslange 1) eine starkere Trennung notwendig wird?. 
Zur Umstellung der Gleichungen (24) bis (43) auf die Berechnung des Einzelfliigels geniigt es, 

folgendes zu andern: 

ctg $6, wird ersetzt durch — tga, 

t wird ersetzt durch 1, 

f(...) wird ersetzt durch j4,(...)—1, 

v(...) wird ersetzt durch 7,(...) , 

sinf, wird ersetzt durch cosa. 


Um Unklarheiten zu vermeiden, sollen die Gleichungen (25) und (26) nochmals ausfihrlich an- 
gegeben werden. Sie lauten fiir den Hinzelfliigel 


M, = fz, sin 6 + Yg, cos 6, N, = [bp, 608 0 — Vp, sin db, (245) 
Mos = 5M [up (2 — 2[t,]) sin 6 + 3 be V_ (z —2[t,]) cos6 + AM, 
A= a= (26E) 


N, = 39 3" fiz @—2[t1)) 0086 — 75 3” bn @ —a[r)]) sind + AN. 
A= Aad 


Hierbei bedeutet das Zeichen »/’: Bei Punkten mit h/l; < 0,4 sind die Werte /iz und ¥,, deren 2[7,] 
auf dem von z auf die Skelettlinie gefallten Lot liegen, auszulassen. — Die Herleitung der Glei- 
chungen (24) bis (43) fir den Einzelfligel erfolgt analog zum Gitter und braucht deshalb nicht 
naher erldutert zu werden. 

Als Beispiel wurden hiernach die Druckverteilungen fiir das Géttinger Profil 387 berechnet. 
Sie sind in Abb. 13 aufgetragen und mit Messungen® verglichen. Dabei wurde die ebene Stroémung 
nach dem Ansatz von L. Prandtl’ auf die Strémung um einen Fligel mit endlicher Spannweite 
umgerechnet. Hiernach gilt 


pes co tes (45) 


1 Siehe FuBnote 4 von S, 351. 

2 Es sei erwahnt, daB sich der zweite Grund dadurch umgehen lieBe, da man beim Gitter anstelle der 
FinfluGfunktionen js und ¢ die Funktionen jig = ju 1/t und ¥¢ = ¢1/t verwendet. Diese brauchten bei der 
Berechnung des Einzelfliigels nur durch ji; und %g ersetzt zu werden. Eine derartige Darstellung der Glei- 
chungen fiir die Gitterrechnung hatte aber zur Folge, das bei der numerischen Rechnung sehr viele zusatz- 
liche Multiplikationen erforderlich sind, um die ben6tigten Werte jig und ?g zu bestimmen. Es diirfte deshalb 
sinnvoller sein, auf diese Méglichkeit zu verzichten. : wy 

8, Prandtl, Ergebnisse der Aerodynamischen Versuchsanstalt zu Gottingen. II. Lieferung. S. 45. Miin- 
chen 1923. 3 aaa ; ~d 

4 L. Prandtl, Ergebnisse der Aerodynamischen Versuchsanstalt zu Gottingen. I. Lieferung. S. 37. Miinchen 
1921. 


25* 


368 W. Richter: Berechnung der Druckverteilung von ebenen Schaufelgittern Ingenieur-Archiv 


Zwischen dem Auftriebsbeiwert c, fiir den Fliigel endlicher Spannweite mit elliptischer Auftriebs- 
verteilung und dem Auftriebsbeiwert cy bei ebener Stroémung besteht der Zusammenhang 


C= eg = UT 6086, = ALF — UG He - (46) 


Im letzten Ausdruck von (46) wurde cosa, durch 1 und sina, durch «,, ersetzt. Dies diirfte 
ebenso wie bei (45) zulassig sein, da beim Einzelfliigel nur kleine Werte von «,, mnteressieren. 


G2 


“Is 


0 Boos see 
SCL Stelle we pe fe ee O 02 OF 06 08 0 O02 OF O06 08 400 02 OF 06 06 Ww 


Abb. 13. Berechnungs- und Mefergebnisse fir einen Einzelfliigel. Auftriebsbeiwert cg in Abhangigkeit vom Zustrémwinkel « fiir den Ein- 
zelfliigel mit b/! = 5 (links); Druckverteilungen cp = (P —Pco)/qoo aufgetragen iiber der Profilsehne fiir drei verschied2ne Zustrémwinkel 
(obere Kurve fiir Ober-, untere Kurve fiir Unterseite des Fliigels). 


Punkte = Messungen der Aerodynamischen Versuchsanstalt zu Gottingen; Kurven = Berechnung. 


Durch Eliminieren von «,, fiir b,,/l,, =0o folgt aus (45) und (46) 


Rey BEd 
nay see 


=. (41) 
1 sae Ti 


als Auftriebsbeiwert fiir den endlich breiten Fliigel mit b// und Anstellwinkel « (im BogenmaB 
einsetzen). 

Zur Berechnung der Druckverteilungen in der Fliigelmitte mu noch der Winkel «,, der zu- 
geordneten ebenen Strémung bestimmt werden. Dafiir ergibt sich aus (46) und (47) durch Eli- 
minieren von ¢, 

Tj —a bl 


<6 = — 
re — jl 


(48) 
Mit diesem Winkel «,, wurden die Druckverteilungen von Abb. 13 berechnet. Die Ubereinstim- 
mung von Messung und Rechnung ist erstaunlich gut. 


5. Zusammenfassung. Es wird ein Verfahren zur Berechnung der Druckverteilung an beliebig 
vorgegebenen ebenen Schaufelgittern bei reibungsfreier, inkompressibler Strémung mitgeteilt. 
Inshesondere ist das vorgeschlagene Verfahren auch anwendbar fiir Gitter mit Schaufelprofilen 
von groBer Wélbung und grofer Dicke. Das neue Verfahren ist ein Naherungsverfahren, das nach 
der Singularitatenmethode arbeitet. Dabei ist ein System von acht linearen Gleichungen zu lésen. 
Das vorliegende Verfahren wurde durch mehrere Beispielrechnungen erprobt, wobei sich durchweg 
befriedigende Ubereinstimmung mit Messungen ergab. Die teils numerische, teils graphische 
Rechnung dauert etwa 50 Stunden pro Gitter. 

Das vorgeschlagene Verfahren ist eine Erweiterung des Verfahrens von H. Schlichting}, fiir 
das eine Berechnungszeit von 20 bis 25 Stunden pro Gitter angegeben wird, auf Gitter mit stark 
gewolbten dicken Profilen. 

Das graphische Berechnungsverfahren von W. Traupel® scheint noch etwas genauere Ergebnisse 


zu liefern, als das vorliegende Verfahren. Jedoch betragt bei jenem Verfahren die Berech i 
nach G. Hubert? etwa 200 Stunden pro Gitter. : St Oe eae 


1 Siehe FuBnote 2 von Seite 351. 
2 Siehe FuBnote 3 von Seite 351. 
3 Siehe FuGBnote 4 von Seite 351. 
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ees A, eee Die praktische Durchfithrung der Rechnung soll fiir das 
er von . 10 gezeigt werden. Dabei wird mit vier Singularitiit i = : 
rechnet. Die Rechnung verlauft folgendermaBen. oc ee 

Kine Schaufel des Gitters wird so groB auf Transparentpapier gezeichnet, daB die Gitterteilung t 
mit der Teilung des Diagramms der Kinflu8funktionen von Abb. 3 iibereinstimmt. Gema® Abb. 8 
werden die Skelettlinie und die t,-Punkte nach Tabelle 1 eingetragen. Sodann werden etwa 20 
einigermaBen gleichmaBig verteilte Konturpunkte eingetragen (siehe Abb. 8). Innerhalb des 
Bereiches mit hjl; < 0,4 miussen diese Konturpunkte so eingetragen werden, daB sie auf den Senk- 
rechten zur Skelettlinie liegen, die durch die t,-Punkte laufen. 


Tabelle 6. Berechnung eines Beispieles: Ermittlung der M- und N-Werte fiir den Konturpunkt A 
von Bild 6 und 8 nach Gl. (26) 


ae Spalte3 Spalte4 Spalte6 
oordinaten ft- und »-Werte Summen von Spalte 3 M- und N-Werte 
nach Abb, 6, 8, 10 aus Abb. 3 nach Tab. 1 nach Gl. (26) 
x,/l = 0,107 3 4 nach Tab. 1 eaitie aus Spalte 4, 5 
Ell = 1,845 5 js <j i =: i 
4 = 140,9° 3u | b ues 2 if M, = 1,994 - 0,632 
een 6 re etn 1 ‘SR 190s oe 
See 5 ; —0,01 y= 8 ey 
hjly = 0,24 0,85 0.997, 229015 |... 10 2 N, = — 1,994 - 0,776 
| 5 Mati th 095 + 0,010 - 0,632 
ae alte 2 155 0.98 0.03 y= 0,994 == 104k 
oordinaten 1.85 0.98 0.05 es —0,010 
? 2 a His Me 2,001 5 0,632 
m-, n-Werte 2,15 0,98 —0,07 lor) ri na eC + 0,107 - 0,776 
pa 53° 2,45 0,975 |-—0,10 ars 1365 = 1,348 
nach Abb. 6 2,80 O08 e OTe 8 a oas N, = —2,001 - 0,776 
cosé = 0,995 3,20 1,00 —0,22 3 = 0.269 + 0,107 2 0,632 
sine = — 0,092 3,55 1,06 | —0,31 ee eee = — 1485 
3 t(z) = 5,55 3,85 Ie eA Te |e ee ee 
nach Abb. 8 4,15 1,29 | —0,55 Spalte5 M; = 2,365 - 0,632 
p= 55> 5 ee rine Costes AM- und AN-Werte =O NAA OHO 
‘nach GI. (18) es oS | ee nach Gl. (27) a BAB a: 
oa ee 5,55 0,03*) | —0,03*) | AM, = 1,873 (0,079 - N, = — 2,365 + 0,776 
nach Gl. (19) Bo 087 4,60 - 0,995 — 0,082 - 0,092) — 0,042 - 0,632 
Bite tb 3 6,15 —1,08 2,45 = 0,133 — 0,167 
nach Gl. (20) 6,45 | —0,96 1,52 | AN, = 1,873 (—0,079 - = — 2,029 
mz = 0,079 6,80. v==0,82.|--. 0,99 - 0,092 — 0,082 - 0,995) 
ns = 0,082 Pls 30,18? 90068 = 0,107 per Ae eam 
nach Tab. 4 F008 i Of OAD feaar = 1,875 (0,022 00030 
'3t, = —0,45 >—3 - 0,995 — 0,368 - 0,092) = 0,135 
nach Gl. (21) = —- 0,023 N, = — 1,269 - 0,776 
m, = 0,022 *) iin baw. dp AN, = 1,873 (—0,022 - — 0,830 - 0,632 
n, = 0,368 da hierbei 3 t, - 0,092 — 0,368 - 0,995) — 0,689 
nach Tab, 4 == 3) 1 (Zz) 1st = — 0,689 = — 2,198 


Fiir alle diese Konturpunkte werden nun die EinfluBfunktionen M und N nach (26) bestimmt. 
Die dazu erforderliche Zahlenrechnung ist in Tabelle 6 fiir den Konturpunkt 4 von Abb. 6 und 8 
ausfiihrlich dargestellt. Fir alle ibrigen Konturpunkte erfolgt sie in genau derselben Weise: 

Spalte 1 von Tabelle 6 enthalt die Koordinaten des Konturpunktes, die aus der Aufzeichnung 
gemaB Abb. 6, 8 und 10 entnommen werden kénnen. Spalte 2 und 5 werden nur bei solchen Kontur- 

punkten bendtigt, deren h/l; < 0,4 ist. Dabei ist 3 r(z) in Spalte 2 der FuBpunkt der Senkrechten 
- gur Skelettlinie, auf der der Konturpunkt liegt. Die m- und n-Werte werden nur fir die Vertei- 
lungen o(r) bestimmt, innerhalb deren Definitionsbereich der Punkt 3 7(z) liegt. Das trifft im 
vorliegenden Fall nur fiir die Verteilungen o,(r) und o,(t) zu, da sowohl 3 t > 3, als auch 3 t > 3 
ist. Die 1M- und AN-Werte werden in Spalte 5 aus den Werten von Spalte 2 nach (27) berechnet. 

Spalte 3 von Tabelle 6 enthalt die fiir die ,-Punkte der Skelettlinie aus Abb. 3 entnommenen 
fi- und ¥-Werte. Hierzu wurde das aufgezeichnete Profil so auf die beiden Diagramme aufgelegt, 
daB sich der Konturpunkt A mit dem jeweiligen Diagramm-Mittelpunkt deckt, die Gitterrichtung 
zur Mittelsenkrechten des Diagramms parallel ist und das Profil von rechts nach links umstrémt 
wird (d. h. gegeniiber den Abb. dieses Aufsatzes um 180° gedreht ist). In dieser Lage des Profiles 
kénnen die ji- baw. ¥-Werte der Spalte 3 in den entsprechenden 37,-Punkten der Skelettlinie dem 
_ jeweiligen Diagramm entnommen werden. 
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Spalte 4 enthalt die gem4B Tabelle 1 gebildeten 
Summen der Werte von Spalte 3. Spalte 6 zeigt BPR i teat 1 ic) — 
schlieBlich die Berechnung der M- und N-Werte FUN es a; 
RBADBDBBDBD op 
aus den Werten von Spalte 4 und 5 nach (26). SHAH HHS 3 
Nach Durchfiihrung dieser Rechnung fiir die ane 28 
. . an 
anderen Konturpunkte werden die Einflu8funk- = Sas 
: : : é —— VY eS 
tionen M und_N, sowie die Funktionen 2 cos 6 Vanridadl & oe % | ep 
. “* ° 
und 2 sin 6 gemaB Abb. 14 iiber & aufgetragen. ea evouowowsdo~ Ba, oa eae 
ap ; : ; NARA 22 > wee 
Die in Abb. 8 nach Abschnitt 2e eingezeichneten PRO One F oes 
> Baa Sens ye FADS 
5 HSSSSS HO] 8 2 OBE+ 
s [APSE ARP) gee ees 
x ln i eg 
o = . 
=~ ao OG 
: weeeeee: | & 28S 
e 9 Ciers 
3148 nosoMNMNSe fl Qua s 
= tHor too Saal > 
BE * | Bae ee 
o 165 SB) 
m3 [le dest etre 
nN |S 
~ 8 onone ons 
NASASNSONSANSASN 
; ; ant NAR Ky 
-——Druckseite——~ Nase ~ Saugseite ———> 3 : concoco 
é ] = — | SDNnAMrKons 
3 a OE oe 
foe | BE se ei ced 
NS cos ~ = | fae ee al 
; See Ros go | BAAS | @ 
- — = pa © BD BD BD BD 
eee aS ela [| Snennee | SS sede 
V | De NS 3 |F =| sageese eee 
2 oo Pl sre] SSotang pe a 
N, Ww INLG|I E oo = ARAN AOS Oe ee 
te oe ee ! Ms hl g|c | | | ++++ be cep 
é or 2 E-€, a & £-E, E-f. EW EE, ‘a oS oF | a SS ie = es 
Abbi Th Auttrsguayider BGéleGtesbtionen’ If aud N ibe der‘ab- 2 = SSssssss SSass| S 
gewickelten Santee zum Beispiel von Abb. 10; # = 1; 2; = Natervallc ic) 7 - So ee ee Sic oe = 
(Abb. 8) zur Aufstellung der in tabell rete des Gleichungssystems <a (os = at it S S Ean 4 | 
: he eSeooona He wt ty tt 
«|S res Sete Sees] 
Py fax} _—_ — 
4 ° <e = ios . Ka¥ako%emy * 
Konturpunkte ¢, und & werden in-Abb. 14 iber- 3 | De = Be ry mae Ue een 
tragen. Sodann werden die Koeffizienten des 3 | ® . eS Se SS es. ae 
Gleichungssystems (33) nach (34) aus Abb. 14 5 8 as ASSARSss $) ESR a le 
entnommen (Mittelwerte der aufgetragenen Funk- a ie aT 7 T Fi a f Fa 
tionen in den einzelnen Intervallen). Das so 8 oS = 
aufgestellte Gleichungssystem ist in Tabelle 7 an- a] 8 BULAe etcebeaetar et 
gegeben. Die letzte Gleichung ist die SchlieBungs- ‘5 |* Se = 
bedingung nach (17a). Zur Vereinfachung der A & 3 tee 
numerischen Berechnung ist beim Gitter der  & | 8 rie 5 Pelee 2 
Faktor l/t von (33) in die Unbekannten J° = ]"*1/t < 3 Sp 
° sats — _a- 
und Q° = Q* I/t einbezogen. Die in Tabelle 7 eben- 5 S Stenenenenenonen el wees 
falls angegebenen Lésungen Qf,...,Q*;['*,...I°3 = (8 Se eS ae ee 
: ; : 2 = 
sind lineare Funktionen von ctg 6,. Man erhalt a laa = SaRSRSLS = S S S 
sie am einfachsten durch Eliminieren der Un- a COTES T ef zr 7 SSSS 
bekannten mit dem gréBten (zur Vermeidung = sane hy Hall 
von Rundungsfehlern) Koeffizienten. Aus der Ge- 3 HO ee Bear 
samtzirkulation nach (15) ergeben sich die Kraft- é ‘apiece Coe Soos|e 
beiwerte nach (6) und der Abstrémwinkel nach (7). 2 SSsssss Lecsed 
In Tabelle 8 werden fiir den Konturpunkt 4 5 = S ba = S 3 = Hi ey 
(Abb. 6) die Geschwindigkeiten wé und w% nach Gegsssdo = eee 
(24) berechnet. Dabei sind die Anteile 0 und £ SB loonans SSSS | o> 
gema (41) getrennt. Diese Rechnung erfolgt a Hie te) ol 
wieder fiir alle iibrigen Konturpunkte in gleicher CSAMAOAIA SSS | & 
Weise. Die Ergebnisse fiir den Anteil 0 (8, = 90°) S| SRAaSke 
sind in Abb. 9 iiber & aufgetragen. Ebenso wie aces if 7 iy ie 
* s . 
dort Awko bestimmt wurde, wird auBerdem aus 
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Tabelle 8. Berechnung eines Beispieles (Fortsetzung): Ermittlung von Geschwindigkeit und Druck 
fiir den Konturpunkt A von Bild 6 und 8 nach Gl. (24, 40, 41 und 43) 


Spalte 1 Spalte 2 Spalte 3 Spalte 4 Spalte 5 
Werte nach Werte nach Tabelle 7 

Tabelle 6 Anordnung gema8 Gl. (24) 

sin 6 = 0,632 1,000 0 0 1 
cos 6 = — 0,776 0 1,000 — 1,000 vate 

M,= 1,268 %= 0,151 | ©,=—0,028 O,=. 0013 | _re,=— 0,396 
M,= 1,348 ®=—0,075| I,= 0,366 Q,=—0,003 | Te,= — 0,076 
M,= 1,595 ®%= 0,095; T= 0,145 Q,=—0,006 | Ie, — 0,018 
M,= 0,135 8 = — 0,171. | Fees. 0,198 2,——0,004 | I'g, = — 0,007 
N, = — 1,541 ee OL025 Pee OsLou cee iy = 0,396 OF =. 0,005 
Nz = — 1,485 — I$, = — 0,366 = — —To,= tae 

2 ; 9% 2, = — 0,075 T?,= 0,076 | @$,=— 0,003 
Ng==— 2,029 | —I%,=—0,145 | Q%—= 0,095 |—r,= 0,018 | 8, = — 0,006 
N,= — 2,198 —Ir%, = — 0,198 eo = — 0,171 ee ay pee 0,001 Qt = — 0,004 
nach Gl. (24): ewe €r 5 tae oe 
oe eee "Ko : 2,080 | why = 0,001 WK = 0,003 wig = 0,003 
der Spalten und 2 lund 3 1 und 4 lund 5 
nach Abb. 9: Aw¢ = — 0,060 analog dazu: Awk g = — 0,002 
nach Gl. (40) und (41) w= 2,020 + 0,001 ctg By (fiir Punkt A) 
nach Gl. (43) fiir Br = 70 90 110° wi AbBATO 

(p= pr)/ar = — 2,60°| -- =3,08 —2,60 bei x./l = 0,7 


der Auftragung von wp, iiber & das Awkg bestimmt. AuBer der Geschwindigkeit wg nach (41) 

zeigt Tabelle 8 schlieBlich noch die Berechnung der Druckverteilung (p — p,)/q, nach (43) fiir die 
drei Zustrémwinkel 6, von Abb. 10. Auch diese Rechnung ist fiir alle itbrigen Konturpunkte in 
gleicher Weise durchzufihren. 


Anhang II. Berechnung der induzierten Geschwindigkeiten der vorgegebenen 
Singularitatenverteilungen. Analog zu (12) wird fiir die induzierten Geschwindigkeiten 
angesetzt: fiir den Einzelfliigel 


uj, = QF Mer SF ie Ver vi, = QF Ver LL? Fer: (49a) 
fiir das Gitter 
: ui, a (Qr LU, f vies V;) l/t > Vir cas (QF es yids TP) I/t 2 (49b) 
fur das Restgitter 
| uy, = (OF te, FI yg) ts vt = Oh vr, — IF eax) Yt. (49c) 
Dabei ist 
2/3 2 2/3 : < 
ie = f onle) fee ate) de, 9, = J 04(e) §(¢— ale) de (50) 
0 


- Analog dazu werden jig, aus [lz und yg, aus Vz, sowie Wp, AUS fiz und YR, aus Vp berechnet. 
Fiir alle Funktionen o,(€) ist 


I(t) = fo,(t) de (51) 
i) 
- monoton steigend von 0 bis 1. Man kann deshalb anstelle von (50) schreiben 


be =f (este) al (0). (52) 


Diese Gleichung (52) und die folgenden Gleichungen (53) und (54) gelten fiir Wg, und [lp,, sowie 
_ fir y,, vg, und vz, analog. 
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Bei der numerischen Rechnung wird das Integral der Gleichung (52) durch eine Summe er- 
setzt. Hierzu wird I,(r) in Z gleichlange Sticke unterteilt und 2(z) jeweils fiir deren Mitte bestimmt. 


Es ergibt sich 


4 ee te 
Mis Zia (2 —2[t,]) (53) 
mit T, aus 
7) 
[ ee) de = “Fz fir, PSUS ee Ze (54) 


Fiir Punkte z, deren Abstand von der Skelettlinie gréBer ist als 0,4 1;, geniigt erfahrungsgemaB 
eine Summe mit Z — 10. Fiir diesen Fall sind die erforderlichen 1,-Werte nach (54) fiir jede der 
ausgewahlten Verteilungen berechnet und in Tabelle 1 angegeben. Fiir die Verteilung o,(7) ist 
dabei zu beachten, daB am Ort der tabellierten 1,-Werte z.T. mehrere Einzelsingularitaten zu- 
sammengefaBt wurden, so daB die zugehérigen Werte ji (z —2[t,]) gemaB Tabelle 1 mehrfach 
in (53) einzusetzen sind. Diese Zusammenfassung wurde gewahlt, um trotz der Unendlichkeits- 
stelle der Funktion o,(t) mit méglichst gleichmaBig verteilten Singularitaten auszukommen. 

Fir Punkte z, deren Abstand von der Skelettlinie kleiner ist als 0,41;, miiBte zur Berechnung 
von ju, nach (53) Z entsprechend erhéht werden. Da das fiir die numerische Rechnung zu mithsam 
ware, sollen Hilfswerte tabelliert werden, mit denen die v- und y-Werte fiir Punkte z in diesem 
Bereich einfacher bestimmt werden kénnen. ErfaBt wurden dabei — ausgenommen im Bereich 
der Profilnase — nur solche Punkte z, die sich auf den durch die 7,-Punkte laufenden Senkrechten 
zur Skelettlinie befinden. In der Umgebung der Profilnase wurden die Hilfswerte fiir Punkte 
bestimmt, die auf der Tangente an die Skelettlinie im Punkt 3 7 = 0,05 und auf den durch diesen 
Punkt verlaufenden Strahlen liegen, die um 45° gegen die Tangente geneigt sind (siehe Abb. 7). 

Zu jeder der Grundverteilungen nach (18) bis (21) wurden fiir deren Anordnung auf der x- 
Achse die Werte gz, und vg, mit hinreichend groBem Z berechnet. Diese Werte sind dem Ver- 
haltnis 1/1; direkt proportional, da die in (53) einzusetzenden Werte “ nach (lla) proportional | 
und umgekehrt proportional zu x — x[1,] = 3 1; (x —t,) sind, gemaf der Definition von 7 in (13). 
Da zwischen den Langen / und /; keine direkte Beziehung besteht, sollen die Hilfswerte in folgender 
Form gebildet werden: 

31; ian oe 10, 
n= ag MEr ~ 10, bez (2 — x[1,])| n, = a VEY i Ve_ (z—-x[t,])}. (55) 


Dabei bedeutet 3”, das der Summand mit x(r,) = x(z) ausgelassen wird, sofern er in der Summe 
vorkommt. Auf Grund dieser Festlegung sind die m- und n-Werte gewissermaBen Korrektur- 
werte, die die Anwendung von (53) auch dann noch erméglichen, wenn der Abstand des Punktes z 
von der Skelettlinie kleiner ist als 0,4 1;. Die Fehler, die durch Vernachlassigung der Skelettlinien- 
kriimmung entstehen, bleiben dadurch hinreichend klein. In Abhangigkeit von h/l, (Abb. 6) 
und 7, sind diese m- und n-Werte in den Tabellen 2 bis 5 fiir die vorgegebenen Verteilungen (r = 1; 
2;...%) angegeben. Dabei ist 7, nach (18) bis (21) durch 7 bestimmt und 1 = x(z)/3 I; nach (13). 

Fiir eine Einzelverteilung auf der Skelettlinie z(r) ergibt sich mit diesen m- und n-Werten 


und y nach Abb. 5 


1 105 3 l 
Mer = 10 =, Ly (2 —2[t,]) + 31 cosy — n, sin y) , 
ee oes I ; (oe) 
ter = 75,51 be @—aEr)) — 54 (m,siny + n, cosy). 


F ir das Restgitter kénnen die Werte wp, und yp, auch fir Punkte z mit h < 0,4 I; nach (53) be- 
stimmt werden, da der Abstand von den Belegungslinien des Restgitters geniigend groB bleibt. 
So ergibt sich fiir die Werte u, und 7, des Gitters gemaB (11) ce 


Hy = Mr, + Me, t/l, Y =VR, + Vg, t/l. 7 (57) 
Hiermit wurden aus (25) die Gleichungen (26) und (27) aufgestellt. 


(Eingegangen am 7. Dezember 1959.) 
Anschrift des Verfassers: Wolfgang Richter, c/o CERN-PS, Genf 23. 
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